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PREFACE 

STIMT7LATBD by a course of jjoat-graduate lectures on the Partial 
Differential Equations of Mathematical Physics which Professor 
E. T. Whittaker gave sixteen years ago in the Mathematical Institute 
of Edinburgh University, one*of the authors of the present work 
(P. B. B.) planned a comprehensive treatise covering the whole of 
this field. Unfortunately, ill health and pressure of other duties have, 
so far, prevented the completion of this scheme. In the meantime 
the subject has been treated from different points of view by Bate- 
man (1932), Courant and Hilbert (1924 and 1938), and Webster 
(1927). 

During the same period there have been.^eat advances in mathe- 
matical physics, especially in the various developments of quantum 
mechanics. As these new theories are still developing rapidly, it 
would perhaps be unwise to attempt at the present juncture another 
general treatise on the mathematics of physics: and, after much 
consideration, we have decided to abandon the original plan and to 
replace it by. the pubhcation of a number of monographs, each 
complete in itself, on various special topics not adequately treated in 
existing books. 

The present monograph deals with the mathematical theory of 
Huygens’ principle in optics and its application to the theory of 
diffraction. No attempt is made to give a complete account of the 
various methods of solving special diffragtion problems. We are 
concerned only with the general theory of the solution of the 
partial differential equations governing the propagation of light 
and we discuss some of the simpler diffraction problems merely 
as illustrative examples. Eor an account of the more technical 
developments of the theory of diffraction we refer the reader to 
the excellent articles by von Laue and Epstein in the Encyklopadie 
der mathematischen Wissensdiaften (Band V, 3. Teil) and that of 
Wolfsohn in the Handbuch der Physik (Band XX — ‘Licht als 
W ellenbewegung ’ ) . 

The standard of knowledge of pure and applied mathematics ex- 
pected of the reader is roughly that of the undergraduate who has 
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completed the compulsory ])art of an honours course and is about to 
take up some specialized study. 

We wish to express here our great indebtedness to Professor 
Whittaker for the original stimulus Avhich led us to this work and 
for his continued interest, encouragement, and advice. We also 
desire to thank the Delegates of the Clarendon Press for undertaking 
this book and their Staff for their unfailing skill in printing it. 

B. B. B. 

E. T. C. 
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THE ANALYTICAL REPRESENTATION OF 
HUYGENS’ PRINCIPLE 

§ 1 . The principle of Huygens 

§1.1. Huygens’ geometrical theory of wave-propagation in optics 
In his Traiti de la Immilre, published in 1690, Huygens discussed 
the process of the propagation of light by the aid of a new principle, 
which has since been generally known by his name. At that timef 
light was regarded as a disturbance in a medium, the aether, of much 
the same general character as sound in air. Huygens supposed that, 
at any instant f — <q, a point-source of light generates a disturbance 
which is propagated into the surrounding medium as an isolated 
spherical wavej which expands with a large constant velocity, the 
velocity of light. This initial disturbance develops into the distur- 
bance at the subsequent instant t = through a succession of states 
at the intermediate instants, and a knowledge of the state at any 
intermediate instant t ~ t' suffices to determine the state at the 
instant t ~ Thus, if we regard each element of the isolated wave 
at the instant i — t' as the centre of a new disturbance, the actual 
effect at the instant f is the resultant of all these secondary 
effects, the actual wave the envelope of the secondary waves. 

The principle, stated in this form, is somewhat vague, and so, 
before we go on to discuss the difficulties of Huygens’ geometrical 
theory of wave-propagation in optics, it is convenient to follow 
Hadamard|l and analyse the principle in the form of a syllogism. 

A. (Major premiss.) In order to determine the effect at the instant 
t = ti of a luminous phenomenon caused by a given disturbance at 
the initial instant t = we may calculate the state at some inter- 
mediate instant t — t' and from that deduce the state at the instant 
t = <1. 

B. (Minor premiss.) If at the instant t — t^, or, more precisely 

t For tho history of the dovolopment of the theories of light, see E. T. Whittaker, 
History of the Theories of Aether and Electricity (Dublin, 1910). 

i An isolated spherical wave, or pulse as it is sometimes called, is a disturbance 
of tho medium localized on the surface of a sphere. 

II See Hadamard’s lecture ‘Le Principe de Huygens’, Bull, de la Soc. rnath, de 
France, 52 (1924), 610-40, or his book Lectures on Cauchy's Problem (Yale, 1923), 
53-6. 
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2 THE ANALYTICAL KEPRESENTATION OF [chap, i 

in the short interval ^ 0 “^ ^ ^ we produce a luminous distur- 
bance localized in the immediate neighbourhood of a point 0, the 
effect at the subsequent instant t = V is localized in a very thin 
spherical shell with centre 0 and radius c(«'— where c is the 
velocity of light. 

yC, (Conclusion.) In order to calculate the effect at the instant 
t = ti due to a luminous disturbance localized at 0 at the instant 
t = Iq, we may replace the initial disturbance by a suitable system 
of luminous disturbances taking place at the intermediate instant 
t = t' and distributed over the surface of the sphere with centre 0 
and radius 

When analysed in this way, Huygens’ principle is seen to involve 
three propositions, and different authors have attached the name 
‘Huygens’ principle’ indiscriminately to any one of the three. In 
the present book we are concerned with proposition G and its 
generalizations involving luminous disturbances which cannot be 
generated by the superposition of spherical waves. 

Proposition A would probably be accepted as immediately obvious; 
it is nothing other than the principle of determinism which runs all 
through classical mathematical physics. Nevertheless, the fact that 
the solutions of the differential equations which govern the propaga- 
tion of light do satisfy this proposition is of considerable interest, for 
it leads to certain remarkable relations connecting these solutions. f 

Proposition A is what the philosophers would describe as one of 
the laws of thought — ^its co*ntrary is inconceivable. Proposition C 
is a physical law capable of very wide generalization. Proposition B, 
on the other hand, has a very special character, since it is a property 
peculiar to certain special types of luminous phenomena. It states 
that an isolated spherical light wave has clean-cut propagation; for 
if such a wave is due to a disturbance localized at the origin and 
acting only during the very short interval tQ—e < ^ effect 

at a distance cT is null until the instant t = tQ-\-T—e and is null 
again after the instant t = ^o+ T. Thus, according to proposition B, 
an isolated spherical wave leaves no residual after-effect. We cannot 
generalize proposition B to cover, for example, two-dimensional 
wave-motions; for, in two dimensions, an initial disturbance always 
gives rise to a residual after-effect. 

t Hadamard, BuU. de la Soc, math, de France, 52 (1924), 241-78; Acta Math. 49 
(1926), 203-44; Journal de Math. 8 (1920), 107-228. 
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To sum up, although the premisses A and B do imply the con- 
clusion C, an argument of this type is incapable of generalization 
owing to the great restrictions under which the proposition B holds. 
We shall prove a general form of proposition (7, and then show that 
JS is a consequence of C for the special type of spherical wave 
considered by Huygens. 

Throughout the book we are concerned only with the generaliza- 
tions of G which are, in effect, governed by the partial differential 
equation of wave-motions. For the still wider generalizations con- 
cerning other differential equations we refer the reader to Hada- 
mard’s Yale lectures already cited. 

§ 1.2. The difficulties of Huygens* theory 

In applying his geometrical theory of wave-propagation Huygens 
encountered certain difficulties which he was able to overcome only 
by making special ad hoc hypotheses. In the first place, he found 
that he could account for the rectilinear propagation of light only by 
assuming that a secondary wave has no effect except at the point 
where it touches its envelope. Secondly, the envelope of the secondary 
spherical waves consists of two sheets, one on each side of the surface 
on which the secondary sources of disturbance lie. It would seem, 
therefore, that one of Huygens’ isolated waves would be propagated 
not only forwards but also backwards. To get over this difficulty 
Huygens had to assume that only one sheet of the envelope is to 
be considered. 

If we wish to avoid making this last assumxition, we must give uj) 
the purely geometrical theory- and have recourse to analysis. To 
illustrate this, let us consider plane waves of sound of small ampli- 
tude. We shall prove in § 3.2 that an initial disturbance is actually 
proj^agated in both directions unless certain conditions are satisfied. 
Only when the initial values and Sq of the velocity and condensation 
are connected by one or other of the relations ^ 

a plane wave which is propagated in a definite direction.! 

A similar conclusion holds for electromagnetic waves in a vacuum: 
in a jirogressive plane electromagnetic wave-motion whose wave 
fronts are perpendicular to the vector n and move in the direction 

t See a letter from Fresnel to Poisson (QJuvres compUtes de Fresnel, 2, 227 : 
quoted by Poincar6, Theorie math, de la LumUre, 1 (1889), 81; Croze, Annales de 
Physique, 5 (1926), 371-439 (380-1)) on this point. For a similar result concerning 
general sound waves, see Love, Proc, London Math, Soc, (2), 1 (1903), 37-62 (54). 
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of n, the electric and magnetic vectors d and h must be equal in 
magnitude and n, d, and h must form a set of riglit-handed orthogonal 
axes. 

§ 2. Huygens’ construction as a contact-transformation 
§2.1. The definition of a contact-transformation 

A surface element at a point P of space is specified by the coordi- 
nates (a:, y, z) of P and the direction cosines (/, m, n) of the normal 
to the element. Let a transformation from the set of variables 
(x,y,z,l,m,n) to the set (x\y\z\V be regarded as turning 
this surface element into another surface element at P' (x\y\z'), 
the normal to which has direction cosines There are oo^ 

surface elements through any given point P; to each corresponds 
in general a surface element through a different point P', since the 
coordinates (x\y\ z') depend on (Z, m, n) as well as on the coordinates 
of P. 

From the equations which define {x\y\z' ,V ,ni\n^) in terms of 
{x, y, z, Z, m, n) it may be possible to eliminate completely the direction 
cosines so as to obtain one or more relations betw^een the coordinates 
of P and P'. 

There are three cases to be considered: 

{a) There may be only a single relation 

^l[x,y,z,x\y\z') = 0 . 

When {x,y,z) are given, this equation, regarded as an equation in 
{x\y\z*), represents a surface. Thus each point P, or, more jjre- 
cisely, the set of surface elements through P, is transformed into a 
surface Qp, 

(b) There may be two relations 

y, z, x\ y\ z') 0, Q..;^(x, y, 2 , x\ y\ z') == 0. 

Then each point P is transformed into a curve. 

(c) There may be three relations, in which case each ]>oint P is 
transformed into a point P'. 

Let us restrict our attention to transformations of the type (a), 
and consider the effect of applying such a transformation to the 
surface elements of a given surface /S. It may happen that, no matter 
what surface jS is chosen, the transformed surface elements build up 
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another surface S'; if this is so, the transformation is called a contact- 
transformation , It can be shown that, in this case, the surface S' 
is the envelope of the surfaces tip corresponding to the individual 
points P of S. 

A contact-transformation of this type is precisely the sort of trans- 
formation which appears in Huygens’ geometrical construction. For 
each point of the wave-front S at the instant t is transformed into a 
sphere of radius c{t' — t), and the wave-front S' at the instant t' is the 
envelope of these spheres. The equations of the transformation are 

x' = x-\-c{t' —t)l, V — I, 
y* = —t)m, m' = m, 

z' = z-{-cif! —t)n, n' — n, 

§ 2.2. The analogy with dynamical systems 

In the modern theory of general dynamics, contact-transformations 
play an important part, since the history of any dynamical system 
may be regarded as the gradual self-unfolding of a contact-trans- 
formation. Consider a dynamical system in which a-re 

the generahzed coordinates specifying the state at the instant t; let 
corresponding generalized momenta. Then if 
and are the values of the coordinates and 

momenta at the instant t', the transformation of the set of variables 

corresponding set of accented 
variables is a contact-transformation. J 

Placing this result beside the statement that the simple geometrical 
construction suggested by Huygens is a contact-transformation, it is 
natural to conjecture that the general ahalytical expression of 
Huygens’ principle ought to involve a contact-transformation, 
generalizing the result just enunciated for finite systems to dynamical 
systems with an infinite number of degrees of freedom. 

Actually Huygens’ ideas have not been developed on these lines. 
The reason for this is that, in the practical applications of Huygens’ 
princii3le, radiation, the disturbance in the medium which we may 
for convenience call the aether, generally proceeds from sources. 

t The name is due to S. Lie. It was suggested by the evident fact that, if two 
surfaces are in contact, so also are their transforms by a contact-transformation. 

t For the theory of the application of contact-transformations to dynamics, see 
E. T. Whittaker, Analytical Dynamics (Cambridge, 1917), Ch. XI et seq. ; Prange, 
Encyc. der math, Wissenschaften, Band IV, Ig, Heft 4. 
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These sources are, from the dynamical point of view, singularities at 
which energy is introduced into the aether; the existence of these 
singularities prevents us from developing Huygens’ principle in the 
way which would be natural if we were dealing with a self-contained 
conservative dynamical system. 

In particular, the solution of a dynamical problem with a finite 
number of degrees of freedom is a solution valid for all values of t, 
both subsequent and antecedent to the instant whereas in radia- 
tion problems we cannot trace tlie radiation back beyond the instant 
when it issued from the source. The formulae which will be useful 
in the practical applications of Huygens’ principle will generally be 
formulae which are valid only at instants t subsequent to some initial 
instant the results obtained by substituting values of t less than 
^0 need bear no relation whatever to the actual phenomenon. 

§ 3. The propagation of sound waves in air 

§3.1. The differential equations of sound waves of small amplitude 

Huygens’ principle in its crudest form takes no account what- 
ever of the phenomenon of polarization, although this phenomenon 
was discovered by Huygens liimself in his experimental work on 
Iceland spar. Until the time of Young and Fresnel, light was regarded 
as a disturbance in a medium analogous to that of sound in air. We 
know now that the propagation of light is of an entirely different 
character from that of sound. To specify a light wave, we need 
to know the three components of the ‘light- vector’, whereas a 
sound wave is specified by a single quantity, the scalar velocity 
potential. There is, then, no precise analogy between the propagation 
of sound and the propagation of light. 

We now proceed to consider the problem of expressing in an 
analytical form the principle of Huygens for a scalar phenomenon, 
namely, the propagation of sound waves of small amplitude. This 
will serve as an introduction to the vector form of Huygens^ 
principle, which is based on the electromagnetic theory of light. 
It will, moreover, provide a justification for Huygens’ principle in 
optics as that subject was understood in the days before Young 
and Fresnel. 

We begin with a brief sketch of the theory of the propagation of 
sound waves of small amplitude. We denote the velocity of the 
medium at the point {x^y^z) at the instant t by q, where is 
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negligible. We suppose that the motion is irrotational, so that q is 
derived from a velocity potential u by the equation 

q = —gradi^, (3.11) 

where grad u denotes the vector with components {dujdx, du/dy, dujdz), 
A knowledge of q alone does not specify the state of the medium at 
[x, y , z) at the instant / ; it is necessary to know in addition the pressure 
p and the density p of the medium. Actually it is more convenient 
to consider, instead of />, the condensation $, defined by 

P “ 

where Pq is the density in the undisturbed state. In dealing with 
sound waves of small amplitude, the square of s is negligible. 

The motion is governed by two equations, a kinematical and a 
dynamical equation. The kinematical equation is the equation^ of 
continmty-f p+div(pq) = 0. 

For sound waves of small amplitude this reduces to 

divq — — s. (3.12) 

The dynamical equation is the expression in vector form of Euler’s 
Dynamical Equations, % namely 

pq+(pq*V)q = — gradp. 

It being assumed that p is a function of p alone, this equation 
simplifies to ij \ 

where the suffix 0 denotes the value in the undisturbed state. From 


this we have q = -c^grad^, (3.13) 

where ~ (dpjdp)^. The constant c has the dimensions of a velocity. 
From equations (3.11), (3.12), and (3.13) it follows that 



where is Laplace’s operator!', 


02 02 0 ^ 
dx^ ^ dy^ ’ 


t See, for example, Ramsey’s Hydrodynamics (1920), 5 (1). 

5 5 5 

t Ibid. 17 (1). The symbol V denotes the vector operator + 
the full-stop denotes the scalar product. 

11 It is the square of the vector operator V which occurs in the preceding footnote. 
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Hence the velocity potential u satisfies the equation 



This is called the equation of wave-motions, any solution of it a wave- 
function, Evidently s and each component of q are wave-functions. 


§ 3.2. Plane waves of sound 

A wave function u which depends only on t and one of the Cartesian 
coordinates, x say, satisfies the equation 

dht I dhi 
dx^ dt^' 


which is called the one -dimensional wave-equation. 


solution is 


u ~ f(ct—x)-\-F{ct-\-x), 


where / and F denote arbitrary functions. 
The particular solution 

u = f(ct—x) 


Its general 

(3.21) 

(3.22) 


is the velocity potential of a disturbance which is propagated parallel 
to the axis of x with velocity c; for the disturbance at the instant t at 
a point of coordinate x is the same as that at the initial instant 
^ = 0 at the point of coordinate x—ct. In other words, the expression 
(3.22) is the velocity potential of a plane sound wave which is propa- 
gated without change of type parallel to Ox with velocity c. The 
constant c is then called the velocity of sound. In virtue of equations 
(3.11) and (3.14), the particle velocity q (which is parallel to Ox) and 
the condensation s are given by 


qz=z cs =f(ct—x)’, 

thus c is not the velocity of the constituent particles of air but is the 
velocity of propagation of the disturbance as a whole. 

In the same way, F{ct+x) is the velocity potential of a plane 
sound wave which is propagated parallel to Ox with velocity — c. 
The general motion with velocity potential (3.21) is the resultant of 
superposing two plane waves travelling parallel to Ox with velocities 
itc. 

The arbitrary functions in the velocity potential 
<f) f{ct—x)-\-F(ct-{-x) 

of the general plane wave of sound can be determined if we are given 
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the values and of the velocity and condensation at the initial 
instant ^ = 0 for all values of x. For we have 

-=f(-x)-r{x), cs,=f\-x)+F\x), 

2/'(-a;) = ?o+C5o, 2F\x) = -Jo+c^o- (3.23) 

These equations determine / and F apart from an unimportant 
additive constant. 

From (3.23) we see that, if and Sq have general values, neither 
/ nor F is identically zero and so the wave is propagated in both 
directions. Only when one of the functions ^oitc^o is identically zero 
do we get a wave of sound which is propagated in one direction. 


§ 3 . 3 . Isotropic spherical waves of sound 
The equation of wave-motions 




1 d^u 


expressed in terms of spherical polar coordinates {r,6,<f}) has the 

^ I 4. ££^^_i 1 d^u 1 dH 

r dr r^ dd sin^0 d<l>^ dt^ ’ 


In particular, if the wave-function u depends only on r and t, it 
satisfies the equation 


d^u 2^ _ 1 ^ 
dr^ r dr dt^' 


(3.31) 


A solution of this equation is the velocity potential of sound waves 
of small amplitude in which the disturbance at any given instant 
is the same at all points of any sphere whose centre is at the origin; 
that is, a solution of (3.31) is the velocity potential of isotropic 
spherical waves. 

Since (3.31) can be written in the form 


d^ 1 d^ 


the function ru satisfies the equation of plane waves, and so the 
velocity potential of isotropic spherical waves is 

« = {f{c.t—r)+F{ct-\rr)}lr, (3.32) 

where / and F are arbitrary functions. 
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The particular wave-function 

« = \f(ct—r) (3.33) 

r 

is the velocity ijotential of isotropic spherical waves diverging from O, 
The constant c is still called the velocity of sound, since it is the 
velocity with Avhich the disturbance as a whole is propagated. But 
the actual mode of propagation differs from that of plane waves in 
that there is a change of type. 

The law of j^ropagation of the condensation s is simple. For, since 

« = ('^• 3 '^) 

the quantity rs is propagated outwards without change of typo AAith 
constant velocity c: hence s diminishes like \jr as the disturbance 
diverges from O. 

The particle velocity q is directed radially outwards and is of 
magnitude . , 

? = ( 3 . 35 ) 

r 

Thus q is in general the sum of two terms: the first, 

(3.36) 

predominates Avhen r is small; the second, 

\ript-r), (3.37) 

predominates when r is large. Thus when an isotropic expanding 

spherical wave diverges from O, the particle velocity does not undergo 
a mere attenuation but actually changes its type. 

It follows from (3.34) and (3.37) that q = cs for a diverging isotropic 
spherical wave of large radius, just as for a progressive plane wave. 
We should then expect that a progressive plane wave and a diverging 
isotropic spherical wave of large radius would have very similar 
characters; but this is not the case, as we now show. 

Let us consider an isotropic diverging disturbance which is con- 
fined to a spherical shell : inside and outside the shell the condensation 
and the velocity are zero. From (3.34) and (3.35) we have 

r^{q—cs) = f{ct—r), 
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SO that /(rj—r) vanishes everywhere except in the shell. 


(3.34), 





Hence, by 


and this vanishes if a and b are respectively less than and greater 
than the inner and outer radii of the shell. From this it follows that 
s cannot be of one sign throughout the shell; in other words a sjylimcal 
wave of positive condensation cannot exist alone:\ This does not hold 
for a progressive plane wave, since the velocity potential is constant 
but not necessarily zero in regions in which the velocity and con- 
densation both vanish. 


§ 3.4. Simple and double sources 

If r denotes distance from a fixed point P', the ex])ression 

« = (3.41) 

is the velocity potential of an isotropic .spherical wave-motion in 
which the waves expand from the (^cntl•e P'. Evidently a and also 
the velocity and condemation are infinite at P', which is therefore 
a singular jjoint of the wave-motion. Moreover, air is flowing across 
a sphere of centre V and radius r at the rate 

and this tends to ^TTf{t) as r -> 0. Hence the motion is characterized 
by the fact that air is being introduced at P' at the ratej 47 r/(^). 
For this reason we say that (3.41) is the velocity potential of a simple 
source at P' and we call f(t) the strength of the source. 

More general wave-functions can be constructed by adding together 
the velocity potentials due to several different simple sources. Let 
us consider, for example, a simple source of strength f{t) at 
P' {x\y\z') and another of strength —f{t) at the adjacent point 
{x' -{-lh,y' -\~mli,z^ -\-nh), where (/,?«.,??) are the direction cosines of 

i‘ This was first pointed out by Stokes in 1849. See Rayleigh, Theory of Sound, 
2 (1878), 99. 

X The function /(/) is not necesssarily jiositive. If f{t) is negative at any instant, 
it merely means that air is then being abstracted at P', A most important case 
arises when /(() = cos7i/, which varies periodically through positive and negative 
values. 
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the line joining the two sources. The velocity potential (3.41) of the 
source at P' is of the form 

Hence the velocity potential due to the two sources is 
u = <l>{x—x\y—y\z—z\t)—<l>{x—x'—lhyy~y'—7nhyZ—z'—nh,t) 


{ 




dx dy ' dz\ 
where the terms omitted are of the order of If we write 


m ) = ^(0 


and keep F fixed whilst we make h tend to zero, we obtain the 
velocity potential 


u == 





dzf r \ c) 


(3.42) 


which we may roughly describe as being due to two very large sources 
very close together. We call (3.42) the velocity potential of a double 
source or, more briefly, a doublet of moment F{t)\ the line through P' 
with direction cosines is called the axis of the doublet. 

For example, the velocity potential of a doublet at the origin, 
whose moment is F{t) and whose axis is the axis of x, is 


u 




From this we see at once that a double source does not emit 
isotropic spherical waves. 


§ 3.5. Poisson’s solution of the equation of wave-motions 

The simple solutions of the preceding section can be generaUzed in 
another way, namely, by considering volume distributions of simple 
sources. If the volume element dx'dy^dz' at {x\y\z') is a simple 
source of ^tvewgth f{x\y' ,z')F{t)dx'dy'dz\ the disturbance at (x,y,z) 
at the instant t is specified by the wave-function 


r 



dx'di/dz\ 


/•2 = (^x—x'f+(y—y')^+{z—z')^. 


where 
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Integrating over the whole volume V containing the sources, we 
obtain the more general wave-function 

u(x,y,z,t) == IJJ dx'dy'dz', 

it being supposed that {x, y, z) is not a point of V, If we transform 
to spherical polar coordinates with {x, y, z) as pole, this becomes 

= JJJ f{x-\-lr,y~\-mr,z+nr)F^ — ^ r sin 0 drd0d<l>, 

(3.51) 

where (Z, m, n) are the direction cosines of the hne from {x, y, z) to 
{x\y\z'). 

In particular, we may suppose that the volume distribution of 
simple sources is active only for a very short interval of time, so 
that F(t) is zero except when — c << Z <C 0; moreover, by multiplying 
by a suitable constant, we can choose F{t) so that 

/'">*= - 4 ^- 

— € 

If we now make c -> 0 in (3.51), we obtain the wave-function 

27r TT 

= A J cZ0 J d0 f{x-\rlct,y~\rmct,z-\-nct) sin0, (3.52) 
0 0 

and this may be written in the more concise form 

u{P;t) = (3.53) 

where Mpj.{f} denotes the mean value of /over the sphere of radius r 
and centre P, 

From (3.53) we have 

Hence «(P; t) -> 0, ^f(x, y, z) 

ct 

aa t-^0. Thus (3.53) is the wave-function which satisfies the initial 
conditions « = 0, dujdt ~ f whenf i — 0. 


t More precisely, this is true when / satisfies certain conditions of continuity. 
In general, the value of u when ^ = 0 and the limit as ^ 0 are not necessarily equal. 



14 THE ANALYTICAL REPRESENTATj'TON OF [chap, i 

If ^ is a wave-function, so also is dujdt since the wave-equation is 
a linear equation with constant coefficients. Hence 




(3.54) 


is a w^ave -function which satisfies the initial condition v ~ g. More- 
over, 

(3.55) 


dt dt 

Now from (3.52) we have 




27r TT 


= I; J / [l^, + «~ + n±\gi^'.y'.z'^meded4, 

r-^ct 

r%ci 

as ^ 0. Since d^Mjdt^ is hounded as ^ 0, it follows from (3.55) that 

satisfies the initial conditions = g, dojdt — 0 when t — 0. 

Finally, by adding (3.53) and (3.54), we obtain tlie wave-function 


u(P;t) ^ 


(3.56) 


which satisfies the initial conditions 




du 

¥ 




when t ~ 0. This formula, which is due to Poisson, f expresses the 
value of a wave-function u in terms of the values of u and c)u/dt at 
some fixed previous instant. 


f Poisson, Memoircs de VAcad. Botf. dea ScA. HI (1819), 121. Other proofs have 
been given by Liouville, Journal de Math. 1 (185G), 1 ; Boussinosq, Comptea rend us, 
94 (1882), 1465. The proof given here cannot easily ho put into a rigorous form: 
we justify it later as a special case of a more general theorem of Kirelihoff. Poisson’s 
formula has been extended to the case of wave motions in a space of constant 
curvature by E. Holder, Leipziger Berichte, 19 (1938), 55-66. 
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That the solution (3.56) holds only when t is positive is of im- 
portance: it arises from the fact that the solution was generated by 
sources acting only when — e <t <0, and these can evidently 
produce no effect when t is negative. 

§ 3.6. Velocity waves and condensation waves 

In the case when u is the velocity potential of sound waves in air, 
the velocity vector q and the condensation s are given by 

q — grad'ifc, 0^.9 = — . 

^ ^ dt 

Hence the function / in Poisson’s solution (3.56) determines the 
initial condensation, the function g the initial distribution of velocity. 

If the air is initially at rest, g is identically zero, and Poisson’s 
solution reduces to its first term 

u = 

The corresponding sound wave is due to an initial condensation in 
a medium initially at rest, and so may be called a condensation wave. 
If the condensation is initially zero everywhere, the solution reduces 
to its second term ^ 

the corresponding sound wave may be called a velocity wave. 

Example, w is a wave function, and U is defined by 
U(r,ct) = 

Prove ttiat lJ(r,ct) --= {(l){ct-\-r)-\-ilj{ct — r)}jr, 

where ^ and i/j are arbitrary functions. Di'duco Poisson’s formula for ii. 

§ 3.7. The verification of Huygens’ principle for expanding isotropic 
spherical waves 

The velocity potential of isotropic spherical waves with centre O 
is of the form , , 

u — ct)~{-~F2{R-\-ct), (3.71) 

where F^ and F 2 are arbitrary functions and where E denotes the 
distance from O. We can determine F^ and F 2 from the initial values 
of u and du/dt, and, as the form of (3.71) shows, an initial disturbance 
splits up in general into two isotropic waves, one converging to 0, 



16 


THE ANALYTICAL REPRESENTATION OF [chap, i 


the other diverging from O. An expanding spherical wave is obtained 
only when the initial values of u and dujdt satisfy the relation 


R du d{Ru) _ 
c dt'^ dR ~ ' 


(3.72) 


We shall consider here the case of an isotropic expanding spherical 
wave-motion in which the initial disturbance is specified by 

« = (< = 0 ), 


where F{R) is non-zero only when jBg < i? < i?i; thus the initial 
disturbance is null except in the shell bounded by the two spheres 
8^ (B, = JBi) and 8^ {R = R^, Then by (3.72) the initial value of 
dujdt is ^ 

P(J8) 


in the shell and is zero elsewhere. We show that Huygens’ descrip- 
tion of the propagation of this disturbance is in agreement with the 
analytical solution given by Poisson, namely, 


where f=-^F'(R), g ^ Lf{R), 

and M denotesf a mean value taken over the sphere 8 with centre P 
and radius ct. (See Fig. 1.) 



Let Rq be the distance of P from the origin. Then if P is outside 
/Si, the disturbance at P is null, except possibly during the interval 

t We omit the double suffix P, ct indicating the centre and radius of the sphere 
over which the mean values are calculated. No confusion will be caused by this. 
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(Ro-M^yc ^ t ^ (i?0-f-^i)/c; for at an instant outside this interval 
of time the sphere 8 of radius ct does not cut into the shell to which 
the initial disturbance is confined, and so the mean values of / and 
g are both zero. Similarly, if P is inside 8^, the disturbance is null 
except possibly during the interval {B 2 —Bq)Ic < ^ < (Bq+B^)Ic. 
It remains to consider Poisson’s solution when 8 does cut into the 
shell bounded by 8^^ and 8^. 

Let 6 be the angle between any radius vector r at P and the line 
FO; we then have 


_ iJ2+;.2_2Po^cos0, B^r^inOdd = BdB. 


Hence 


M{g} = ^ J J 9^ sin 0 d6d(l> = i J gr sin 0 


de 


1 f ^ « 
n* 


where the limits of integration will be found later. This gives 


Similarly, 




2B, 

Substituting in Poisson’s formula, we have 


o' 


( 3 . 73 ) 

2i?(, \ c 8t f 2Ro \ c 8t j ' ' 

There are tliree cases to be considered, viz. (a) P outside 8^, (6) P 
inside 8^, (c) P between 8^ and 

(a) P outside 8^\ {B^ > P^). 

(i) ct < Po“^i- Then = 0 since 8 does not cut into the shell. 

(ii) Bq—Bi cct <Z Bq—B^, In this case, we have 


and so 

4617 


B[ = Pi, Pg = Pq— c^. 


u{P\t) 


F(B^-ct) P(Pi) 
PQ 2 Pq 



18 THE ANALYTICAL REPBESENTATION OF [chap, i 

This agrees with the known value 

only if F{R^ = 0 , i.e. if F{R) is continuous at i? == R^, We shall, 
in what follows, assume that the initial data are continuous, f so 
that F{R^) -= F{R^) = 0 . 

(iii) i?o“^2 < < i?o+^2- Then R[ -= R^, J?' ^ iJg, and so 

u(P\t)-==Q. 

(iv) JZQ+iZg <ct < Rq-\-R^. Then R[ — R^, iZ' = cI—Rq] again 
u{P;t) = 0 . 

(v) < cZ. In this case u(P;t) = 0 since S does not cut 
into the shell. 

Thus Poisson's formula verifies Huygens’ geometrical construction 
in the case when P lies outside the shell, provided that the initial 
data are continuous. 

(b) P inside 82; (Rq < iZg)- 

(i) ct < R.^—Rq* Then u = 0, since 8 does not cut into the 
shell. 

(ii) R2—Rq <ct < Ri~Rq. In this case we have 

R2 ~ ■^ 2 > ™ — Fq, 

since F(R^ — 0 by hypothesis, we have u{P;t) — 0. 

(m)X R^-RQ<ct <R2+R^ Then R[ = R,] R2 = R^. and 
u(P;t) — 0 . 

(iv) R2+R0 <ct < Ri+Rq. Then iZ; = iZ^, iZ^ == ct—R^^; hence 
u{P;t) = 0 . 

(v) Rq+Ri < ct. In this case, the effect at P is null, since 8 does 
not cut into the shell. 

Poisson’s formula does not give rise to a returning wave, as the 
crude geometrical form of Huygens’ principle would. Thus in the 
present analytical formulation there is no need of a special hypothesis 
discarding the inner sheet of the envelope of the secondary wave- 
fronts. 

f I'or a discussion of wave-motions with discontinuities at the wave fronts, see 
Love, Proc. London Math. Soc. (2), 1 (1903), 37-62, 291-344 ; Rayleigh, ibid. 2 (1904), 
266-9. The present example is taken from the first of these papers. 

t We assume hero that Rq > i(Ri— i? 2 )' It fhe argument needs slight 

changes. 
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(c) P lies between /S, and (i?2 < jRj < ^i)- 

A similar argument shows that 

u(P;t) = ±F(Eo-ct) 

Hq 

until the instant t = (ii^— i?2)/c, and that the disturbance is there- 
after null at P, 

The disturbance we have just considered is confined at the instant 
t to the shell bounded by the sj)heres R ^ R = R^-\-ct, In 

particular, if R^—R^ is very small, we have in effect a solitary 
spherical expanding wave of the type considered by Huygens. The 
initial disturbance is confined to a very thin layer on the surface of 
the spliere At the subsequent instant t the disturbance is confined 
to a layer of tlie same thickness on the surface of the sphere 
R R.^~\-ct, which is the outer sheet of the envelope of spheres of 
radii ct, whose centres lie on Thus Huygens’ geometrical construc- 
tion, with its restriction that only one sheet of the envelope is to be 
considered, is justified by Poisson’s analytical solution of the equation 
of wave-motions. 

Moreover, Huygens’ statement that a secondary wave is of effect 
only at the point where it touches the envelope also follows. For 
wlien the layer is very thin, Poisson’s integrals for !/{/} and M{g} 
are extended over a very small area on S near the point where 
OP cuts ASg*, and if P is in the wave-front at the instant t, the 
secondary wave, with centre Pg radius ct touches the envelope 
at P. 

In a similar wayf we can justify Huygens’ construction for an 
isolated spheri(‘al converging wave or an isolated plane wave; but 
only if the initial values of the velocity potential and condensation 
are suitably chosen. In other words, to justify Huygens’ principle 
for isolated weaves, we must have recourse to analysis and take into 
account the dynamics of the medium in wdiich the wave-motion 
occurs. 

Example. Apply Poisson’s solution of the ecpiation of wave-motions under 
th(^ initial conditions 

“ ^ < ^i)' =0 (R> R < Rt), 


+ Sco Croze, Annalcs de Physique, 5 (1926), 370-439 (383-98). 
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where jP(JR) is continuous at and jRj* Discuss the nature of the solution 
so obtained. (Love.f) 

§ 4. Huygens’ principle for monochromatic phenomena 

§4.1. Fresners extension of Huygens’ principle 

In his memoir on Diffraction, which won the Paris Academy’s prize 
in 1 8 1 8, Fresnel J made an important extension of Huygens ’ principle, 
in that he replaced Huygens’ isolated spherical waves by purely 
periodic trains of spherical waves and made use of the principle of 
interference. On this theory, light ought to appear, not necessarily 
on the envelope of the secondary waves, but at every point where 
these secondary waves reinforce one another; on the other hand, there 
should be darkness wherever these secondary waves destroy one 
another. In this way, Fresnel was able to account, not only for the 
rectilinear propagation of light of very short wave-length and the 
laws of reflection and refraction, but also for certain diffraction 
phenomena. 

Fresnel, at this stage, still regarded light as a disturbance in an 
aether analogous to sound in air. He did not yet realize that the 
phenomenon of polarization made such a theory untenable. (His 
elastic-solid theory of the luminiferous aether, in which light con- 
sisted of transverse vibrations, dates from 1821.) Accordingly, we 
shall discuss here the theory as applied to sound waves of small 
amplitude; the application to the scalar theory of light is merely 
a matter of changing the terminology. 

Let us consider sound waves of a very general character, which 
are generated by making a certain part of the medium execute forced 
vibrations, not necessarily ‘monochromatic’. Such a wave-motion is 
evidently much more complicated than that which arises when a 
portion of the fluid is initially disturbed from its equilibrium state 
and the disturbance is allowed to propagate itself freely; for a wave- 
motion of the former type is the result of superposing an infinite 
succession of wave-motions of the latter type. It seems likely that 
in a wave-motion produced by the forced vibrations of a portion V 
of the fluid, the effect at the instant t at the point P outside V would 
depend on the velocity and condensation at each point $ of F at 

t Loc. cit. 44. 

X M6m, de VAcad, 5 (1826), 339; reprinted in Frosnol’s Q^uvrett completes, 1, 247. 
For the history of this development, see E. T. Whittaker, History of the Theories 
of Aether and Electricity (1910), 113-16 and Croze’s paper already cited. 
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the instant t—PQjc, This idea, which is implicit in the work of 
Fresnel, t attained a rigorous analytical formulation in Kirchholff’s 
integral theorem. 

In the case of ‘monochromatic’ forced vibrations the matter is 
somewhat simpler. Let us consider for definiteness the case con- 
sidered by Fresnel, namely that of a ‘monochromatic’ source of 
expanding waves. If the source is at the origin, the wave-motion is 
characterized by the velocity potential 

u = ( 4 . 11 ) 

T 


or, rather, by its real part.J In this motion the sphere S, whose 
equation is r — is a wave-surface; for all the particles of the fiuid 
on S at the instant t are being affected by the disturbance which 
left the source at the instant t — ^q/c. 

The disturbance which reaches S at the instant t continues to be 
propagated outwards and reaches the wave-surface of radius 
at the instant , v , i . f . {q 




It seems natural to regard the vibrations of the particles on S' at 
the instant t' as being due to the vibrations of the particles on S at 
the instant t, instead of being due to the pulsations of the source 
at the instant t'—rQjc, 

This is, in fact, the assumption which Fresnel made. He supposed 
that each element of the wave-surface S acts as a secondary source 
which is sustained by the displacement and velocity given to the 
particles of the surface element by the primary wave. The resultant 
effect is produced by the interference of thes^ secondary waves. 

If Fresnel’s theory is to provide a valid extension of Huygens’ 
principle, the secondary sources must produce not only the correct 
effect outside S but also a null effect inside S, Fresnel believed that 
this could be done only by taking account of the dynamical effect 
of the condensation and of the velocity in the primary wave at each 
point of S; that is, by supposing that each secondary source emits 
a ‘condensation wave’ and a ‘velocity wave’. Thus an analytical 
formulation of Huygens’ principle which involved only one of these 


t See Croze, loc. eit. 398-405. 

i As is usual in this sort of work, it is simpler to work with complex wave-functions, 
and then to take real parts at the end of the analysis. 
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two types of secondary wave would give rise to a non-null effect 
within S and would have no physical justification. 

Unfortunately Fresnel was unable to carry out this programme 
completely. He considered only the case when the wave-surface S 
is of large radius compared with the wave-length A. Tn this case (cf. 
§ 3.3) the velocity q and the condensation s at ])oints of S are con- 
nected approximately by the relation 

q = cs. 

He thought that it would then suflice to attribute to the elements 
of 8 the same velocity as is communicated to them by the primary 
wave in order to obtain the correct eflect at ])()ints outside 8 and at 
a large distance from it. Actually considerations of this nature lead 
to incorrect results xinless one makes, as Fresnel did, the following 
additional assumptions: 

(i) the elements of 8 execute vibrations Avdiose amplitude is to tlie 
amplitude in the primary wave as is 1 : A; 

(ii) the elements of 8 arc oscillating a quarter of a period ahead 
of the primary wave. 

The necessity for these two additional assumptions*!* led many to 
regard Fresnel’s theory merely as a convenient means of calculation 
which lacked any sound physical basis. 

Various attempts have been made to overcome the phase difficulty. 
For example, (^euyf replaced each clement (18 by a small pulsating 
sphere (of surface area 2(18), which emitted monochromatic isotropic 
waves. This does give the correct phase; but as it also gives rise to 
an effect inside 8 which is not null, it cannot bo accepted as a satis- 
factory formulation of Huygens’ principle. 

Nevertheless, it is the case, as we shall see in § 4.4, that a careful 
analytical formulation of Fresnel’s original idea, that each clement 
of 8 gives rise to a condensation wave and a velocity wave which 
can be determined on dynamical ininciples, leads to correct results 
without any arbitrary assumptions regarding the phase and ampli- 
tude of the secondary sources. Moreover, the formula obtained in 
this way is in complete agreement with Helmholtz’s analytical 

•f Frefincl also assumod that tho offoct at a point P outside S is due only to the 
secondary so\irccs on the part of S visiblo from P, and not to the secondary sources 
all over *S'. This assumption makes litllo difforenoo when tlio wave-length is small. 

} Contptes renting. 111 (1890), 910-12; Annalea de Physique ( 6 ), 24 (1891), 145. 
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formulation of Huygens’ principle for monochromatic phenomena, 
which we now discuss. 


§ 4.2. Helmholtz’s formula 

The velocity potential of sound waves of small amplitude satisfies 
the equation of wave-motions 






If the disturbance is ‘monochromatic’, u will be of the form 
where v is independent of t and is a solution of 


(V2+P)i; z:. 0. (4.21) 

This equation reduces to Laplace’s equation when i = 0; and the 
formula of Helmholtz,*]' which we are about to prove, then reduces 
to the well-known Green’s equivalent layer formula in the theory of 
attractions. 

We start with Green’s identity 


Jjf J = J J (*" ^-^1;) (4.22) 


where H is a closed surfiicc bounding the volume V and djdv denotes 
differentiation along the outward normal to S, This identity is 
certainly validj when v and iv and their first- and second-order partial 
derivatives are continuous within and on S, or, as we shall say in 
future, when v and w are regular within and on 8. If v and w satisfy 
(4.21) and the prescribed conditions of continuity, (4.22) becomes 


in-s- 


dS = 0. 


dv\ 


In particular, if w = where r denotes distance from the fixed 
point P, we have 


provided that P and the singularities!! of v lie outside 8, 


t Journals. Math. 57 (1859), 7. 

X Hero, and throughout this book, we tako the simplest sufficient conditions for 
the validity <.)£ Green’s identity. Lighter sufficient conditions can be obtained as in 
Kellogg’s discussion of the Divergenco Theorem {Foundations of Potential Theory 
(Berlin, 1929), Ch. IV). 

II By a singularity of v wo moan a point at wliich v or one of its first or second 
partial derivatives is discontinuous. 
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If P lies inside S but the singularities of v are still outside, equation 
(4.23) no longer holds, since w becomes infinite at P. To avoid this 
difficulty, we apply Green’s identity (4.22) to the volume bounded 
externally by 8 and internally by the sphere o*, of centre P and small 
radius e. This gives 



a 



a 



a 


since the integral over 8 is independent of c. But since v and dvidr 
are continuous at P and the element of surface on a is e^midd9d<f)y 
the value of this limit is — 47r?;(P). 

We have thus proved the following theorem of Helmholtz: 

Let V be a solution of == 0 whose first- and second-order 

partial derivatives are continuous within and on a closed surface 8^ 


where r is distance from a fixed point P and djdv denotes differentiation 
along the outward normal to 8, Then the value of /(P) is 47ti?(P) or 
zero according as P lies inside or outside 8, 

In this theorem the sources of the disturbance specified by the 
* monochromatic’ wave-function ve-^^^ lie outside a certain closed 
surface. In most of the applications all the sources lie at a finite 
distance and the effect at a distant point is to be expressed as an 
integral over a surface containing the sources. Accordingly we must 
modify the theorem to cover the case when all the sources lie inside 
a closed surface. 

The function v is now a solution of = 0 whose first- and 

second-order partial derivatives are continuous on and outside a 
closed surface 8. We choose R so large that the sphere 2, whose 
equation is r = P, encloses 8; and we denote by the volume 
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bounded externally by S and internally by S. Then if the normal 
V is drawn into the value of 


J{P) = 




dv d 

T dv dv 


.S’ L 



dS 


is ~^7Tv(P) or zero according as P is or is not a j)oint of 1^. 

The normal v at jDoints of S is the ordinary outward normal. 
Hence, in the notation of (4.24), 5'^^ ^ y . { v ^ 

T{P)— JJ — tTujj ^dco — JJ da>, 

where dco is the solid angle subtended by an element of the sphere S 
at its centre P. Now make i? -> oo. The term 



dco 


tends to zero if y -> 0 as r -> oo, uniformly with respect to the polar 
angles 0 and <j6; in ])artigular, it tends to zero if 

|n’l < asr-->oo. (4.25) 


The term 


tends to zero if 


//'•“Kl’-Hl-/” 


as r 00 


(4.26) 


uniformly with respect to 6 and Condition (4.25) is called by 
Sommerfeldt the ‘Endlichkeitsbedingung’, condition (4.26) the ‘ Aus> 
strahlungsbcdingung’. We shall return later to point out the im- 
portance of (4.26); at tlie moment, we merely remark that, when it 
is satisfied, the wave-function ve is the velocity potential of a 
system of expanding waves. 

When (4.25) and (4,26) are satisfied it follows that the value of 
I(P) is — 47ri;(P) or zero according as P lies outside or inside the 
closed surface S containing all the sources. Hence we have the 
following theorem: 


t Jahresberichi der D,M,V. 21 (1912), 309-53 (326-34). 
B 
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Let V be a solution of = 0 whose first- and second-order 

partial derivatives are continuous outside and on a closed surface 8, and let 

\rv\ <K, 


'■(S “***") 


uniformly with respect to 6 and (f) as r oo. Let 

>Sf 

where r is the distance from a fixed point P and djdv denotes differentia- 
tion along the o\dward normal to 8. Then the value of I{P) is —47 tv{P) 
or zero according as P lies outside or inside 8, 

Example, (fy and ift arc two solutions of = 0 which satisfy the 

following conditions: 

(i) (f> and its first- and second -order partial derivatives are continuous within 
and on the closed surface S; 

(ii) ijj and its first- and second-order partial derivatives are continuous 
outside and on S, and i/j satisfies Sonimerfold’s conditions at infinity. 

Prove that, if cf) — i/j on S , the value of 




is 4i7T(l>{P) or 47 t^(P) according as P is inside or outside S. Show also that 
if dify/dv ~ difildv on S, the value of 




is ~-47 t^(P) or — 47ri/f(P) according as P is inside or outside S. 

§ 4 . 3 . Restrictions on the use of Helmholtz’s formula 
Helmholtz’s formula 




expresses the value of the solution v of 

+ = 0 , 


(4.32) 


regular within a closed siirface S, at any point P inside S in terms 
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of the values 0 and Y taken by v and dv/dv on 8. We can, however, 
express v(P) in terms of <f> alone, by the equationf 

where O denotes the Green’s function with singularity at P. 

We see from (4.33) that a knowledge of the boundary values of v 
alone on 8 determines v inside 8, and, moreover, determines v in 
general uniquely. In particular, we can find from (4.33) the boundary 
values of dv/dv on 8. Thus a knowledge of O alone determines T, 
and in general does so uniquely. Conversely, it can be shown that 
a knowledge of T alone determines 3> and, in general, does so uniquely. 
Hence the functions d> and Y in (4.31) are related and cannot be 
assigned arbitrarily and inde|)endently of each other. 

Of course, when O and Y are arbitrarily assigned, the expression 
on the right-hand side of (4.31) does satisfy the differential equation 
(4.32), since it is obtained by the addition of particular solutions of 
the form and its normal derivative. But these arbitrary values 
of O and Y are no longer the boundary values of v and dv/dv. 

There is, however, an exceptional case to which the theory of 
Green’s function is not api)lieable, and in this case a knowledge of 
one of the functions and Y does not determine the other. To see 
how this arises, let us suppose that u ~ v cos kct is the velocity 
potential of sound waves of small amplitude, so that v satisfies 
equation (4.32). Then if is a rigid boundary, the air inside 8 
possesses certain normal modes of vibration, corresponding to a 
sequence A'g,.- of values of k, the eigenvalues of the problem. 
The velocity potential u = cos kj. ct corres})onding to one of these 
normal modes satisfies the boundary condition dVj.ldv ~ 0 on 8. 
Hence the problem of finding the solution of 

{V^+kf)v ^ 0 , 

regular within 8, given the boundary values of dvjdv on 8, is an 
indeterminate problem; for, if V is one solution, 

V V+Av^ 

is also a solution for all values of the constant A, Similarly, there 

t SeePockels, Vber die part wile Dfjfcrentialglcichunff (A + k-)u — 0 (Leipzig, 1891), 
p. 280. The equation (4.33) does not provide an analytical formulation of Huygens’ 
principle since, as we shall see later, it does not give a null effect outside S. 
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exists another sequence of values of Ic such that the problem 

of finding the solution of 

(V2+t;2)v == 0, 


regular within S, given the boundary values of v on S, is also an 
indeterminate problem. 

To sum up, except when the constant k is one of the eigenvalues 
and i*', a knowledge of one of the boundary functions O and T 
in (4.31) determines the other uniquely, and together they determine 
V uniquely. 

These exceptional cases do not occur in the exterior problem. It is 
true that there are solutions of the ccpiatioii (7^~\-k^)v — 0 , for certain 
special values of k, which are regular oiitside S and satisfy the 
boundary condition v = 0 (or dr/dv = 0 ) on S and the condition of 


finiteness 


Ircl < K 


(4.34) 


at infinity. But all these solutions (eigenfunctions) re])resent standing 
waves and do not satisfy Sommerfekrs Aitsstraklungsbedinguvg, 
namely that /o v 

( 4 - 35 ) 


uniformly with rcs])ect to the ])olar angles 6 and (/► as r co. 

By imposing the additional (condition (4.35), we can apply the 
Green’s function argument to the exterior problem and assert that, 
for any value of A’, a knowledge of one of the boundary values O or 
T*, taken on S by a solution v of (4.32), regular outside S, and its 
normal derivative, determines the other uniquely, and that together 
they determine v by means of the formula 






d*S'. 


§ 4 . 4 . The physical meaning of Helmholtz’s formula 
Let u ™ 

be the velocity potential of sound waves in air due to 'monochro- 
matic’ sources all lying inside a certain closed surface S. Helmholtz’s 
formula expresses the value of at a point P outside S as due to 
a certain distribution of simple and double sources over S. Our proof 
of this was of an entirely analytical character; we now give an 
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alternative proof by means of the dynamical arguments suggested 
by Fresnel. t 

Let Q be a typical point of 8, at a distance from P; let v be the 



Fig. 2 


unit vector drawn along the inward normal to 8 at Q. The particle 
velocity at Q lias a component in the outward-normal direction 

dv 

Hence air flows across the element d8 at Q at the rate 

dv 

In addition to jiroducing a flux across d8, the sources within 8 
also cliaiige the ])ressure at Q from the equilibrium pressure p to 
p+Sp. If p is the e(piilibriuni density of air, 8p is given by 

Sp — /?c%, 

wliere .s* is the condensation, and so 

8p == —pikever^^^’^. 

Hence there is a thrust on the area d8 at Q oS magnitude 

—pikcve~^^^'*d8. 

Now siq)])osc that the sources and all the air inside 8 are 
destroyed. In order to get the effect outside 8 specified by the 
velocity potential we must introduce sources over 8 with the 

following properties: 

(a) Air is created at each element d8 at the rate 


dv 
— € 
dv 




t (Vozo, Annates dc PhysiqtiCy 5 (1926), 370 (408-11). Seo also Larmor, Proc, 
London Math. Soc. (2), 19 (1921), 169-80. 
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(6) A force dS perpendicular to dS acts on the air 

in contact with dS, 

The creation of air at dS makes this element a source of strength 


47r dv 


e-^^'dS; 


the corresponding velocity potential at P is 


(4.41) 

47rri dv 

It is well knownf that a concentrated force acting at Q 

along the outward normal to S gives rise to a sound wave of velocity 
potential ^ 

4:7rpkc dv\ T*! J 


Hence when the force is 
potential at P is 


given by (6) above, the resulting velocity 


47T [ Vi I 


(4.42) 


The total effect of the sources on dS is obtained by adding the 
velocity potentials (4.41) and (4.42). Of these, the first is determined 
by the normal velocities of the air particles on dS and so represents 
the velocity wave. The second depends on the condensation in the 
primary wave at Q and so represents the condensation wave (§ 3.6). 

The effect of all the sources on S is obtained by adding (4.41) and 
(4.42) and integrating over S; this gives 




s 

which is Helmholtz’s formula in the case when P lies outside a closed 
surface S containing all the sources. By a similar argument we can 
show that Helmholtz’s formula also holds when P lies inside a closed 
surface S containing none of the sources. 

Fresnel believed that, if all the sources were inside 8, the secondary 
sources on each separate element dS would produce a null effect at 
points within 8. This is evidently not the case. The effect inside 8 
produced by the action of all the secondary sources on 8 is however 
null, as can be seen by the following simple argument. 

Suppose that P is a point inside 8, Introduce a diaphragm G 


f Lamb, Hydrodynamic (Cambridge, 1916), 496. 
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dividing the space inside S into two regions, of which one contains P 
and the other contains all the sources. The rim of G divides S into 
two areas 



ff 


Em. 3 


/S' and S", as shown in the figure. We now apply Helmholtz’s 
formula to the two surfaces S'-\-G and /S"+(7. This gives 

„ 1 r r dv d 

( 4 . 43 ) 


«' + C 


= e; J J I'-v;- - a - v^- r. i) 

^**4* C* 


Evidently on C the outward normals to the closed surfaces S'-\-C 
and S"+C are in opposite senses. Hence, when we add (4.43) and 
(4.44), we obtain 


Iff 0 (emri-ct)\\ 

I V w) 

/Sf 


dS; 


in other words, the effect inside S is null. 

Apart from the fact that the secondary source on each separate 
element dS does not produce a null effect inside S, we now see that 
Helmholtz’s formula presents an analytical form of Huygens’ prin- 
ciple in the sense understood by Fresnel. It is not merely a convenient 
analytical formula, but has a sound dynamic^-l basis; the secondary 
sources are real in the sense that their amplitude and phase are 
determined dynamically by the velocity and condensation in the air 
at points of the surface S, 

If it were possible to find a distribution of sources on S, different 
from that of Helmholtz, which would give the correct effect at points 
outside S and a null effect inside S, such a distribution would have 
no physical meaning whatever, since it could not be obtained by 
considering the dynamics of sound waves in air. Actually no such 
alternative distribution exists. For if there were such a distribution, 
we could obtain by subtraction another distribution which would 
give a nuU effect everywhere. But a distribution of sources on S 
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giving a null effect everywhere would itself be null, contrary to our 
assumption. Hence Helmholtz’s distribution of sources is the only 
one with the required property. 

§ 4.5. Retarded values 

It is often convenient to make use of the idea of the 'retarded 
value’ of a function. If ^ is a function of the coordinates {x,y,z) 
of a variable point Q and of the time t, say 

^ = 4>{x,y,z,t), 

and if r is the distance of Q from a fixed i)oint P, we write 

bf>] = 

and call [<^] the retarded value of (f>. 

§ 4.6. Stokes’s diffraction formula and the theory of Fresnel 

We now derive from Helmholtz’s integral an important approxi- 
mate analytical formulation of Huygens’ principle, due to Stokes; 
this approximation is valid when the wave-length is very small 
compared with the other distances involved. 

Let us consider in the first instance the case of ex])anding mono- 
chromatic isotropic s])herical waves Avith complex velocity potential 

= . 

r 



If we take S to be the wave-surface r ■== t-q, it follows from Helm- 
holtz’s theorem that the value of 

477 JJ \ r, Sro\roj drX~^\~n 
8 
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is u{P\ t) or zero according as P lies outside or inside S, This integral 
transforms into 

I = — L f f cos ^) - 1 + 1 cos dS, 


where 


dr. 


cosi/f cos LPQn ~ — — - 


dTfi 


Now, by hypothesis, the wave-length ^Trjk is small compared with 
Vq and and so l/r^ and l/r^ are negligible compared with k. This 
gives the approximate formulaf 

I = _ JL r f — (1+cos.A) dS, (4.61) 

47rJJ 

s 

which is substantially FresneVs diffraction formula. 

^ikiXQ—ct) 

But since 

on /S, we have there 
du 
dn 

since we are neglecting 1 /r^ in comparison with k. Hence 

^ik(ro+ri-ct) 


U = 


^ik(rQ-ci)j j \ ^ ^ik{ro~ct) 

= [ilc 1 z= i]{ , 

U \ ^o/ »-o 


’^1 = a*- 

dn\ 


Substituting in equation (4.61), we have 

j_ rr [dui i+cosi/j 

47r J J [a?ij Ti 


u{P-t) 


dS, 


(4.62) 


when P is a ])oiut of the region into which the wave-surface P is 
progressing. 

The case of a monochromatic progressive plane wave-motion is 
rather more difficult. We start with a complex velocity potential 

U = e^k{x-ci)^ 

but have to take real parts before proceeding to a certain limit. We 
shall denote the real part of U by u. 


t Note that, if A is the wave-length, the contribution of the element dS is of the 


form 


^ikra—iiri ^ik{rx—ct) 

(1-fcos^) dS, 

2Aro ' ^ ri 


At first sight this seems to imply, as Fresnel believed, that the secondary sources 
oscillate a quarter of a period ahead of the primary wave. 

4617 y 
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Without loss of generality we may take the wave-surface S to be 
the plane a; = 0; if P has coordinates there are two cases 

to be considered according as is positive or negative. We apply 
Helmholtz’s theorem, taking as the surface of integration 8 the 
boundary of the volume in which R and x'^O: 8 consists of 



x^poaiiive Fio. 5 



a plane area 8^ and a portion of a sphere 8^, as shown in the two 
figures. In either case Helmholtz’s integral is 


/ = 


1 r /* a 

J J i rj dv 

«,+s, 





d8, 


and this has the value U(P;t) or zero according as is positive or 
negative. 

The integral over is 




ii/f? 


^ f) ] 

dx I j 


d8. 


If we neglect l/r^ in comparison with k, we obtain the approximation 


where square brackets again indicate a retarded value. 

The integral over 8^ is 

/. = * f r l'fl‘ dS 

injj I r, at «•, 
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Now the integrand has a constant value on all the circles in which 
S 2 is cut by planes parallel to S^, and so we can carry out the integra- 
tion by dividing into zones bounded by these circles. The area 
on ^2 bounded by the two circles x ^ x ^ is 27 rRh^, 

and 7*1 r- R on /Sg ; hence 


^ik(R+xi-rt) 

2R 

fjik(R+xi-H) 

I 


41. 

J e^^i{ik^—ikR-{-l} 

-fi 

= e^-kiR 


Combining these formulae for 7^ and I 2 , we have the approximate 
formula 

*Vi 

We cannot make R tend to infinity in (4.63), since the second term 
on tlie right-hand side does not tend to a limit. This difficulty can be 
overcome by taking real parts and choosing special values for R. 

[f we write u ~ cos lc(x—ct), it follows from (4.63) that 


47r J J Ti [Sm 


+ '^^^cosA-(i?-cO, 


(4.64) 

if P is a point of the region into which the wave-surface S is ad- 
vancing. If we give R the special value ct+(m+^)TTlk, where m is 
a positive integer, (4.64) becomes 

Hence, if we make 7ri tend to infinity, 

-S’ 

where the integration is, in a sense, over the whole wave-surface S 
but has to be calculated by the special limiting process indicated. 

We have thus proved that, for monochromatic plane or spherical 
waves with velocity potential u, 

where P is a point of the region into which the wave-surface 8 is 
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progressing. This formula is usually known as Stokes's diffraction 
formula , For the application of the formula to the approximate 
solution of diffraction problems, we refer the reader to the standard 
text-books. 


§ 5. Wave-motions in three dimensions 
§ 5 . 1 , Kirchhoff^s formula 

In the earher part of this chapter we have given the two analytical 
formulations of Huygens’ principle for sound waves associated with 
the names of Poisson and Helmholtz. Poisson’s formula justifies the 
principle in the case of an isolated wave due to an initial disturbance: 
Helmholtz’s formula holds in the case of a ‘monochromatic’ distur- 
bance. We now show that these formulae are particular cases of a 
general theorem due to Kirchhoflf, concerning sound waves of any 
structure and origin. 

We have seen that, if 


is a ‘monochromatic’ wave-function with no singularities within or 
on the closed surface S, then the value of u at a point P 
within S at the instant t is given by 




pikir~cl) 

^ ^ dS. 

r dn) 


where r is the distance from P to a typical point (x, y, z) on S and 
djdn denotes differentiation along the inward-drawn normal to S, 
The value of the integral is zero when P is outside S, 

Introducing the ‘retarded’ values J of the various functions in- 
volved, we find that Helmholtz’s formula can be written as 


or, finally. 


since — = —ikcu and therefore 
dt 


.¥j ” 


ikc[u\. 


t It was given by Stokes in his memoir on the ‘ Dynamical Theory of Diffraction ’ 
(Trans, Camb, Phil, ^Soc. 9 (1849), 1; Math, and Phys. Papers^ 2, 243). 
t See § 4.5. 
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The formula we have just obtained is due to Kirchhoff. It ex- 
presses u[P\t) as a surface integral in which the period 'Zn/ikc) does 
not occur explicitly, and so it is true for any period. Now an arbitrary 
function of t can be expressed as a sum of periodic constituents by 
means of a Fourier series or a Fourier integral. It follows that, since 
Kirchhoff’s formula is linear in u, it holds for any solution of the 
equation of wave-motions, not merely for solutions corresponding to 
monochromatic disturbances.! 

We have thus obtained the following theorem 

Let 2 i(x, y, z, i) be a solution of the equation 


= i ^ 

bt^ 


v^hose partial derivatives of the first and second orders are continuous 
within and on a closed surface aV, and let {x^,yi, z^) be aimnt within N, 
Then 







1 br 'dll 1 du 1 
cr bn bt r dn_ j 


where r is the distance from {x^yVxyZf) to a typical point of S, b/bn 
denotes differeritiation along the inward normal to S, and square 
brackets indicate retarded values. If, however, the point lies 

outside S, the value of the integral is zero. 

This theorem can also be applied when u has no singularities in 
the volume V bounded internally by a closed surface S and externally 
by a sphere 2 whose equation is oci^-\-y^-\-z^ = R^\ if we denote by 
bjbn differentiation along the normal drawn^into V, the value of the 
sum of the integrals over S and 2 is or zero according 

as {x^, y^, Zj^) is or is not a point of V.\\ 

Now make i? oo; Ave then obtain a theorem valid when u has 
no singularities outside a closed surface S for all values of t from 
—00 up to the instant under consideration, provided only that the 


t We do not go into a detailed application of the theory of Fourier series or 
integrals, as it is possible to give a direct proof of Kirchhoff’s formula without using 
Helmholtz’s formula. See § 6.2. 

I See Kirchhoff, Berliner Sifzungsber. (1882), 641; An7ialen der Phys. 18 (1883), 
663 ; Vorlesungen u. math. Phys. 2 {Optik), 23. The theorem has been generalized by 
W. R. Morgans in Phil. Mag. 9 (1930), 141, to cover the case when S is & moving 
surface. 

II It must not be a point on S or 2. 
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integral over S tends to zero. This last condition is certainly satisfied 
if If behaves like /(d—JB)/i2 for large values of R, where /(cO and/'(r^) 
are bounded near t ~ —oo. The result so obtained is: 

Let u(Xj y, z, t) be a wave-function which has no singularities oiitside 
a closed surface S for all values of t front — oc wp to the instant under 
consideration, and which behaves like f{ct — R)/R at a large distance R 
from the origin, where f{ct) andf\ct) are bounded near t = — oo. Then 
if is a point outside 8, Kirchhoff's formula holds, provided 

that djdn ineans differentiation along the outward normal to 8. If 
is inside 8, the value of the integral is zero. 

These two theorems constitute Kirchhoff’s integral formulation of 
Huygens’ principle. As in the case of Helmholtz’s formula, the 
boundary values of u and dujdn on 8 cannot be arbitrarily assigned 
independently of each other. (See § 4.3.) 


§ 5.2. A direct proof I of Kirchhoff ’s formula 

Let S' be a closed surface within and on which the function v(x, y, z) 
and its partial derivatives of orders one and two are continuous, and 
let r denote the distance from a fixed point F outside 8. 

Then, by Green’s transformation, we have 


»S' V 

where V is the volume bounded by 8 and v the normal to 8 directed 
out of F. 

Let u{x, y, z, t) be a solution of 




1 d^u 


which has no singularities inside or on 8, and let us take 
V = u[x,y,z,t—^^ -= [ m ] 

in the formula (5.21). Differentiating with respect to x, we have 


dv 

dx 



(5.22) 


t This proof is somewhat similar to one published by Gutzmor, Journal f. Math. 
114 ( 1895 ). 333 . 
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and therefore 


dv r du\ \ dr \ du 

dv dv c dv dt 


From (5.22), we have 


(5.23) 


Ydx^\ c 0a;[0a;0^J’ 

“ [dxdtjc'^iWy 

and so 

dh) r d'^n\ 2 dr F d^u 1,1 (^^Y\ ^^^1 ^ [ ^^1 

dx^ idx^ J c dx j [fix) [ J ^ L J 

Adding the three equations of this type, we obtain 



1 du d^r 
c dt 2^ dx^ 


where 2 denotes summation over the three variables x, y, z. 

The integrand of the volume integral in equation (5.21) is therefore 


1 ^ 2 ^ 2 Xd^u 2 du 2^x—x^‘d^u' 

r ^ cV [ dt^ cr^ dt c 2i r^ dxdt ' 


(5.24) 


If we could transform the expression on the right-hand side into 
a divergence, we could turn the volume integral in equation (5.21) 
into a surface integral by Green’s theorem. To do so, we proceed as 
follows: 


d d Ix—xAXdu'] x—x^ d Xdii] 

dx\ /L^^J 

__ 1 2{x-~xy)^Xdu'\ x—xAd'^u'\ {x—x{)‘^Xdhi'\ 

Adding the three equations of this type, we obtain 

2 d lx — 1 1 I — ^’iF d'^u^ 

r2 \ji\]^7^\^\7r\Wy 2, r^~ydxdt\ 

Thus (5.24) becomes 
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and therefore 






c j J J ^ dx\r dx 
c J J r dv I dt \ 


[S]) 


dV 


(5.25) 


Substituting in (5.21) from (5.23) and (5.25), we obtain 

which is Kirchhoff ’s formula in the case when F and all the singu- 
larities of u lie outside S. 

If, however, P lies inside 8 and all the singularities of u lie outside, 
equation (5.21) has to be replaced by 

S a Vi 

where is the volume bounded externally by 8 and internally by 
the small sphere a with centre P and radius €. A repetition of the 
previous argument gives 

S a 

Jf we now make e tend to zero, remembering that djdv ~ —djdr 
on a, we obtain 

s 

which is Kirchhoif’s formula in the case when P lies inside 8 and 
all the singularities lie outside. f 

The proof of Kirchhoff’s formulae in the case when all the singu- 
larities lie inside the closed surface 8 now proceeds as in §5.1. 


§ 5.3. The theorem of determinacy of the solution^ 

The formulae of Kirchhoif enable us to express the value of u at 
a point P on one side of a closed surface 8 and at the instant i, in 
terms of the surface values of u and its first partial derivatives at 

t The minus sign is duo to the difference of definition of the normals v and n. 

X Love, Proc, London Math, 80 c, (2), 1 (1903), 37-62 (42-3). 
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§^0 

previous instants. Actually the data are redundant, inasmuch as a 
knowledge of the surface values of u alone for all values of t is 
sufficient to determine those of dujdt, and is, in fact, sufficient, with 
a knowledge of the initial values of u and dujdt, to determine u 
throughout the whole region of space in which it satisfies the equation 
of wave-motions and the prescribed conditions of continuity. This 
statement constitutes a part of the theorem of determinacy of the 
solution of the equation of wave-motions; it follows from it that, if u 
and dujdt are given everywhere initially and u is given for all values 
of t on the surface S, dujdn can have only one definite value at any 
point of 8 at any given instant. The surface values of dujdn would 
also suffice in the same way for the determination of u\ this consti- 
tutes the other part of the theorem of determinacy. 

To prove this theorem, let 


=*//II(£)+(5F+(IT+?(tF 


dxdydZy 


where integration is over the volume F, bounded by 8y in which the 
wave-function u satisfies the conditions of continuity. Then 


d\V 

dt 


rrr n dud'^u dn dht 

J J J 


dxdydz 





by (h’ecn\s theorem. Tntegrating with respect to t from the initial 
instant ^ = 0, we have 

0 s 

Now if there were two wave-functions satisfying the same initial 
and boundary conditions as well as the prescribed conditions of 
continuity, their difference would be a wave-function u for which 
1^0 vanishes, and, moreover, either dujdn or dujdt would vanish on 
8 for all values of t, according as the surface values of dujdn or n are 
given. It follows that, for such a function, W is zero for positive 

4617 ^ 
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values of t. But since W is the integral of a function which is never 
negative, the integrand 



is zero everywhere in V. From this it follows that u is constant in 
V. Being initially zero, u must then be zero for all positive values of t. 
Hence there are not two wave-functions which satisfy the same initial 
and boundary conditions and the same conditions of continuity. 

To extend the theorem to space outside a closed surface, it is 
necessary to prescribe the asymptotic behaviour of the wave- 
functions under consideration. 


§ 5.4. The physical interpretation of Kirchhoff^s formula 

The velocity potential due to a simple source of strength f(t) is 


“ = -A‘ 


;) - ;[/«)]. 


Hence the velocity potential at a point P outside a closed surface S 
due to a distribution over S of simple sources of strength f{t) i)er 
unit area is ^ ^ , 

«= IJ (5-41) 

where r is the distance from P to a tyih(^al point of S and / may 
depend on the position of the element diS as well as on t. 

Again, the velocity potential due to a doublet of strength F{t) 
whose axis is directed along the unit vector n is 


^_/i\ 

erdn dt ^ ^dn\r)' 


Thus the disturbance at P outside S due to a distribution of doublets, 
directed normally to aV and of strength F(t) per unit area, is speci- 


where djdii means differentiation along the outward normal to /S; 
F may also depend on the position of the element dS, 

It is readily seenf that in crossing a sheet of doublets directed 
normally to 8 of strength F per unit area, there is a discontinuity 


t Cf. Larmor, Proc, London Math, Soc. (2), 1 (1903), 1-13 (6-7). 
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in u of amount AttF and no discontinuity in dujdn', whereas in crossing 
a sheet of simple sources of strength / per unit area, there is a dis- 
continuity —477/ in dujdn and no discontinuity in u. Thus if we wish 
to have a distribution of sources and doublets on S such that the 
disturbance shall change suddenly from a null one just inside 8 to 
that specified by u and dujdn just outside 8, then the distribution 

must consist of sources of strength — ^ ~ per unit area and normally 

directed doublets of strength m/ 477 per unit area; and, by (5.41) and 
(5.42), these produce outside 8 a di.sturbance specified by 



1 c)r^d 
cr dn 


ldtydu\rj^ ^ 


This expression must represent u outside S, which is precisely 
Kirch ho ff’s result. Hence Kirchhoff’s integral formula asserts that 
the. disturbance outside a closed surface S due to real sources inside 8 
is the same as would be produced by a fictitious distribution of sources 

over 8 of strength — — area, together with norirmlly directed 

doublets of strength u/47r per unit area. 

The proof of Kirch ho If ’s formula given in §5.2 and the argument 
of Larmor which we have just sketched both deal with the analytical 
theory of solutions of the equation of wave-motions. In the case 
when the wave-function w is the velocity potential of sound waves 
of small amplitude, the dynamical argument of § 4.4 is applicable 
with slight changes and leads at once to Kirchhoff’s formula. Thus 
Kirchhoff’s formula must be regarded as the analytical formulation 
of Huygens’ principle for sound waves in air.^ 

Kirchhoff’s distribution is not the only distribution of secondary 
sources on the closed surface 8 wliich gives the same effect at a point 
P outside 8 as is given by the primary sources within 8. For we can 
superpose on it any distribution which gives a null effect outside 8, 
such as the Kirchhoff distribution corresponding to any set of 
primary sources outside 8. But the distribution of secondary sources 
obtained in this way does not give a null effect inside 8 and so, if 
this distribution had any physical meaning, an expanding wave -front 
would appear to be propagated not only forwards but also backwards, 
which is impossible. 

Kirchhoff’s distribution is, in fact, the only distribution of 
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secondary sources which gives the correct effect at points outside S 
and a null effect inside. For if it were not, it would be possible to 
construct, by subtraction, a distribution of secondary sources which 
would have a null effect everywhere. 

The analysis of this chapter still leaves open the question of the 
validity of Huygens’ principle in Optics. For, as we have already 
remarked, a theory which regards the propagation of light as 
analogous to that of sound is quite inadequate, since a luminous 
disturbance is characterized by a vector, not a scalar. A very close 
representation of the physical propagation of light was attained by 
Stokes in his memoir on the 'Dynamical Theory of Diffraction’, in 
which the aether was regarded as an elastic solid whose state is 
specified everywhere and at every instant of time by the strain or 
displacement and the velocity of each element of volume: each 
element of volume thus disturbed is regarded as a source from which 
a secondary disturbance is propagated, and the law of the disturbance 
in a secondary wave is found on purely dynamical principles. The 
]*esultant disturbance at any point at a subsequent instant is then 
expressed as a volume integral. Whilst it is true that this volume 
integral can be reduced mathematically to a surface integral, similar 
to that which occurs in Kirchhoff ’s theory, the exact physical repre- 
sentation of the propagation of light according to the elastic solid 
theory is lost in the reduction. 

We shall not refer further to this theory as it has been replaced 
by the electromagnetic theory developed in Chapter III. 


Ex. 1. u is a wave -function whose singularities all lie on one side of 
a plane >S'. Prove that the value of w at a point F on the other sidt^ of S can 
be expressed in the forms 


(i) 


U(P',t) = 


involving only simple sources, and 



(/aV, 


(ii) 




involving only doublets. Show also that neither of these distributions give 
a null effect at a point on the same side of S as the singularities of ii, 
[Consider Kirchhoff ’s integral for P\ the image of P in S.] 


[i(] -- 


Ex. 2. Prove that, if 
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wlioro r is tlio distance from (Xi,yi,Zx), then 

Heiico prove if u is a wave' -function, the value for du/cW. at 

is o})taine(l either by differentiating the Kirchhoff int(‘gral for u with respi'ct 
to or by forming tlio Kirchhoff integral for dujdx. 

Ex. 3. I^y means of Kirchhoff ’s formula prove that, if k is a solution of 
th('. (equation of wave-motions, 





(IS, 


where S is th(* sphere of centre P and radius ct and ( )<, indicates that the 

exprc'ssion insidt' the brackets is to he civaluattMl at the instant t * 0. 
( Poinear6.) 


Ex. 4 . A spherical (expanding wave is defined by the initial values 
F(n) dn cF'(R) 

V’ 

By means of Poincare’s formula (Ex. 3), show that at a point P at a distance 
from O and at tlie instant t. 


ii 


UR 

1“ 2RR^ct 



Apply this to the case wlam F(R) 0 ('xcept wlien R^ < R < R^, and 
deduce the results of § 3.7. (Love.) 

Ex. 5. Dediico Poisson’s formula 

u(P;t) 

C f 

for a wave-function satisfying the initial conditions 




from Poincare’s formula (Ex. 3). 

§ 6. Wave motions in two dimensions 
§6.1. The equation of cylindrical waves 

A wave-function u which does not involve one of the Cartesian 
coordinates, z say, represents a disturbance which is the same in all 
])lanes perpendicular to Oz. Such a disturbance is usually called 


t This theorem shows that Kirchhoff’s formula is consistent. It is of importance 
in electromagnetic theory, where we may apply Kirchhoff’s formula either to the 
components of the electric and magnetic vectors or to the scalar and vector potentials 
and then use the formulae 

d — —ia — grad<^, h — curia; 

c 

by the theorem, both methods give the same results. 
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a cylindrical wave, since the wave-fronts are cylinders with generators 
parallel to Oz\ the corresponding wave-function satisfies the equation 
of cylindrical waves 

This differential equation is also called the equation of two-dimen- 
sional wave-motions, since it occurs in problems, such as the problem 
of the normal vibrations of a membrane, which are actually two- 
dimensional. 

Let us consider the solution of the Avave-eqiiation 
bH . dH , bht 1 

in the case where the initial conditions 

u = g{x,y), ('■-=<•) 

do not involve the variable z. If the solution were 

u == ^{x,y,zj), 
involving the variable then 

u — (f>{x,y,z-\-h,l) 

would also be a solution for any value of the constant h\ and this is 
imijossible since there is only one wave-function which satisfies the 
conditions (6.13). Hence a wave-function which satisfies initial 
conditions not involving the variable must be a cylindrical wave- 
function. 

From this it follows that a method of determining a wave-function 
in three dimensions under given conditions can also be used for 
finding cylindrical wave -functions. This method is called by Hada- 
mard the Method of Descent, for it involves descending from three 
spatial dimensions to two. 

There is, however, an important difference between wave-motions 
in two and three dimensions. To see how this difference arises, let 
us consider first the three-dimensional wave-motion under the initial 
conditions 

M = y,z), f{x, y, z) (t = 0), 

bt 

where / and g vanish whenever 

7J2 
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the initial disturbance is thus confined to the region R ^ a. By 
Poisson’s formula 

u{P-,t) = tM{f}+l(tM{g}), 

we see that the dinturbance at a point P at a distance R > a from 
the origin is null until the instant t = {R~a)[c, and that it is again 
null when t > (i2+a)/c. Thus the disturbance at P lasts for a time 
2a jc; in other words, we have a clean-cut disturbance. 

Let us compare this result with the cylindrical wave -motion 
defined by 

^ y), fix, y) (t = 0 ), 

where / and g vanish whenever 

^2 ^ > a^; 

the initial disturbance is thus confined to the cylinder r ^ a. An 
application of Poisson’s formula shows, as before, that at a point P 
at a distance r > a from the axis of z the effect is null until the instant 
t = (r— a)/c. But after this instant the disturbance at P is never 
null; for the sphere with centre P and radius ct always intersects the 
cylinder r — a, no matter how large t may be, and so the mean values 
Avhich appear in Poisson’s formula are not zero. Thus, if the distur- 
bance in a two-dimensional wave-motion is initially confined to the 
region ^ head of the wave is sharply defined; but 

instead of having a sharply defined rear, such a wave-motion pos- 
sesses a ‘tail’. This characteristic of a two-dimensional wave-motion 
is sometimes called diffusion. 

In most cases this residual after-effect would be expected by the 
physicist. For in such problems the sourees'of disturbance are line 
sources extending to infinity in a three-dimensional space: even 
though such line sources may act only for a finite time, the distur- 
bance from each element of a line source travels with a finite 
speed and so the whole disturbance can never pass any given point 
completely. There are, however, some problems, sucli as the problem 
of the vibrating membrane, in which physical intuition would not 
lead us to expect diffusion. 

From the purely mathematical point of view, the phenomenon of 
diffusion is of great interest. As Hadamardf has shown, it also occurs 


t See his Yale Lectures on Cauchy's Problem (1923), 236. 
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in more general partial differential equations and its occurrence is 
intimately associated with the properties of Hadamard’s ‘elementary 
solution’. 


Ex. 1. Provo by moans of Poisson’s formula that Iho cylindrical wavo- 
function which satisfies the initial conditions 


^ ^uldt^f(x,y) (t =■ 0) 

rt 


u{x,y,t) = i J J/(x-+,-cos<?, 

0 0 ^ 
ct In 


;• d$ dr 


Ex. 2.t Prove that a line source of strength f{t) along Oz generates 
expanding cylindrical waves specifi('d by 

t—rjc 

1 r x/o. 

“ stt J 

' 00 

whore provided that/'(0 0 wlum t < —a. 

If the source acts only during the interval 0 < ^ < P, verify tiiat tlie head 
of the wave is propagated cleanly and that there is a ‘tail’ 

T 

“ "" 2rQ{cH^-r^} j 

0 

when ct — r is large compared with cT, 

Ex. 3. Tho fimction u satisfies the differential etjuation 

cihi 1 dhi 
dX^ Ct’^ 


and tho initial conditions u ^ <l)(x), du/a = wlien t 0, where ff> and i/s 
arc non-zero only in the small interval x' c x r x", Discuss the behaviour 
of u at tho point x = a*o as t varies, and show that there is a residual after- 
effect 



x' 


(Hadamard.) 


§ 6.2. Weber’s solution for ‘monochromatic’ cylindrical waves 
Although it is possible to determine cylindrical wave-functions by 
^^PPlyiiig the method of descent, it is usually more convenient to 
discuss the equation of cylindrical waves, 

d'^u 1 8^u 


t See Lamb, Hydrodynamics (Cambridge, 1916), 291. 
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on its own merits as a two-dimensional problem, without introducing 
a third Cartesian coordinate z. All the formulae concerning wave- 
motion in three dimensions have their two-dimensional analogues, 
but the proofs of the two-dimensional formulae are usually much 
more difficult. 

We prove first a formula, due to Weber, which is the analogue of 
Helmholtz’s formula for ‘monochromatic’ wave-functions. Helm- 
holtz’s theory def)ended on the existence of a monochromatic 
spherical wave-function To extend his theory to two 

dimensions, we must find a monochromatic cylindrical wave-function 

qi = 

where is a function of r = alone. Since the equation of 

cylindrical waves in plane polar coordinates is 

it is evident that v is a solution of 


dh) I dv , ^ 

I?+-rdr+^^=<’’ 


and this is Bessel ’s equation of order zero and independent variable kr. 
The first solution of this equation 


^ nlnl 
0 

has no singularity for any finite value of r. It follows that the 
cylindrical wave-function is not the velocity potential due 

to a source*!* at a finite distance, and so is unsuited to our purpose. 
The coin|)lete solution of Bessel’s equation is 


V = AJQ(kr)+BYo(kr), 

where A and B arc arbitrary constants. Here YQ{kr) is Weber’s 
function! defined by 


Y,{kr) = ^log{lkr)J,(kr)-± J 

7T 7T ^ 




n\n\ 


{lkrf». 


j" The reader who regards cylindrical waves as existing in a three-dimensional 
space will understand by the word ‘source’ a uniform distribution of centres of dis- 
turbance along a lino parallel to Oz and extending to intinity in both directions. 

t See Watson, Bessel Functions (Cambridge, 1922), 64, on the various definitions 
of Bessel functions of the second kind. 

4617 w 
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where 


= 1+1+1+...+A- 


and y denotes Euler’s constant. Weber’s function has a singularity 
at the origin neir which it behaves like (2/7r)logr, but has no other 
singularity at a finite distance. The cylindrical wave-function 
is then the velocity potential of a source at the origin, and 
plays the same part in Weber’s theory as did in Helm- 

holtz’s work. 

The most general cylindrical wave-function of period 2TTl(kc) is of 
the form 

where v does not depend on t and is a solution of 

(V2+i> 0. (6.21) 

Now let r be a closed contour bounding the region D in the {x,y) 
plane. If v and w are two functions whose first- and second-order 
partial derivatives are continuous within and on F, Green’s trans- 
formation gives 


/ ~ JJ {vVhc—wVh') dxdif, 


( 6 . 22 ) 


where djdv means differentiation along the normal to F drawn out of 
D, If V and w are both solutions of equation (6.21), this formula 
becomes /• / o ^ \ 


li-i 


In particular, by taking w ~ Fy(i’rj), where denotes the distance 
from F (x^,yi)y we obtain 




ds = 0 


(6.23) 


provided that F and the singularities of v lie outside F. 

If, however, F lies inside F and the singularities of v lie outside, 
equation (6.23) no longer holds since YQ{kr^) has a logarithmic singu- 
larity at P. In this case we apply Green’s transformation to tlie 
region bounded externally by F and internally by a circle g ol 
centre F and radius €. Formula (6.22) becomes 

/ ~ (v^'^w—wVh)) dxdy. 
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Proceeding as before, we obtain 


K 




ds 


a 

= lii/ * 


since the integral round P does not (iepend on e. Now on a the 

dominant parts of Y^{kr^ and — Po{^^i) respectively (2/7r)log€ 

or^ 

and 2/(776); moreover v and its first derivatives are continuous at P . 
It follows at once that 


lini f 

c^O J 


4r(P). 


We have thus proved the following theorem: 

Let V be a solution of the equation 

irhose partial derivatives of the first and second orders are continuous 
within and on a closed curve P; aiid let 


where i\ is the distance from a fixed point P and djdv means differentia- 
tion along the outward normal to P. Then / = 0 or 4r(P) according 
as P lies oidside or inside P. 

This is Weber’sf analogue of Helmholtz’s theorem. There is a 
corresponding theorem in the ease when v has continuous partial 
derivatives of the first and second orders everywhere outside a 
closed curve P. This is most easily obtained by applying the previous 
argument first to the case when v has no singularities in the annulus 
bounded internally by P and externally by a circle whose equa- 
tion is r ~ i?. The value of 


I + / {'’£’'•<‘"''-’"<‘■'■‘>'3 

r s 

t Math. Ann. 1 (1869), 1-36. 
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where djdv means differentiation along the normal out of the annulus, 
then turns out to be 4i;(P) or zero according as P is or is not a point 
of the annulus. The required result follows by making R-^co, pro- 
vided that the integral over S tends to zero; and a sufficient condition 
for this is that v should behave like Y^ikr) for large values of r. The 
details, which involve a use of the asymptotic jwoperties of Y^{kr), are 
left to the reader. 


§ 6.3. Volterra’s analogue of Kirchhoff’s formula 

In Weber’s formula of §6.2 we may evidently replace Y^(z) by any 
other Bessel function of order zero which behaves like (2/7r)log z near 
3 = 0. In particular we may replace it by 

Y^{z)—iJ^{z) = ~iH^^\z), 

where H''^^(z) is a Hankel functionf of order zero: this gives 

4f(P) = ds. 

r 

Now it is known that, when z is real,t 


-iWTO = 

0 1 

Hence, if is a cylindrical wave-function in which v does 

not depend on /, 

Iftp-D— J_ r f gWofl-cC (pHne-ri) 

1 ’ ^ 1 1 

where P lies inside the closed contour F which contains no singu- 
larities of but when P is outside F, the value of the integral is zero. 
If we put r^d = ip, we obtain 


ds 


u(P t) — — r i— f ^ v~ { ^ 


27rJ ie«/J SvJ »‘i)j 


t Watson, Bessel Functions, 73. 


t Ibid. 180 (10). 
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where, if / ^ we write 

dv dx dv dy dv ’ Sv dr I dv ’ 

the variables x, y, and being regarded as independent. By a slight 
change of notation, this becomes 


u{x^,y^,t) 



(6.31) 


The equation (6.31) expresses the ‘monochromatic’ cylindrical 
wave-function u — integral in which there is no explicit 

mention of the period 27rl{kc), and so it holds for any period what- 
ever. Moreover, as it is linear in u, (6.3 1 ) also holds for any cylindrical 
wave-function obtained by adding ‘monochromatic’ solutions of 
different frequencies. But as an arbitrary function of t can be 
expressed as a sum of periodic constituents by means of a Fourier 
series or integral, the equation (6.31) holds for all cylindrical wave- 
functions. This result, which is due to Volterra,*]* is the analogue of 
Kirchhoff’s theorem for two-dimensional wave-motions. 

By a similar argument, we can show that Volterra’s formula holds 
when (^i,yi) lies outside a closed curve F containing all the singu- 
larities of u, provided that v denotes the normal drawn inwards and 
provided that u behaves suitably at infinity. 

The phenomenon of ‘diffusion’, to which attention was drawn in 
§6.1, is evident from Volterra’s formula (6.31 ). For whereas in Kirch- 
hoff’s formula (§ 5.1 ) the value of a three-dimensional wave-function 
at P at the instant t depended on the disturbance at points Q of 
a surface S at the instant t—PQjc, the value of a cylindrical wave- 
function is expressed in terms of the disturbance at points Q of a 
curve r at the instant t—PQjc and all previous instants. 


§ 6.4. Note on the proof of Volterra’s formula 

Whilst it is possible under certain conditions to justify the apj)eal 
made in § 6.3 to the theory of Fourier series and integrals, it is 
most desirable to give a direct proof of the analogue of Kirchhoff’s 
formula. VolterraJ has ])roved directly a more general theorem 
which includes the result of § 6.3 as a particular case. 

t Acta Math, 18 (1894), 161. Voltorra’s formula has been used by Soinnierfeld 
{Zeits. f. Math. ii. Phys. 46 (1901), 11-97 (§9)) in bis discussion of the diffraction of 
X-rays. 

X Rend, dei Lined (6), Ig (1892), 161, 265; Stockholm Lectures (1912). 
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In this more general theorem, which can be regarded as tlic 
analogue of the Riemann-Green theory of one-dimensional wave- 
propagationf, (x^y^t) are interpreted as Cartesian coordinates in a 
three-dimensional space. On a certain surface a we are given the 
values taken by a cylindrical wave-function n and its first partial 
derivatives. J Then 




1 r r c (dll dx du dy 1 du dt\ 

277 J J — t^'^—r^}\dx dv dy dv dt dv] 

<^1 


^ ^ ff (^-^1 dr 1 an 

CTi 


where = {x—Xi)^+(y—yi)^] 


in this formula o-^ is the area cut out of a by the cone r = 
and d/dv denotes differentiation along the normal to a. 

The surface a in this theorem has to satisfy certain conditions, 
which we shall not attempt to discuss. We may, however, remark 
that, in deducing the formula of § 6.3 from the general theorem, we 
take <T to be a cylinder whose generators are parallel to Ot and which 
has the curve F as cross-section. Volterra’s proof of the more general 
theorem just enunciated and a later proof due to Hadamard arc 
both too difficult to give here. 

Quite recently Professor Marcel Riesz has discovered an elegant and 
simple method for dealing with problems of this type. We conclude 
this chapter by showing where the difficulties in the work of Volterra 
and Hadamard lie and how Riesz has been able to avoid these 
difficulties by making use of the theory of analytical continuation 
of a function of a complex variable. For the sake of simplicity we 
restrict ourselves to the initial value problem, Avhen Volterra’s surface 
(7 reduces to a plane per 2 )endicular to the axis of f. 


§ 7. Marcel Riesz's solution of the equation of cylindrical waves 
§ 7.1. An analogy with potential theory 

Let be a solution of Laplace’s equation 


d'^ii 


(7.1 r) 


t For an account of this theory see, for example, Goursat, Courft (V AnahjHc^ 3 
(1923), 137-55. 

t We recall that a knowledge of u and one non -tangential partial derivative 
determines the other two partial derivatives of the first order. 
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which has no singularities within or on a closed surface S, and let 

1 


V = 


where P (^Tq, Zq) is a fixed point inside S. Then if V is the volume 
bounded externally by 8 and internally by a small sphere 8q with 
centre P, we have 


JJJ d^dydz = 0; 


from this it follows by Green’s transformation that 

dv 




»S' + Ko 


where djdv denotes differentiation along the normal out of V. If we 
now make the radius of 8q tend to zero, we obtain 

When this equation, which is a special case of Helmholtz’s formula 
of § 1.2, is applied to gravitation, it states that the potential at a 
])oint P inside 8 due to matter outside 8 is the same as the potential 
due to ‘Green’s Equivalent Layer’ of matter and normal doublets 
on 8, Alternatively we may regard (7.13) as providing the solution 
of the boundary value problem for Laplace’s equation in which the 
values*!* dujdv are given on 8, It is the latter interpretation 

of the Green’s equivalent layer theorem which interests us here. 

If we replace x and y by ix and iy respectively, u becomes a 
solution of the equation of cylindrical waves 

^ . chi d-u bhi 

and V becomes the j)articular solution of this equation 

By analogy with potential theory, we should expect that an applica- 
tion of (Jreen’s transformation to the identity 

{uL{v)—vL(;u)} dxdydz 0 (7.16) 


ijj 


t Tlio boundary valuos of u and dnjdv aro not independent. 
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would enable us to find u in terms of the boundary values of u and 
dujdv on a surface 8, But if we attempt to carry out this process, we 
immediately meet difficulties which do not appear in potential theory. 
In the first ])lace, the function v defined by (7.15) is real only 

< (z-Zof. 

Accordingly we take the volumef V to be bounded by ])art of the 

^ (z— 2o)“ (7.17) 

and by the area which this cone cuts out of 8; then v is real every- 
where in V. 

The more serious difficulty is that, when we choose V in this way 
and apply Green’s transformation to (7.16), we obtain an integral 
of the required form over a part of 8 together with an integral over 
the cone (7.17), and the latter integral is infinite since v and its 
derivatives are infinite on the cone. 

There are two classical ways of avoiding this difficulty. VolterraJ 
replaced v by its integral with respect to Zq] in fact, he wrote 

( . ^^0 =cosh-i- l^“^ol 

in the identity (7.16). This function has no singularity on the cone, 
but has a line of singularities on the axis of the cone. By cutting out 
the singularities on the axis by means of a small coaxial cylinder 
and then applying Green’s transformation, Volterra obtained a 
formula for ^ 

J •u(Xf„y„,Zo)dzo 

in terms of the boundary values of u and du/dv, and readily deduced 
the required expression for u(xQ,yQ,ZQ), The values of u and dujdv 
on the cone do not appear in the solution on account of the properties 
of Volterra ’s function v. 

Hadamard,l| on the other hand, used the identity (7.16) with 


V{(2-2o)2-(x-a:o)2-(?/-2/o)2}’ 


t This should bo compared with the corresponding stop in the Riemann- Green 
theory of hyperbolic equations in two independent variables. 

% Rend, del Lincei (5), Tg (1892), 161, 265. See also Volterra’s Stockholm Lectures 
(1912). 

11 See, for example, his Yale Lectures on Cauchy's Problem (1923). An excellent 
account of his method is given in Courant and Hilbert, Methoden der math. Physik, 2 
(1938), 438 -42. 
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and applied Green’s transformation. He did not try to avoid the 
occurrence of divergent integrals, but developed a new method of 
])icking out the ‘finite part’ of a divergent integral. 

Professor Marcel Rieszt has shown how all the difficulties 
of Hadamard’s method disappear if we introduce an additional 
complex parameter a. The real part of a can be chosen so large that 
Green’s transformation is immediately applicable; the determination 
of the finite part of Hadamard’s integral is then replaced by the 
process of continuing analytically an analytic function of the com- 
plex variable a. We give below an account of Riesz’s method as it 
applies to a simple yet typical problem in the theory of cylindrical 
waves. 


§ 7.2. Integrals of fractional order 

The repeated integral of order 7 i 


X X 1 a"«-i 

J J ... J f{x„) dx„dx,^_-^ ...dx^ 

a a a 

can 1)0 easily transformed into the simple integral 

J ( 7 . 21 ) 

a 

The latter expression has a meaning for non-integer values of n 
provided that R1 n > 0, and so i3rovides what is usually called the 
Riemann-Liouville definition of the integral of f{x) of fractional 
order n. Professor M. Riesz’s method of solving the equation of 
cylindrical waves depends on an extension of this idea to functions 
of more than one variable. 

b^or simi)licity, we shall consider here the following problem 

To find the value at the point and at the instant 0)0/ 

the ^solution of the equation 




^ S'. 0. 


d^u 

dx'^ dy^ 

given the values of u and dujdt everywhere at the instant t ~ 0, 


( 1 . 22 ) 


t So far, Professor Riesz lias only published a very brief aceouiit of his work in 
Compte.srendii8 du congr^s int. des math. (Oslo, 1936), 2, 44-5. We attempt to give hero 
a simple account of his ideas by elementary methods, and refer the reader to Professor 
Riesz for more elegant methods applicable to more general problems. 

t The extension to more general problems is more complicated but involves no 
essentially new ideas. 

4617 T 
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Actually it is more convenient to regard {x,y,t) as rectangular 

Cartesian coordinates: our ])roblem is then to find the value of the 

solution of T / \ 

L{u) =f(x,y,t) 


at a point P {Xf^y^Jo), given the values of u and dulBt on tlie plane 
t = 0. The definition of the fractional integral approj)riate to this 
jrroblem^ is then as follows: 

The odh integral of u is 

■-= 27 Tf(L^l) JJ/ (7.23) 

JJ 

where T = (t-tff-{x-xf)'^~(y-yff 

and D is the volume bounded by the plane t ^ 0 and the cone F 0. 

If we assume merely the integrability of 2 ^, we can easily show that 
the function I^u so defined is an analytic function of the complex 
variable a, regular when Rla >► 1. In the application we have in 
view, we shall make further restrictions on the nature of u and 
shall be able to continue I^u analytically into the wider doman 
Riot > — 1. 

The integral on the right-hand side of (7.23) is, in form, similar to 
the Riemann-Liouville integral (7.21), but the reason why l^u is 
called the ath integral of u lies deeper than mere similarity of form . 
It can be shown that, for quite a wide class of functions u, 

and so the operator P is the inverse of the differential o 2 )eratorJ Lq. 
Moreover, we write /“ rather than /„, since tlie parameter (v obeys 
the law of indices 


Because of these two i)roperties we can regard i“a as being a 
generalization of the Riemann-Liouville integral of fractional order, 
though this aspect of the function does not concern us here. 

Ex. 1. The function F[X,Y, Z), where < Z^, is defijied by the 

equation 

F(X,Y,Z) = JIJ (z^-x^~y^)v{(Z-zf-(X-xf-(Y-yfYdxd!idz, 
t See § 7.G. 

J The suffix in Lq indicates that Xy y, tin. L are replaced by Xq, 2/o» respectively. 
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where U is tiio volume specified by the inequalities 

< 2^ (X-a:)2+(r— 3/)2 < (Z-zY, 0 z ^ Z. 

Prove that F{X, Y,Z) = F{0, 0, a), 

whore a = +^{Z^-X’‘-Y^). 

Hence show that 

F{X,Y,Z) 

' r(2p+2q+5) 

Ex. 2. Uso Ex. 1 to prove that 

2/o» ^o) “ Z/o» ^o)* 

§ 7.3. A transformation of I°^u 

We have already remarked that, if u is integrable, the function l^u 
is an analytic function of a, regular when R1 a > 1. But if we assume 
that u has continuous first partial derivatives, we can continue I^u 
aiialytically into an analytic function, regular when R1 a > — 1 . As 
a first step in this analytical continuation, we need the following 

Lemma. // Rla > 1 and u is continuousy 

I°‘u(Xo,yo,to) = 2rrr(l~l) J/J “(•*’’^’^ 0 — ' dxdtjdz, (7.31) 

V 

where V is the hemispherical voluyne defined by 

zz {x—Xo)^+{y—ijo)^+z^ < z > 0. 

Under the conditions of the lemma, we can wi ite the triple integral 
(7.23) defining as a repeated integral, namely 

^0 

/“it(a;o,2/o,<o) == ^^(1— Tj J j j dxdy, 

0 2: 

where i] is the area in which 

(.t:-Xo)2+(y-;y„)2 

If we make the substitution t ~ t^—r, we obtain 
/«M(a;o,»/o,<o) 

h 

0 i:' 

where S' is the area 


{x—x^Y+iy—y^Y < r^- 
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This equation can be written in the form 


to 

I°^u{xQ,yo,t„) = — J dr j dxdy. 


where 




We now regard {x.y.z) as Cartesian coordinates in a new space. 
The area is then the projection of the hemisphere S specified by 

(x~XQY-\-{y--y^f-\-z^ = r^, z 0 , 

on the plane z ^ 0, and the element of area on this hemisphere is 


(IS = -dxdy. 
z 


Hence 




0 /S 


^<X .2 

-rf - dS 
r 


jva-2 

-r )- — dxdydz, 
r 


which proves the lemma. 


§ 7.4. The analytical continuation of I^u 

We now transform the triple integral (7..‘U) to polar coordinates 
defined by 

X ™ x^-\-r sin d cos y = «iR ^ sin z — r cos 0 

and, for brevity, we write 

u ™ u{xQ-{-rm\d cos<f), y,)+?*siu0sin0, 

Then we have 

ftt 27 r \tt 

I^u{x^,yQ,ty^) --- J \ \ ur^ Hill 0 do, 

0 0 0 ( 741 ) 


the representation as a repeated integral being valid since u is con- 
tinuous. We also make the additional assumptions that the three 
partial derivatives u^, Uy, and are continuous. 
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If K1 o£ > J , we may integrate by parts, to obtain 
j a sin 0 (16 ~ ^ [mcos“- 10J*"+ j M,gcos“-i0 

irr 

1 1 r 

— “^^“(■^o>2/o>^o ■“'■)+ j (Mj.,cos^+'w^,«iii^)>'cos“0 

Oi 1 (X X f 


Hence we liave 


fo 

’‘u(x„,yoJo) -= J u(u:o,t/„,(o—r)r'^-^ dr 


(^^^cos (l>+Uy Hin. ^)r“cos“0 drd6d<f). 
A further integration by parts gives 

1 

r(a+ r((x-f- lj / Z/o» ^0™^’)^“ + 

0 

(?7j^co8 (fy+u^^sin 0)r«cos“e drdddtj). 

We liave thus proved that, if u and its first 'partial derivatives are 
cofitinuous and if Rlcv > 1, Eiesz's fractional integral of order oc can 
be written in the form 

?/q, ^q) 

1 

= f(^^“K.2/o,0)/S+ J u,S.x„,>ut^-r)r<-dr -|- 

0 

/o 27T Xtt 

27rr(a) J J J drdddtf). (7.42) 

0 0 0 

The function /“m, as originally defined, is an analytic function of 
the complex variable a, regular when Rla > 1. The expression on 
the right-hand side of (7.42) is, howeAmr, an analytic function regular 
when Rla > —1, and so this equation provides the analytical 
continuation of I°‘u across the line Rla = 1. When u and its first 
partial derivatives are continuous, we define /“m by equation (7.42) 
in the domain Rla > —1. 



62 


THK ANALYTICAL REPRESENTATION OF 


[chap. I 


In particular, we have 

to 

I'>u{xo,yo,to) = u{xo,yo,^)+ J VipoyyoJo—r) 

0 

= M(a;o.2/o.<o)- 

The trawiforviation i® is then the identical transformation. 
Example. Provo that 


(7.43) 


§ 7.5. Riesz’s solution of the initial value problem 

The problem is to find the value of the solution of 


L(u) ~ 


d^u 

^2 


8^u 8hi 
^2 “^2 




( 7 . 51 ) 


at the ‘point-event’ (a‘o,?/o,fo)> given the values of u and duldt when 
t = 0. We suppose that tg > 0, so that the problem is an initial 
value problem; the case when <0 < 0 can be similarly treated but is 
of less importance. 

As in § 7. ‘2, we tvrite 

r = {t-t,r-ix-x,r-{y-y,r, 

and we denote by D the volume in (x, y, t) sjjace bounded by the 
cone r == 0 and the plane < = 0. If (/, m, n) are the direction cosines 
of the outward normal to S, the boundary of D, we have, by Green’s 
transformation, 


J/J — vL{u)} dxdydt 




(7.52) 


s 

a result which is certainly valid when u and v possess continuous 
partial derivatives of the first and second orders. In this identity 
we take u to be the required solution of (7.51) and 


V = r(«-w/2. 


where Rla > 5 to ensure the validityf of (7.62). Since 
L(v) = a(a-l)ri“-3)/2. 


t The validity of (7.62) could be proved for smaller values of Rla. There is no 
point in doing so as we shall ultimately use analytical continuation. 
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this gives 

{cc{oc - 1 dxdydt 

"" JJ yo)+»^{<-<o)}+ 

N 

+r<“-W4“+m — -71— )] dS. 

\ dx dy 

Now S consists of two parts, a portion of the cone F — 0 and the 
area specified by 

{x-x^f^iy-y^Y < <i, < = 0; 

and the integrand of the surface integral vanishes on the cone. 
Hence 

JJJ dxdydt 

D 

r<“-W2 ^ -(a- 1 dxdy, 

ho 

or, in Riesz’s fractional-integral notation, 

y.„- 

^ (7.53) 

Position in can be specified by polar coordinates (p, <f>), where 
X — a;o+pcos^, y — f/o+psin^ (0 < p < 0 < ^ ^ 2n). 

The given initial values of w and dujdt are then functions of p and <f>, 

u-= lJ{p,<t>), ^=^V(p,<j>) (< = 0 ) 

say. With this notation, (7.53) becomes 

= 2^+1) JJ 

or 

I-u = JJ (<g-p^)<“-«'^Fpdpd^ + 
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By § 7.4, I^u is an analytic function of a, regular when R1 a > — 1, 
since, by hypothesis, ii and its first partial derivatives are continuous. 
The three terms on the right-hand side of (7.54) are evidently also 
regular when Rla > — 1. Hence, by analytical continuation, the 
equation (7.54) which was proved under the assumption Rla > 5 is 
valid when R1 a > — 1. In particular, when a == 0, we have 

72/+^[jJ (il-p^~)-iVdxdy+^Jj {tl-pYiUdxdy , 

*■ i: 5: 

and so, by (7.43), 

J) 

J J ^ ^ ” dxdi/ 4 - 

+ ^ JJ M(^o(^o—(‘^—^of—(!/-!/o)^)-^tixdy. (7.55) 

L 

This is Volterra's formula expressing the value of a solution of L(u) = f 
in terms of tlie values of u and dufdt when t = 0. When / ~ 0, (7.55) 
reduces to a special case of the general theorem enunciated in § (>.4, 
viz. the case when ex is a plane. 

We cannot simplify the last term on the right-hand side of (7.55) 
by differentiating under the sign of integration, since that would give 
rise to a divergent integral; but by a slightly more elaborate pro- 
cedure we can express (7.55) in a form which does not involve the 
differential coefficient of a double integral. Let us write 

so that U is the value of u on the boundary of S. Then equation 
(7.53) becomes 

r I* — 27r , 

+ JJ {C,-IMI-P^r~^nu-V)pdpd<f, +q fUd<l>. (7.50) 
^ 0 
Now suppose that U possesses a continuous partial derivative with 
respect to p ; then 

V—U = —(to—p)U^(T,<f>), 
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where p < t ^ and so the term involving U—U in (7.56) is an 
analytic function of a, regular when Rla > —1. Putting a — 0, we 
obtain 

“(*o>yo>^o) ~ ^ j j j + 

D 

+ ^JJ {il—{^—^of~{y—yoY}-^dxdy — 

S 

^ **’ 0 )^ iu Z/o)^} ^ dxdy -|- 

i: 

2rr 

+±1 v^. 

0 

This transformation of Volterra’s formula is a particular ease of a 
general formula of Hadamard.f 

§ 7.6. The advantages of Riesz’s method 

The first and most obvious advantage of Ricsz/s method is that 
it avoids the awkward limiting processes which appear in the theories 
of Volterra and Hadamard, by introducing an arbitrary complex 
parameter ol and using analytical continuation. 

The second advantage lies much deeper. We remarked in § 6.1 
that there is a striking difierence between the solution of the equation 
of wave-motions in two and in three dimensions; for in three dimen- 
sions wave-propagation is clean-cut, whereas in two dimensions it 
is diffused. More generally, diffusion always occurs in sf)ace of an 
even number of dimensions, but may or may not occur in space of an 
odd number of dimensions. This differeiu’C de])ends ultimately, as 
Hadamard has shown, on the different characters of the elementary 
solution of the equation of wave-motions in spaces with an even or 
odd number of dimensions. Hadamard ’s theory of the finite part of 
a divergent integral can be a])plied in any sjDace with an even number 
of dimensions. But when the number of spatial dimensions is odd, 
Hadamard ’s theory docs not apply; we have to use either Hadamard ’s 
method of descent or else the idea of the logarithmic partj of a 
divergent integral. 

t See Hadamard’s Yale Lectures on Cauchy's Problem (1923), 208, equation (60). 
t See, for example, Courant-Hilbert, Mefhoden der math. Physik, 2 (1938), 443-8, 
for a discussion of the three-dimensional case. 
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Professor Marcel Riesz’s method, on the other hand, does not 
involve any distinction between spaces of an even and odd number 
of dimensions. The solution is, in fact, independent of the number 


of spatial dimensions. 

Let us suppose we wish to find the solution of 

T, ^ _ 3 % dhi dhi 8 % _ r, ,, 


^^m-1 


at (x5,a:§,...,a;®j_^,<®), given the values of u and its ‘conormal’ deri- 
vative on a certain hypersurface 8 in space. The 

function /“« is now defined to be 


I«u{xlxl...,t^) = JJJ u(xi,X2,...,t)T^^-"'^l^dx^dx^...dt, 


where 


m-l 

r = {t-ty- 2 {^r-4)^ 

r-l 


and D is the volume bounded by 8 and the hypcrcone P — 0. The 
constant -H,„(a) is equal to Ja— ^m). 

The function so defined possesses all the properties of the 
simpler function (7.23). It is regular when Rla > m — 2, and can be 
continued analytically into Rla > —1 when u is sufficiently well 
behaved. Riesz then proves the identity 


= /<»+2/- 


f 

(a+‘^) J 


+ 2 - m)/2 


^p(a+2-m)/2’ 


where S is the part cut out of S by the hypercone and v denotes the 
outward ‘conornial’. By analytical continuation, this identity holds 
when a = 0 and provides Riesz’s ‘solution invariantive ’ of the 
generalized equation of wave-motions. 

Finally, we wish to draw attention to the work of Dr. M. Mathisson, 
who has discovered a new method of solving the initial value problem 
for the equation of wave-motions and other equations of normal 
hyperbolic type.f In his most recent work he has been able to give 
a criterion for the existence or non-existence of diffusion without 
having to find the elementary solution of Hadamard. 

Ex. 1. When Rla > 1, the function J'^u is defined by 

= 2jrr(a- l ) JJJ 

I) 

in the notation of § 7.2. Show that, if u and its first partial derivatives are 
I Math. Annalcn, 107 (1933), 400-19; Coinptes rendus, 208 (1939), 1776-8. 
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contimioiis, J^u oan b(> continued into an analytic function regular when 
Rla > —1, and that, when this continuation is carried out, 

J^u{Xq, 2/o» ^o) ~ ^{^ 0 * Vo* ^o)* 

[Express J^u a,?i a power series in Jc,] 

Ex. 2. Prove that 


cosh k'^T 


is a solution of 


"" vr 


o‘u 

'bfi' 


Ex. 3. Prove that the solution of the initial value problem for the equation 

d^u a^u ahi 


+ 


+ 


I ® rr , \ coshtvPo 

S a. JJ - 


v/here is the value of F when ^ = 0. 
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HUYGENS’ PRINCIPLE AND THE DIFFRACTION OF LIGHT 

§ 1 . Kirchhoff’s theory of diffraction 

§1.1. Kirchhoff’s application of Helmholtz’s formula 

According to geonietrictal optics, rays of light are straight lines. 
With this assumption, light from a point-source incident on a non- 
reflecting opaque screen gives a sharply defined shadow, which it is 
convenient to call the geometrical shadow. Actually it is observed 
experimentally that tlie light is propagated up to the screen as if the 
screen were absent, but that, beyond the screen, light enters the 
geometrical shadow. This phenomenon, which violates the laws of 
geometrical optics, is known as the diffraction of light. 

It was Fresnel who first discovered the real cause of diffraction, 
namely the mutual interference of the secondary waves emitted by 
those parts of the original wave-front whicli were not obstructed 
by the diffraction screen. But, in order to obtain satisfactory results, 
Fresnel had to make somewdiat arbitrary assumptions on the nature 
of the secondary Avaves.f Most of the difficulties of Fresnel’s theory 
were overcome by Kirchhoff, who used Helmholtz’s formulation 
of Huygens’ })rincii)le for monochromatic phenomena.*! We give 
here a brief critical account of Kirchhoff ’s work, |1 and refer the reader 
to the standard works on physical optics for a detailed discussion 
of special diffraction problems. 

Before we can apply Helmholtz’s formula to monochromatic dif- 
fraction problems, we have to overcome tw^o rather serious difficulties: 
w^e have to cxf)ress mathematically what is meant by a point-source 
of light, and we have to give mathematical conditions for a screen 
to be opaque and non-reflecting. 

In the first place, the notion of a point-source of monochromatic 
light is a somewhat idealized one. On the electromagnetic theory of 
light it would be reasonable to regard it as a Hertzian electric 
oscillator or as a Fitzgerald magnetic oscillator.*t*'|* Kirchhoff was, 
how^cver, w^orking with the older theory which regarded the aether 

t Ch. I, § 4.1. t Oh. I, § 4.2. 

II For Kirchhoff’s own account of his theory, see his Vorlesungen uber math. Physik, 
2 (Optik), (Leipzig, 1891). 

tt See Whittaker, History of the Theories of Aether and Electricity, 345-7, 360-2. 
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as an elastic solid, and he assumed that the displacement vector in 
the aether is of the formf 

^ ^ik(r~cl) 

r 


when monochromatic light is emitted by a point-source at r = 0, 
A being a constant complex vector. The intensity of the light is then 
measured by A • A/r^, where A is the complex vector conjugate to A. 
He then applied Helmholtz’s formula to each component of the 
displacement vector separately. 

Tollowing KirchhofF, we shall assume in this chapter that mono- 
chromatic light from a point-source L can be represented by one or 
more wave-functions of the form 


Ua 


v^e 


-ikct — 


^ik{r~ct) 


when there is no diffracting screen, r being the distance from L, 
We have to find the corresjjonding wave-function u == when 

diffraction occurs. 

Let us suppose, then, that monochromatic light from a point-source 
L defined in this way is diffracted by an opaque non-reflecting body 
whose boundary is the closed surface S. Let Sq be a small sphere with 
centre L; then if P is a point exterior to 8 and 8q, Helmholtz’s 
formula gives 


4:7rn{P) = 



^ik(ri-cl) 

7*1 dv 


a /eiWi-c-m I 


dS, 


where is the distance from P to a typical point of a surface of 
integration and djdu denotes differentiation along the inward normal. 
This formula would give the value of u everywhere if we knew the 
boundary values of v and dvjdv on 8 and and it is in the determina- 
tion of these boundary values that our second difficulty arises. 

Actually the surface Sq causes no trouble; for it is evident that, 
when the radius of 8^ is small, the screen has no effect at points of 
Sq and so v — dvjdv — dvjdv there. There is, however, a serious 
difficulty in that we know very little about what happens on S; 
roughly speaking, all we can say is that the part of 8 invisible from 

f Actually Kirchhoff used the real part of this complex vector. The use of the 
complex vector simplifies the analysis. 
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L is very feebly illuminated. Kirclilioff made the following assump- 
tions, which, though reasonable, are nevertheless quite arbitrary: 

(i) V ~ Vq, dvjdv = dvjdv on /S^, the part of S visible from L\ 

(ii) V ~ Q, dvjdv “ 0 on /Sg, the part of 8 invisible from L. 

In words, he assumed that there is no change in the light on 8 -^ and 
that 8 ^ is quite dark. In this way, Kirchholf obtained the formula 


^7TU[P) = 

(See Fig. 6.) 



Ti dv 



( 18 . ( 1 . 11 ) 


Before discussing the validity of Kirchholf ’s assumptions we 
introduce a modification of his formula (1.11) which is often more 




convenient to apply. Let us suppose for sim])licity tliat the tangent 
cone from L to 8 has contact along a single connected curve F, which 
divides 8 ^ from 82 * F may be either a simple closed curve or an arc 
which extends to infinity at both ends; for definiteness wc consider 
only the former case, and leave the modifications when F extends 
to infinity to the reader. 

We now construct an unclosed surface S with F as rim with the 
following properties: 

(i) 2 and 8 ^ form a closed surface S'; 

(ii) L lies inside S', P outside. 

(See Fig. 7.) The surface 8 q can be continuously deformed into S' 
without crossing L or P, and this deformation does not alter the 
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value of the integral over aSJ,. Thus in formula (1.11) wo can replace 
Sq by S'. Taking into account the opposite directions of the normal 
V to S and S' on their common part, we find that 


4r7TU(F) 


JJl 


g{k(n-cl) 


dv 




dS, 


( 1 . 12 ) 


where v denotes the normal from the ‘dark’ to the ‘illuminated’ sidef 
of S. 

In the sim])lest type of diffraction problem, light passes through 
a hole in an opaque non-reflecting ])lane screen. The surface S can 
then be taken to be a jfiane area bridging the gap in the screen. The 
extension to more complicated plane scireens can be left to the reader. 


Example. Provo that Kirohhoff’s formula (1.11) can also bo written in the 


form 


where p is tho normal from the illuminated to the dark side of 


§ 1.2. Criticisms of Kirchhoff’s theory 

All we know from experiment is that the part of the screen invisible 
from the source is very slightly illuminated; Kirchhoff assumes it 
to be ])erfectly dark. This assumption makes the boundary values 
of V and dv/^p on the screen discontinuous across the curve F, which 
is certainly not the case physically. It might further be objected 
that Helmholtz’s formula was proved on the basis of continuous 
boundary values, so that it is not evident without further investiga- 
tion that the function w(P) does satisfy Kirchhoff ’s boundary condi- 
tions. In point of fact, it does not. For PoincareJ has shown that 
Kirchhoff s boundary values of v are incompatible with his boundary 
values ofW dvjdv. 

To prove this, let Q be a point within S, and let r^ denote distance 
measured from Q. Then, since and have no singularities 

outside the closed surface 


//-//{ 

V 


gifrrj 



dS = 0. 


( 1 . 21 ) 


t This is the most convenient way of describing the sense of v. Of course, as S is 
not a screen, both sides of are, in fact, illuminated. 

t Thi’orie matMmatiqiie de la liwiUre, 2 (Paris, 1892), 187-8. 

II Wo recall that a knowledge of the boundary values of v suffices to determine 
those of dvjdv, (Ch. I, § 4.3.) 
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Similarly, if is the wave-function satisfying Kirchhoff’s 

boundary conditions, 


^ s% 

But by the boundary conditions, the latter equation reduces to 

V 

Hence (1.21) becomes 






0 . 


( 1 . 22 ) 


This is, however, impossible. For we can construct a diffraction 
problem in which 8^ plays the part of S in the analysis of the pre- 
ceding section and, by (1.12), the right-hand side of (1.22) should 
be, not zero, but the non-zero value of at Q in this new 

problem. 

In spite of this difficulty, Kirchhoff’s formula does give results 
which agree very well with experiment. Some authors explain this 
by saying that the formula gives the first step in an accurate solution 
by successive approximations. The next step would be to obtain 
a second approximation by substituting in Helmholtz’s formula the 
boundary values given by Kirchhoff’s first approximation; and so 
on indefinitely. If this iterative process were a rapidly convergent 
one, the success of the method would be explained; but the iteration 
has never been carried beyond the first stage and the convergence 
never discussed. 

Other authors, notably Kottler, regard Kirchhoff’s formula as an 
accurate solution, not of a boundary value problem, but of a ‘saltus 
problem’. A brief account of this work is given later. 


§ 1.3. The Fresnel-Kirchhoff approximate diffraction formula 

It turns out to be difficult to apply the accurate Kirchhoff formula 
(1.12) to special diffraction problems.. However, in most problems, 
the wave-length A is very small compared with the other distances 
involved, and we can then use a simple approximation which we now 
obtain. 
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If i/ is a point-source of light characterized by the wave-function 
^mr-ci)j^^ Kirchhoff’s diffraction formula can be written in the form 


4tiTu(P) 



7*1 


dS, (1.31) 


where is the distance from i to a typical i)oint Q of S and djdn 
denotes differentiation along the normal n from the illuminated to 
the dark side of 2. (See Fig. 8.) Since 


/3 

^ -- cos 6^1, 
dn 


where 6q and 0i are the angles the normal n makes with LQ and PQ 
respectively, the equation (1.31) may be written as 



Fig. 8 


Now, by hypothesis, the wave-length "Zn/k is small compared with 
Vq and and so we may neglect IjvQ and Ijr^ in comparison with k 
in (1.32). This gives 

1 C C 

u{P) = — ii(cos0i— <;os0o) d/S, (1.33) 

i: 

a result more usually written in the form 

m(P) = :^ f f cos^i)— . (1.34) 

s 

This api^roximation is usually called the Fresnel-Kirchhoff diffraction 
formula. (Cf. Ch. T, § 4.6.) 

As Larmorf has remarked, the formula (1.34) ‘puts in evidence 
the factor of attenuation (ror^)-^ and the phase depending on the 
path Moreover, it is a consequence of Kelvin’s Principle of 

t Proc. London Math. Soc. (2), 19 (1919), 169-80 (174). 

4617 T 
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stationary Phase that, when k is large, the important i)art of the 
integral (1.34) arises from the part of the range of integration near 
whi(4i the phase of is stationary, and so the Fresnel-Kirch- 

hoff formula ‘shows that it is only the elements hS, that lie near the 
line joining tlie source to any point, which xjroduce the disturbance 
at that point’. 


§2. Maggi’s transformation 

§ 2 . 1 . Helmholtz’s integral over an unclosed surface 

We have seen that when Kirchhoff applied Helmholtz’s formula 
to diffraction problems, he was led to an equation (1.12), which 
ex])resses a component of the light-vector as an integral of Helm- 
holtz’s type, extended over a certain unclosed surfa(*e il. We now 
show that such an integral can be expressed as a line integral over P, 
the rim of and carry out the transformation in detail in the two 
most important cases which arise. The resulting expression as a line 
integral provides the most suitable way of apjdying Kirchhoff ’s ideas 
accurately. 

We recall that the value taken at a point P on one side of a closed 
surface aV by a monochromatic wave-function u of period 2TTl{kc) duo 
to sources on the other side of S is given by 


s 


^ikn 

ri 


du_ d 

d7i dn\ 7*1 ]] 


dS, 


where d/dn is differentiation along the normal out of the region con- 
taining P. But if P and the sources are all 
on the same side of S, the value of the 
integral is zero. It follows that, if we deform 
S continuously, the value of the integral is 
unaltered, provided that, in the deformation, 
S docs not cross P or any of the sources. 

Now consider the case when there is only 
one source at the point L, Take any simple 
closed curve P in space and through it draw 
two unclosed surfaces a?j andNg, each having P as rim; we suppose that 
8^ intersects the segment LP but that 8^ does not. Then, by the 
formula of Helmholtz, 



dS. 
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Now keep fixed and continuously deform in such a way that 
it never crosses L or P. Then since the integral over 
altered, the value of the integral over /Sg does not change and so it 
depends only on the form of the rim F and not on the actual shape 
of S^. This means that it must be possible to express the integral 
over S 2 as a line integral round F, say 



^ikri ^ 

U — 

dn dn 


Hi 



(a^ dx dy +a. dz). 


If we denote by n' the unit vector normal to aS\ drawn in the 
opposite sense to n, we have, by Helmholtz’s formula. 



= ~47rw(P)+ J {a^dx +ay dy+a^dz), 
r 


To sum up: if u is a 7nonochromatic ivave-f miction of period "iTrjfkc) 
due to Cl source at L, and if S is an unclosed surface with F as rim, 
then 


>'>’ ' ' 1 ’ ( 2 . 11 ) 

where e = 0 if the segment LP does not cut S, e -- 1 if it cuts S once. 
The normal n is su|)|)()sed to be drawn from the 'illuminated’ to the 
‘dark’ side of S. 

In the difiraction x)roblem which will concern us most, we have 
u — where Vq is distance from the source L. 

We shall find an explicit formula, due to Maggi,‘j' for the vector a 
in this case. We may evidently suppose that the segment LP does not 
cut /S; then, dropping the time factor the formula (2.11) can 

be written in vector notation as 


JJ B*ndAV = I* Q,-tds, (2.12) 

H r 

I The present proof is that of Kottler, Ayinalen der PJiys. 70 (1923), 405-56. For 
a geometrical proof see Rubinowicz, ibid. 53 (1917), 257-78. 



76 


HUYGENS’ PRINCIPLE AND THE 


[chap. II 


where t is the unit vector along the tangent to V and 


pi/cri 


a = a/i+aj+a^k, (2.14) 

If we adopt the usual convention that the directions of t and n are 
connected by the right-hand screw law, an application of Stokes’s 
theorem gives 


JJ B*n = JJ curl a *11 dS, 


s s 

and so a is a solution of curia = B. 


(2.15) 



Suppose that the distance between L and P is 2/. Then we can 
choose our axes so that L is (0, 0, —/), P (0, 0,/). As our problem 
is a ‘ bipolar’ problem, it is most convenient to introduce spheroidal 
coordinates (A, /x, v) defined by the equations 

X ^/sinhAsin/iCosv, y =/sinhAsin/xsinv, 2 : —/ cosh A cos / a. 

The surfaces A = constant arc confocal ellipsoids of revolution with 
L and P as foci, /x, constant confocal hyperboloids of revolution 
with the same foci, and the surfaces v = constant are planes througli 
the line LP . We make the orthogonal curvilinear coordinates 
(A,/Lt, v) quite definite by the conventions 

A ^ 0, 0 ^ 77, 0 V < 277. 

In this system of coordinates we have 

Tq LQ — /(cosh A-f-cos/i), = PQ --/(coshA— cos/x). 


If we denote by r the position vector OQ, drjdX is a tangent vector 
to the curve /x = constant, v = constant. Denoting its magnitude 
by h., we have , . , 

^ ^ = /v(cosh2A— cos V) ; 
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and so we may write 




where is a unit vector tangent to n — constant, v = constant. 
Similarly we define unit vectors and ig by 

Sr _ , . _ I, • 

dfJL^ 

where = />/(cosh3A~ cos^/i), == /sinh A sin [jl. 

The three vectors ig, are mutually orthogonal vectors, though 
they change their directions in space as the point at which they are 
drawn varies. Any vector A at a point Q can be resolved into com- 
ponents A^, A 2 , A 3 in the directions of the vectors ig, and £3 at Q; 
these components are defined by 

A = Aiii-f-A 2 i 2 +A 3 i 3 . 

The differential of r is given by 

AiiidA+Agigdju-f Agigdi/. 

Since the magnitude of dr is ds, the element of length in space, it 
follows that ^ 

By well-known formulae, we now have 

^ Ailri SA\ra/~ro ^Irrjj’ 

and similarly for the other components. The equations (2.15) to 
determine a^, a^, Og are then ■' 

= M.{«^ ip, 

{«-■ -- !(«-)), ,.n, 


since and are independent of v. (2.18) gives 
„ i _ „ 7. _ 


( 2 . 18 ) 



78 


HUYGENS' PRINCIPLE AND THE 


[chap, it 


where ^ is an arbitrary function of A, fx, y. We can, however, take 
(/> to be identically zero; for if not, a would involve an added term 
grad^, which would disappear on integration round P. Hence 
= a 2 = 0. 

To simplify the remaining equations, we observe that, ii 

0(r) 

then 

sin LL — sinh A , sin jjl - -- sinli A ; 

^ d\ dfi d\ diJi 

moreover, (sinh2A+sinV)<A(^o)^('^i) f- 2 /A:/(M)shA 

The equations (2.16), (2.17) then reduce to 

1^(03/13) = /sin 




Integrating, we obtain 


t-Tj) 

aj/tg = /sin V'/'K)'/'(»i) "=/«”! V > 


and so 


sin a 

a = ~~ N ^3’ 


(2.19) 


itilu. <1. — ^3^3 — . X *3’ 

siniiA 7*0 

It remains to put this in a form independent ol our s[)ecial choice 
of coordinates. 

We have 

I Ldx . .du ... dz\ . . , . 

r -iKini^+jcosv. 

h^X dv dy dvj 

Thus ig is perpendicular to the ])lane LPQ, and so is a constant 
multiple of x r^. But since 

To = r+/k, Ti = r-J'k, 

we have 

Tq X = — 2/r xk — 2/“sinh Asin/x(isinv~j cosv), 


2/2 sinh A sin /a 


and so 



§2] DIFFRACTION OF LIGHT 79 

Substituting in the formula for a, we obtain 

Again, it is easily seen that 

2/2sinh2A -- r^r^+rQ-Vj^, 

pikiro^Ti) f* y -p 

and so, finally, a == ® ^ - . (2.110) 

^0^1 r^ri+To-ri 

We have thus proved that, if /S' is an unclosed surface with rim F, 

J J 1 ri 0n\ To I r„ dn\ Jj ~ ] rori+ro-Ti 

'' (2.111) 

provided that the segment LP does not intersect S. The modification 
Avhen LP does intersect /S' is evident. 


Example, p is a fixed unit vector and u — exp(iA;r«p). Prove that 


u = lim ' 


ifr/’ 


when) 

Hence show that 
p\kr, 


/— >00 ^ 0 ®* 

To = r+/p. 


,s 

^ r oxp{ ik(r .p -I- T-i)} ( pxrd’t ^ 

.1 rt n+p-Fi 

r 

provided that tlie vector through P in the din'ction — p does not intersect S, 


§ 2.2. Geometrical optics as a limiting form of physical optics 

There is a general theorem due to Kirchhoff*!* which states that 
geometrical optics is a limiting form of physical optics. More pre- 
cisely, the diffuse boundary of the shadow in diffraction phenomena 
becomes the sharp shadow of geometrical o])tics as the wave-length 
of the light tends to zero. Kirchhoff ’s proof was based on his theory 
of diffraction as outlined in § 1.1, and depended on a transformation 
of Helmholtz’s integral and the use of what would now be called the 
Riemann-Lebesgue lemma. J But, owing to the generality of his work, 

t V orle^umjcn u, math, Phys, 2 {Opiik), 35. See also Encyc. der math, Wissen. 
V 3 . 437. 

X We refer here to the ‘Hilfsatz’ on p. 33 of Kirchhoff 's Vorlesungen. Actually the 
‘Hilfsatz’ is not correctly stated — the additional assumption that dFjd^ satisfies 
‘Dirichlot’s conditions’ is made in the course of its proof. 
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it is not easy to present KirchhofF’s argument in a satisfactory form.t 
In the present section we show how Maggi’s transformation of 
Helmholtz’s integral enables us to prove Kirchhoff’s result under 
conditions which are sufficiently general for most purposes. 

Let us suppose that L is a point-sourcej of light which can, in 
some way, be specified by the wave-function 


This light is diffracted by a non-reflecting opaque body with bound- 
ing surface 8. By (1.12) the effect at a point P is specified by the 
wave-function 


u 


1 r r d -<•'>! 

/ 


dn 


a 


dS, 


s 

where i] is an unclosed surface whose rim is the curve F along which 
the tangent cone from Lto S touches S, and djdn denotes differentia- 
tion along the normal || drawn from the ‘illuminated’ to the ‘dark’ 
side of 



Applying Maggi’s transformation, we obtain 

u = J a-t d.9, (2.21) 

r 

where € = 0 or 1 according as P is or is not in the geometrical shadow. 

f The difficulty is to show that the ‘Hilfsatz’ or the Riemann-Lebesgue lemma 
is applicable. 

J The proof when the light consists of plane waves {L at infinity) follows similar 
lines and is omitted here. 

II The sense of the normal in formula (1.12) has been reversed. 
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To prove KirchhofF’s theorem, we have to show that the line integral 
in (2.21) tends to zero as the wave-length 27r/(A;c) tends to zero, i.e. as 
k -> 00. 

Using the spheroidal coordinates of § 2.1, we have 


t 


dr 

ds 


7 • dA , 7 • M/JLfe , 7 • U/K 


djji 
‘ ds 


dv 


and so a*t 


dv 




7 Wl' 

^ds /(cosh^A— cos^/x) 


pik{2f cosh A -f/) ^ 

d6** 


Hence the integral to be discussed is 


J a*t ds 


/ 


• 2 

g2iA:/coshA ^ 

/(cosh^A — COS^/x) 


( 2 . 22 ) 


r r 

the time factor being omitted. 

We now make the following assumptions: 

(i) P is not on the boundary of the geometrical shadow. This im- 
plies that A > 0 on r, and so the integrand in (2.22) is a continuous 
function of A and /x which can be differentiated as often as we please. 

(ii) On no arc of F is A constant. This means that no arc of F lies 
on a spheroid rg+ri = constant. 

(iii) F is a simple closed curve| on which A, v are continuous 
functions of the -arc s, which can be differentiated as often as we 
please. 

(iv) F can be divided into arcs on which A is strictly monotonic. 
By (iii), dX/ds is of one sign on such an arc and vanishes only at its 
ends. 


Of these assumptions, (i) and (ii) are necessary for the truth of the 
theorem, whilst (iii) and (iv) are conveniently simple sufficient 
conditions. 

Consider an arc y on which cosh A increases steadily from a to b; 
by (i), a > 1. The integral along y is 


/ — J a • t d6’ — J 


^2iA;/coahA 


sin^/x 


dv 


/(cosh^A— cos‘^/x) ds 


ds 


J* gOiw 

a 


sin^/x 


dv ds 7 
:j- du. 


/(ix^— cos‘^/x) ds du 


t The case when F is a simple curve extending to infinity is ruled out here for 
simplicity. The reader will see that the result still holds in this case provided that 
the integral along F is absolutely convergent. 

4617 M 
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where u = cosh A, = 2i/. We write this in the form 


h 



We must next discuss how s depends on u. 

Take any two numbers a' and 6' such that a < a' < 6' <6. 
Then, by the inverse function theorem, s is a steadily increasing 
function of u in the interval a' ^ u ^ b% and can be differentiated 
as often as we i)lease. Hence, by an integration by parts beginning 
with fc' 

i j 

tt' 

we see that the contribution to I of the interval (a', 6') tends to zero 
as n -> 00 . 

To deal with the interval (a, a') we measure s from u ~ a and 
observe that, by (iv), an expansion 

where m is a positive integer greater than unity, holds in an interval 
0 ^ we can choose a' so that a' < reversion of 

series we deduce an expansion 

5 = 2 b^(u—ayl»\ 

r=l 

valid when a ^ u ^ a', the real positive mth root being understood 
here. It then follows that 


(l>{u) 


ds 

du 



u—a 

r=l 


c^{u-ayl'«. 


We have now to consider the behaviour of 


a QQ 

f Y cJu—oY^"^ --- 
J ^ u-a 


as n-> CO, 

We can divide the infinite series into two parts. To the sum of 
those terms for which r > m, our previous argument depending on 



§ 2 ] 


DIFFRACTION OF LIGHT 


83 


integration by parts can be applied, and it turns out that they have 
no effect when n->ao. If r < m, we have 


a: 


yvlm-l 

^inaMv ^ 

jflm 


CO 

^ina r 

j 


pina f y. 

dv = 1 — 

flHm 


on using a well-known integral valid when 0 < r < m. Hence 


n 

J du 


du -> 0 


as a -> 00 . A similar argument shows that the integral over (h', b) 
also tends to zero. Thus, finally. 


lim r a*t ds ~ 0 

A :- >00 


for ea(‘h finite arc y of F having the property assumed under 
condition (iv). 

We have thus proved tliat, under quite general conditions, 


u “ et/n— 




as the wave-length 27 t/(A'c) tends to zero. Tn other words, the limit 
of the wave-function u has a sharj) discontinuity on the edge of the 
geometrical shadow; in the shadow there is absolute darkness. 

By our assumption (ii) we ruled out of consideration the case when 
tliere is an arc*, of F on which A has a constant value. This excluded 
case is of some experimental interest. Let us supi^ose that y is an 
arc on which A has the constant value a. Then 



and this certainly does not tend to zero as h -> oo. Thus, when an arc 
of F lies on one of the spheroids rg-f 7*1 " constant, there is always 
illumination at the point P, even if it is in the geometri(*al shadow, 
no matter how small the wave-length is. 

This result agrees with experiment. For Fresnel observed that. 
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when there is a source of light on the axis of a circular screen, there 
is always brightness on the axis beyond the screen, no matter what 
the wave-length may be. 


§ 3. Diffraction by a black half-plane 

§3.1. The diffraction of plane waves of monochromatic light 

The simplest typo of diffraction problem is that in which plane 
monochromatic light is incident on a plane black screen, i.e. on a very 
thin non-reflecting opaque plane surface. In spite of its simplicity 
a problem of this type is of considerable experimental interest: for 
it can be realized by means of a spectrometer in which the coUimator 
and telescope are both focused on infinity and are separated by the 
diffraction screen. We consider here the simplest problem of this 
tyjDe, namely that in which the black screen is of infinite extent and 
is bounded by a straight edge. 

Let us suppose that the incident monochromatic light is specified 
by the wave-function ^ _ g-mx+ct) 


or rather by its real part. The wave-surfaces are all perpendicular 
to Ox and are travelling with velocity — c parallel to Ox. If P is 
a point at which .rcosa-f-ysina < 0, where — < a < .Iv, Helm- 
holtz’s formula gives 


u(P) = 


1 r r I d l^\ ^1 

s 


dS, 


(3.11) 


where integration is over the whole plane S whose equation is 

a:cosa-l-ysin(\ — 0 


and n is the unit vector — i cos«— j sin a normal to S. 

If, however, a black screen covers up the part of S on which y is 
negative, the wave-function U specifying the effect beyond 8 is 
obtained, according to Kirchhoff’s a.ssumption, by putting 

u = dujdn = 0 


on the screen. Hence (3.1 1) then becomes 


U(P) = 


■*” J J i n / 

•S'l 


e’*’’* du\ 
dnf 


dS, 


(3.12) 


where 8^ is the part of 8 on which y > 0. (See Fig. 12.) 

We now apply to (3.12) Maggi’s transformation in the form 
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appropriate to plane waves. As the source L is the point at infinity 
on the x-axis, the formula of § 2.1, Ex., gives 


U{P) = — ~ J (3.13) 

r 



when P is in the geometrical shadow; in this formula F is the straight 
edge of the screen (the axis of ;:), t the unit vector tangent to F, and 

>i(ri-i.ri)' 

If P has coordinates we have 


a = 


IXFi. 


when Q ^ point on F; hence 

i.r^ = — Xi, i X Ti = — j(z— 2j. 

But since n~ — icosa— j sina, the right-hand screw rule gives 

(ixri)-t = -y^. 


t ™ k, and so 
Hence (3.13) becomes 

00 

U(P)== — !- -w.rfz. 


(3.14) 


It is more convenient to use, not tlie cartesian coordinates of P, 
but its cylindrical coordinates {p,<f>,Zi) with Oz as axis. This gives 
Xi = pcos<f>, yi = psin<f>, rf = p^+{z—Zj)^, 
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If we substitute these values in (3.14) and then change the variable 
of integration to r, where 


we obtain 


U{P) = 


z — ^i+psinhr, 



^ikpcoshr-ikct ^ 

cosh T-|- cos 0 


(3.15) 


The formula (3.15) holds only when P is in the geometrical shadow. 
If the incident light falls directly on P we must add to the right- 
hand side the term u{P), Hence 


U{P) = j 


^ikp cosh T ikrf 


sin<^ 


cosh T-f- cos 


(It, 


(3.16) 


where e = 0 or 1 according as P is or is not in the geometrical 
shadow. We shall write equation (3.16) in the form 

U (3.17) 

where u* is written for €U, the wave-function according to the laws 
of geometrical optics. The term then represents the effect of 
diffraction. 


§ 3.2. Some analytical transformations and approximations 

As we have just seen, when plane monochromatic light is incident 
on a black half-plane, the wave-function u^, where 

00 

r (3.21) 

277 J cosh T-|- COS 9 

0 

gives the correction which has to be added to the wave-function u* 
of geometrical optics in order to account for the diffraction of light 
into the geometrical shadow. The formula (3.21) is rather difficult 
to apply as it stands. In the present section we obtain two trans- 
formations of the formula which are much easier to handle, and we 
deduce approximate formulae valid when p, the distance from the 
edge of the screen, is either very small or very large compared with 
the wave-length ^TTjk. 

The function is evidently a periodi(j function of </> of period 277. 
We consider then only the range —77 < <^ < 77; the part of this 
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range in which — tt < 0 < oc—^n forms the geometrical shadow. 
The two formulae we shall prove are 


00 + ini 


«« = j e^‘ 

2 it J 


^ik(pGOHh^-rl) 


sin<^ 


and 


cosh ^+cos0 

kp 


dC 




where denotes the Hankel function 


(3.22) 

(3.23) 



To prove (3.22) we consider the contour integral 
r gi/f/3coshf sin0 

J cosh 5 + cos ^ 

c 

taken round the closed contour C of the figure. Since the integrand, 
regarded as a function of the complex variable has all its singu- 
larities on the imaginary axis, the value of the integral is zero. 
Hence 

R R-\-ini ^ R-\-\ni 

r _L r _sin^ I* 

J J cosh^+cos^ J cosh^+cos^ 

OR 0 

(3.24) 

But on the side of C parallel to the imaginary axis, we have 
^ = Jt-j-irj, and so 


jR-hiTri 


I 

i 

J 


kp cosh ^ 


sin<^ 


cosh ^ -f- cos ^ 


d^ 


<C I 


|sin<j&| 


|cosh(i?+i7^) I — |cos (f) I 


in 

, f IsilK^I , 
J sinhi?— |cos^l 



88 


HUYGENS’ PRINCIPLE AND THE 


[chap, h 


Hence J 0 as E-> oo. 

li 

Making R tend to infinity in (3.24), we obtain 

00 ^ 00 I - Itti 

C ^ikpcoxH -.di^ f 

J cosh COS 9 J cosh^+^*GS9 

0 0 

the path of integration in the latter integral being the dotted curve 
in the figure which is asym])totic to the line r] ~ ^tt. Tf we substitute 
in (3.21), we find that 

00 + Ini 

— _L r pik(p cosh C-~cO ^in^ 

277 J cosh^ + cosc^ 

0 

which is equation (3.22). 

The simplest way of proving (3.23) is to show that satisfies the 
differential equation 


+ i cos ~ \ sin 

d{kp) 

We cannot find du^^jd{kp) from (3.21), since differentiation under the 
sign of integration gives rise to an integral which does not converge. 
However, from (3.22) we have 


d{kp) 




00 4 - Ini 

J* pikip cosh C-rt) 
0 


i sin (f) cosh ^ 
cosh ^+cos(^ 


dt 


provided that this integral converges uniformly with respect to p; 
this is certainly the case since the integrand behaves like 

^-fepsinh^ 


where ^ == |+^7ri and | is large and positive. Hence we have 

<x>-\-ini 


du^ 


8{kp) 

Using the formula! 


+ i cos (f)U^ ~ j <f) d^- 

277 J 


0 

oo4-i7ri 


m J 


t Watson, Bessel Functions, 180 (8). 
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we obtain the required differential equation 


(3.25) 


kp 


^ g- ikct 

2,7 ’ 


From thi.8 it follows that 

y' 

But, by (3.21), 

“""(0) = - f e dr^.-- 

Ztt J cosh T + COS 0 i 

0 

j)rovi(lcd that — tt <C ^ tt. Hence vve have 

1 r 1 

== - — I'-') iTTsin^ J e'f™«^.&(,i)(^) (3.23) 

*■0 J 

if — TT < ^ < 77, a formula due to Dr. Erdelyi. 

The ecjuation (3.23) enables us to find very easily an approximate 
formula for valid when p is small com 2 :)ared with the wave-Jengtli, 
i.e. when kp is small. For 

-- Jo{^)+iYo(0 = -log^+0(]) 

77 

as ^ -> 0, and so 

9 * 

r = :^kp\og(kp) + Oik-p). 

J 77 


Hence 


n 


1 

277 


er<^<\<f>—i sin <j> kp log(kp)+0{kp)], (3.2(5) 


when kp is small. 

To find an asymptotic formula for valid when kp is large, we 
go back to (3.22). Making the substitution 

cosh^ = 

we have 


g-iiw /• 

2,7 J 


^{Arpcosli^ 


sillc^ 


cosh ^+cos^ 




^ik{p-cl)-\lTTi 

277 


— oc 

/ 


e-^P« . sin<A 

( 1 + cos (f> -j- iv) ^/( 2 y + i y “) 


— ^£^ik{p-d)-\-\TTi I p^-kpv 


277 


'J 


sinc^ 

( 1 + cos ^ + iv)yj{2v + iv"^) 


dv. 


N 
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the rotation of the path of integration through a right angle being 
justified by Cauchy’s theorem. j* It follows by Watson’s lemma J that 


yli }__cik{p~ct)^riTTi r c-kpv 

277 J (l + cos^)^(2y) 

0 

as kp -> 00 . Hence 


dv 


— - 77 ^ . tan l<j) €}kpHfri-ikcl (3.27) 

2i^^2i7TfCpj 

when p is large compared with the wave-length, provided that 

— 77 <C 77 . 

In conclusion we observe that can be expanded as a Fourier 
series of the formH 


1 

«« .= (3.28) 

n=l 

where the up]^er or lower sign is taken according as <j> is positive or 
negative and ij — 0, J, or i according as 0 < |^| < I 77 , 1^1 ~ i77, 
or ^77 < |<^1 <77. The function is a ‘cut’ Hankel function, 

that is, it is the function 


Hm = 

deprived of terms involving negative powers of 

Although this expansion leads at once to the approximation (3.26), 
its real interest lies in the occurrence of the ‘cut’ Hankel functions 
in a physical problem. The function is a solution of the equation 
of cylindrical waves, and so we should expect the Fourier series for 
to be in terms of the normal functions for cylindrical waves. 
These normal functions are, as is well known. 


cos 


(3.29) 


where = Jn{h)±’^Y„{kp), 

the upper or lower sign being taken according as m = 1 or 2: more- 
over, when kp is large, tliese functions behave Hke 

(,y ‘ co>- 

t The argument consists essentially in considering the integral round the complete 
boundary of the fourth quadrant of the circle |t;| == R, land then making R tend to 
infinity. 

t W^atson, Bessel Functions, 236. 

II Copson and Ferrar, Proc. Edin. Math, Soc, (2), 5 (1938), 159-68. See also Watson, 
ibid. 5 (1938), 174-81 ; Erd61yi, ibid. 6 (1939), 11. 
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SO that (3.29) represents an expanding or contracting cylindrical 
wave-motion according as m = 1 or 2. 

The expansion (3.28) consists, however, of a term representing a 
plane wave-motion together with terms of the form 

h^^\kp)^in ncj) 

which cannot be formed from the normal functions (3.29) and do 
not satisfy the equation of cylindrical waves. A phenomenon of this 
type was first noticed by Whipplef in his work on diffraction by 
a wedge. 


§ 3.3. The diffracted waves 

So far, we have shown that the diffraction of plane monochromatic 
light by a black half-jilane is characterized by the wave-function 

U =z 


where represents the effect according to the laws of geometrical 
optics and is the correction term which takes account of diffrac- 
tion. We now consider the function U in more detail. 

Although 16 * is discontinuous across the geometrical shadow, is 
continuous, since has a discontinuity which compensates the 
discontinuity of For we have 

lim lim 'a* = 0; 

^->7r — 0 TT+O 

also, by (3.23), 


lim tt® = — 

— > 77 — 0 


lim ^ 

77 j 0 


Hence lim U — lim U = 

^-> 77—0 ^-^— 77+0 

which proves that U is continuous across r the boundary of the 
geometrical shadow. There is thus no sudden change from light to 
darkness, but a gradual change. 

Moreover, in the geometrical shadow, U ™ and so the intensity 
of illumination is measured by 


when p is large compared with the wave-length. Since the intensity in 
the incident light has measure unity, the illumination of the geo- 
metrical shadow is very feeble, but there is nowhere absolute darkness 


t Proc, London Math. Soc. (2), 16 (1917), 94-111 (104). 
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since is never zero in the shadow. The intensity gradually 

decreases as we move farther into the shadow. 

When p is large compared with the wave-length 

i tan 
2^J{27Tkp) 

so that the diffracted wave is a non-isotropie cylindrical wave propa- 
gated outwards from the edge of the screen. This means that the 
edge of the screen should appear to he luminous when observed from 
points in the geometrical shadow, a result which can be verified 
experimentally. Outside the shadow the effect is masked by the 
incident light. Actually diffraction by a black half-plane is purely 
an edge effect; for neither u* nor depends on the angle of incidence 
cx and so LI is the same for all black screens having the same edge 
and the same geometrical shadow. 

Lastly, we observe that outside the geometrical shadow the two 
wave-motions specified by and interfere and produce inter- 
ference fringes — the diffraction i)attcrn. But when p is large com- 
pared with the wave-length, the anuditudc of is very small except 
near — tt. Thus the diffraction ])attern is observable only near the 
boundary of the geometrical shadow. 

The actual diffraction pattern (*an be calculated by a (careful 
approximation to near the boundary of the shadow. In sj)ite of 
the defects of the theory, the results agree quite well with the 
experiments. 

§ 3.4. The diffraction of spherical waves by a black half-plane 

The i)roblem we have just discussed is the limiting case of the more 
general problem of the diffraction of monochromatic spherical waves 
by a black lialf-planc, whicb we shall now consider. The incident 
light is generated by a point-source L and is characterized by the 
wave-function , 

u =z 

^0 

where is distance measured from the source. 

We choose axes of coordinates so that L lies on the axis of x and 
the straight edge of the screen lies along the 2 :-axis. Then the screen 
covers the half-])lane a:cos a+y siiia — 0, y 0, where the angle a 
lies between On Kirchho fit’s theory, the effect at a point P 
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beyond the plane of the screen is characterized by the 
function 


V{P) 


J \ '>'1 Snj 


dS, 


wave 


where >Sj is the half-plane a: eos ^ ~ y ^ in this formula 

is the distance measured from P and n is the unit vector 


normal to 


—i cos a—j sin oc 



If we a])ply Maggi’s transformation (§2.1) to Kirchhotf s integral, 
we obtain , ^ 

U{P) ^ I a*td6% (3.41) 

r 

where V is the straight edge of the screen and c = 0 or 1 according 
as P is or is not in the geometrical shadow.^ In this equation 


f* V T* 

a 1 


(3.42) 


t-k, 

^0^1 »o>'i+ro-ri 

Let L and P he at distances and respectively from the 
edge (){ the black screen, so that their coordinates are (/Sq. 0, 0) and 
(pi cos <f>, pj^ sin <f>, Zj). Then 

To = —ip„-\-kz, Tj -= — ipiCos<^-j/)isin^+k(3— Cl), 
since x — y — 0 on F. Hence 


“ PoPi(“OS^ 4 + S(3 — Si), 

To X Fi icpisin<^d-jpp(2— Cl)— jc^i cos^+kpopisiiK^, 

(roxrj-t popisin^, 
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and so (3.41) becomes 


V(P) = J 


giA:(ro+ri -ct) gjj2 ^ 

»’0»‘l+P0/®lCOS<^ + z{z— Zi) 


dz, 

(3.43) 


We now make the substitutions 


2 ! = poSinhTQ, z— 2 i = pisinbrj, 
so that Tq = pflCOshTfl, = pj coshri. 

It follows from (3.43) that 

1 r pik(R~ct) gijY A 

U(P) = eu{P)-— , u 

47r J po Pi cosh To cosh Ti cosh (to+ti)+ cos 0 

(3.44) 

where i? — ^o+^‘i- Lastly, we take a new variable of integration 


T = To + Ti 


sinh“^ — + sinh~^ ^ 

Po Pi 


-^1 


It is easily shown that 

dr R 

dz Po Pi cosh To cosh ti 

and that = Po+Pi+^i+2poPiCosh t. 

Thus, finally, oo 


^ik{R~rl) 

R cosh T+ cos 0 


(3.45) 


This formula, which should be compared with the corresponding 
formula (3.16) for plane waves is due to Rubinowicz.t 
Equation (3.45) is usually written in the form 

U r= 


where u* is the wave-function according to the laws of geometrical 
optics and is the correction term which has to be added to account 
for the diffraction of light into the geometrical shadow. As in the 
case of plane waves, neither u* nor depends on the angle a, and 
so U is the same for all black screens having the same straight edge 
and the same geometrical shadow: more briefly, diffraction by a black 
half-plane is an edge effect. 


t Annalen der Phya. 53 (1917), 2.57. The proof given here is due to Kottler 
(ibid. 70 (1923), 405). 
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§ 3.5. The diffracted wave and Fresnel's integrals 

When monochromatic sj)herical waves are diffracted by a straight 
edge, the diffraction wave is, as we have just proved, characterized 
by the wave-function 

^ikiR-d) gJn ^ 

cosh T + cos 0 ’ 


= -1 f 
277 J li 


(3.51) 


where — pg+Pi+2Ji+2poPiCoshT. 

We now consider j’ the asymptotic behaviour of where and p^ 
are large compared with the wave-length ^Trjk, 

Since k is very large compared with l/p^ and l/p^, the important 
part of the integral (3.51) is due to that part of the range of integra- 
tion near which the phase of is stationary, and so is due to the 
part of the range of integration near r = 0. Now when r is very 
small, we have 

correct to the second order in r, where 

-^1 = V{(Po+Pi)^+^i}* 


It follows that 


u 


B ^ 


2 ^ikiRi-d) 

277 


where 


(Ri-ct) r 

Bi J T 

K^poPi] 


^+r’ 


(3.52) 


— 2 COS Itf). 


We have now to discuss the behaviour of 


m = j p.53) 

0 

for large positive values of q. 

By rotating the path of integration in (3.53) through 45° by means 
of Cauchy’s theorem, we obtain 


F{q) 


<x 

J 






2* 


t The subsequent analysis is a modification of a similar investigation made by 
Whipple {Proc, London Math. Soc. (2), 16 (1917), 103). The main difference is the 
rotation of the path of integration in (3.53) which simplifies the justification of what 
follows. 
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Hence we have 

F{q) = 


CO 



0 


/ 




dr 




00 00 

= j dr j dr) 

0 q 

As this repeated integral is absolutely convergent, we may invert 
the order of integration, to obtain 


CO oo 

F{q) = J drj J e-v'r' dr 

q 0 

00 

^ 4 


Vtt 


where /(|) denotes the complex Fresnel integral 


(3.54) 


i 

Finally, if we substitute for F{q) from (3.54) in (3.52), wo find that 

^ ^ tan li> e~«W^’+l-)/(#), 

zVtt ill 

when the wave-length 'Irrlk is small comj3arcd with the distances 
pQ and Pi of the source and the point of observation from tlie edge 
of the screen. 

The intensity of the diffracted wave is measured by 


\u 


B\2 


^7rlt\ 


\mw 


A qualitative discussion of the variation of intensity can then be 
carried out geometrically by means of the curve (known as Cornu’s 
Spiral), whose parametric equations are 


Z = J 008 (^ 2 ) d^, y = J sin(|2) df; 

f f 
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for the distance from the origin to the point of parameter qifs on this 
curve is \f(qi/f) 1 . 


Ex. 1. Show that the equation (3.45) can bo written in the form 


U(P) - €U{P)- 


2 r Qik(R ct) 


l! 

Ri 


7T J R^-r^ ^{(R^-Rl)(R^~Rl)} 


dR, 


where R{ = (p„ |•pl)“ + z^. Ri = (P(,—piY-\-z\y f PL- 

Ex. 2. is the wave-function of a line of simple sources through the 
point L paralh'l to Oz, Prove that 

wlioro rl^ ■- pl \ p\—2p^pycr>iy<j>. 

Hence show that, if th(^ light from this line-source is diffracted by tlie black 
half-plane .rcosoc-f ^«ina - - 0, ?/ < 0, then 


CT) 

TJ(P) = a,(n-\i ( dr, 

J cosh T -[- COS 

— 00 

where R'^ = po + 2po cosh r. 


§ 3 . 6 . The asymptotic behaviour of the diffracted wave 

An alternative method of obtaining the asymptotic behaviour of 
the diffracted wave when p^ and p^ are large compared with the 
wave-length depends on a transformation of the result of Ex. I 
above. This example shows that 


where 


u 


^ _2 F PoPisinc^ 

ttJ ~R'^-f^\j{(R^-Ii\){R^-Rl)} ’ 

Jii 

— iPo + 

r= PL = Vfpg-fpf— 2popicos^+2j2}. 


(3.61) 


In equation (3.61) we make the substitution E = i?i(l-fT) and 
then apply Cauchy’s theorem to rotate the path of integration 
through a right angle. This gives 




— - PoPi sin (f> e' 


00 

^•/) r . 
J 


(,ikRiT 


(iq_^)2_^2}^{2T+T2}V{^f(l+T)2-i2|} 


PqP^ sin (f) e*A:(J?i-tO+l 7 ri 


X 


GO 




( 1 + m)^—r^},J{2u+iu^}.,J{Rl{ 1 + iu)^—Rl} ' 


4617 


o 
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When the wave-length is small compared with pQ and kB^ is large; 
the asymptotic expansion of can then be obtained by applying 
Watson’s lemma. f The dominant term is 


u 


B 


2poPiSin<^ g/A-(7?i-^/)+l7ri r 

J ^(2u) 

0 

__ Pi ^ ^ikiRi-cinini f ( '^ \ 

7T(Rl-ry(Rl-Bl) aJ \lkRJ 


Hence follows the approximation 


Aj(S‘nkRi poPi) 


(3.62) 


This result does not hold when F is near the boundary of the geo- 
metrical shadow, on which <j> — ±7r. 


§ 4. Kottler’s theory of diffraction by a black screen 

Instead of regarding Kirchhoff ’s solution of the problem of diffrac- 
tion by a black screen as a fairly accurate first approximation to 
a boundary value problem, Kottler J has shown that it is the rigorous 
solution of a somewhat different problem; in fact, it is the solution, 
not of a boundary value problem, but of what we may term a ‘saltus 
problem’. 

Let us consider the case of diffraction by a thin black screen in 
the form of a surface S bounded by a rim F; for simplicity, we take F 
to be a simple closed curve and leave the reader to make the 
appropriate modifications when F either extends to infinity or con- 
sists of several closed curves. Wo suppose that one face, say, of 
this screen is illuminated, and that the other, is dark. We shall 
regard and/S_ as constituting a single degenerate closed surface S' . 

If the undisturbed light is specified by the wave-function then, 
according to Kirchhoff, the effect when the screen is present is 
specified by a wave -function u satisfying the boundary conditions 


u = Uq, 


du __ BUq 
dn dn 


on 


0 on S_. 
dn 

Kottler uses a different definition of blackness; he says that S is a 


t Watson, Bessel Functions^ 236. 


t Annalen der Phys. 70 (1923), 405-56. 
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black screen if u and dujdn are discontinuous across S, the dis- 
continuity being defined by 



the actual values of u and dujdn on the two sides of S' are unknown. 

As an example of this idea, let us suppose that the incident light 
is monochromati(i light due to a simple source at the point L, the 
undisturbed wave-function being 

The wave-function u will be of the form where v is independent 

of t and satisfies the following conditions: 

(i) V is a one-valued function of the coordinates of the point of 
observation P, and satisfies the equation 


" 0 , 

excel )t when P is at L or on S'. 

(ii) V becomes infinite at L, its i3rincipal part near L being Ijr^. 

(iii) V has a discontinuity across the degenerate closed surface S' 
specified by 




u, 


'dvy 


F, 


where U and V are known functions, and n is the normal from 
S^^ to S^. 

(iv) V behaves like on a sphere of large radius r. 

Let B be the volume bounded externally by a sphere S of large 
radius R, and internally by the degenerate closed surface S', the 
small sphere o-q with centre L and radius e„, and the small sphere cr^ 
with centre P and radius Then by the usual aiiplication of Green’s 
transformation, we have 




where n denotes the normal drawn into D. 

We consider separately the four terms on the left-hand side of 
equation (4.01). By hypothesis (iv), the integral over S tends to zero 
as a -> 00 . The integral over o-jl tends to —4:7tv{P) as -> 0, where 
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v{P) denotes the value of v at P. By hypothesis (ii), the integral over 
oTg has the limit, as Cq 0, 

pikLP 

47r = 47rz>u(P). 



Hence 

4,.(P) _ 4«v(i')+ jj £) <*■"'' 

S' 

It remains to consider the integral over the two faces of S\ 

Now 


S' S+ S - 

== r f L h_i 1 ^''']+^,- /"'’‘"A 

J J I Sn+\rij d?i_\ Ti ) ri dn 




dS 


by condition (iii). It follows from (4.02) that 


HP) ^ MP)- 


477 



s 



dS, 


(4.03) 


where n is the normal to 8, drawn from the illuminated to the dark 
side. 
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111 the problem of diffraction by a thin black screen having the 
form 8, Kottler’s hypothesis is that 

U ~ Vq, V = dv^dn. 

Hence, by (4.03), his solution of the problem is 




Equation (4.04) is the formula Kirchhoff gave for diffraction by 
a black screen (§ 1.1, Ex.). But whereas Kirchhoff ’s formula pur- 
ported to give a solution of a boundary value problem in which the 
boundary data turned out to be incompatible, Kottler derived the 
same formula rigorously from a different definition of ‘blackness’. 
It is impossible to give a satisfactory physical definition of a thin 
black screen; Kottler ’s work shows us what analytical definition of 
‘blackness’ gives rise to Kirchhoff ’s formula. 

It should be noticed that it is not possible to apjily Kottler ’s 
definition of ‘blackness’ to a thick screen. A further assumption is 
needed, namely that the illuminated part of a thick screen behaves 
like a thin screen and that the shape of the dark j)art is of no 
importance. 
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HUYGENS’ PRINCIPLE FOR ELECTROMAGNETIC WAVES 


§ 1 . Huygens’ principle and Maxwell’s equations 
§1.1. The formulae of Larmor and Tedonef 

Kirchhoff’s formula gives an analytic representation of Huygens’ 
principle for a field specified by a single scalar potential such as 
the velocity potential of sound waves in air. But the original and 
still the most important applications of the principle relate to fields 
of radiation which, owing to the property of polarization, cannot be 
represented by a single scalar potential. To deal with radiation, it is 
in fact necessary to have recourse to the electromagnetic theory of 
light. Whilst it is true that Kirchhoff’s formula can be apj^lied to 
each of the components of the electric and magnetic vectors, this 
does not constitute a valid formulation of Huygens’ principle as it 
possesses no physical interpretation. We shall therefore consider the 
analytical formulation of Huygens’ principle for an electromagnetic 
field in free aether, regarded as a single entity. 

Let the electric and magnetic forces be given vectorially by 

it being supposed that d and h are measured in electrostatic and 
electromagnetic units respectively. Then Maxwell’s eejuations of the 
electromagnetic field in free aether are 

ccurlh=-d, divh-=0, \ 

ccurld — ~h, divd — 0, j 


where c is the velocity of light. It follows from these equations that 
each component of d and h satisfies the wave equation 






As we have already remarked, each of these six components can be 
represented by a Kirchhoff integral, dei^ending on the boundary 

t We are greatly indebted to Professor E. T. Whittaker for drawing our attention 
to the importance of the work of Larmor and Tedoue, and for allowing us to make 
free use of his lecture notes in § 1 of this chapter. 
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values of the component under consideration. But it turns out, as 
we shall sec, that the secondary disturbances obtained in this way 
are not solutions of Maxwell’s equations, and so such a representation 
of the electromagnetic field is not satisfactory from the physical point 
of view; each secondary source in Huygens’ principle ought to give 
rise to an electromagnetic wave. We must then consider, not the 
equation of wave-motions satisfied by each component of d and h, 
but Maxwell ’s equations as a whole. 

Consider an electromagnetic field specified by the vectors d and h 
which liave no singularities in the region D bounded by the closed 
surface S. We wish to find a surface distribution of electric and 
magnetic charges, etc., on H which will give rise to the actual electro- 
magnetic field in /) and a null field on the other side of 8. This way 
of approaching the j)roblem was first suggested by Larmor.f 

At a point P on 8 we can resolve d into a tangential component 
d^ and a component d,j^ in the direction of the inward normal unit 
vector n; similarly for h. Our proposed distribution must be such 
that the corresponding components of d and h are discontinuous 
across 8. Now in order to produce a zero field outside and a tan- 
gential component just inside, there must be an electric currentj 
flowing at P ol strength hfj {4:7r); vectorially the current in the 
element d8 at P is 

I = -^n X = -^n X h. (1-12) 

477 477 

This current gives rise to a surface distribution of electric charge 
Avhich proves to be just sufficient to produce the required dis- 
continuity in the normal component of d. 

The discontinuity in d,,^ in crossing 8 P will be produced only 
by the part of 8 very near to P, and this neighbourhood of P can 
be regarded as a plane in the usual way. We take Oz to be in the 
direction of n, Oy to be the direction of h^. Then 

h jhy+kh^, n = k 
so that I = —ihJ(^ 7 T), 

f Proc. London Math. Soc. (2), 1 (1903), 1—13. For an interesting but entirely 
different method of attacking this problem, see Love, Phil. Tram. (A), 197 (1901), 
1-45. 

X Tins follows by the usual typo of argument. The discontinuity arises only from 
the distribution on the element dS at P, and for points near P this element behaves 
like an infinite piano. 
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The increase of charge per unit time on the surface element dxdy 

is then i 07 i 

a dxdy = - — ^dxdy = - — d^dxdy 
47t dx 47TC 

by Maxwell’s equations. Hence the charge per unit area is 

(T = — d. = — d-n. (1.13) 

47rc " 47 tC 

Thus the effect of the current is to produce a charge d.n/( 477 -c) per 
unit area distributed over /S, the charge being measured in electro- 
magnetic units. But a surface charge of density d-n/(47r) in electro- 
static units produces a discontinuity d*n in the normal electric force, 
i.e. it produces the required discontinuity in d,j. 

Since Maxwell’s equations are invariant under the transformation 

d->h, h->— d, 


we require a magnetic current in 6' of density 


K = — --n xd 

47r 

(1.14) 

1 

and a magnetic charge r “ —Ji^n 

(1.15) 

per unit area to produce the required discontinuities in the normal 
component of h and the tangential component of d. 

The field at a ])oint Pq (xo, y^, within S, due to the electric current 
I and electric charge a, is given by the equationsf 

di = -grado^-^A, 

c 


hj — curloA, 


where O and A are the scalar and vector potentials 



(1.161) 

s 

(1.162) 


t Tho operators cutIq and gradg refer to variation of the position of the point 
yo. *o)- 
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where r is the distance from Pq to a tyiiical point P {x, y, z) of S, and 
square brackets denote retarded values. (See Ch. I, § 4.5.) 

Similarly the field at P^, due to the magnetic current K and magnetic 


charge r, is given hy 


1 . 



h,= 

-gradoT-iB, 

c 



dj — 

— curio B, 


wheref 


Si 

(1.171) 


B == 

1 

X 

(1.172) 


s 


The field duo to the electric and magnetic currents and charges 
distributed over S is 

d -- di+da, h h^+h^ 

and so is given by 

d = —gradgO — ^A— curloB, (ld8) 

c 

h = —gradQ^F — -B + curlo A. (1-19) 

c 

Equations (1.18) and (1.19), where O, 'F, A, and B are defined by 
(1.161), (1.162), (1.171), (1.172) as integrals over the surface S, pro- 
vide the required analytical expressions for Huygens’ principle in an 
electromagnetic field. We shall call these formulae the Larmor-Tedone 
formulae. 

It might be expected that, from their mode of derivation, the 
secondary waves emitted by the surface elements according to the 
Larmor-Tedone formulae would be electromagnetic waves satisfying 
Maxwell’s equations, but this is however not the case. The Larmor- 
Tedone formulae suffer from precisely the same defect as Kirchlioff ’s 
formula did. We return to this point in § 2.2. 

§ 1 . 2 . The complex form of the Larmor-Tedone formulae 
If we introduce the complex vector J 

q = d+ih, 

t Professor Max Born calls 'F and B the scalar and vector antipoteniials. 
t Bateman (Electrical and Optical Wave Motiom (Cambridge, 1915), p. 4) points 
nnt that if wo use a complex time factor with a complex electromagnetic vector, 

it is necessary to take q — in order to avoid confusion in equating real and 

imaginary parts. Maxwell’s equations then become 

ccurlq = ipiq, divq == 0. 


4617 


P 
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Maxwell’s equations take the simple form 


ccurlq = iq, divq = 0. 

(1.21) 

The use of this vector enables us to write the 
formulae very compactly; for, if we write 

Larmor-Tedone 

a = (D-f-iT, C = A+iB, 


we have the following theorem : 


Let S be a closed surface within and on which the cornjplex electro- 
magnetic vector ({(x, y, z, t) and its first partial derivatives are continuous. 
Let r be the distance froyn Pq (^o> 2/oj ^o) ^ typical point P of 8 and 

let n be the unit normal vector at P, drawn inwards. Then, if 

Q — gradofl— -C+icurloC, 

c 

(1.22) 

whare - « 

4^0-= JJ 

s 

(1.23) 

47TtC — JJ [“><q]^> 

(1.24) 


the value of Q is q(xQ, Zq, t) or zero according as lies inside or 
outside S. 


§ 1.3. Tedone's proof of the Larmor-Tedone formulae 

The way in which the analytical formulation of Huygens’ principle 
for an electromagnetic field was obtained in § 1.1, though valuable 
as affording insight into the physical meaning of the expressions 
which occur, can hardly be regarded as sufficient from the logical 
point of view. The following analytical proof is based on one due 
to Tedone.f 

Let the direction cosines of n be (A, /x, v). We shall denote by djdx 
a partial differentiation with respect to x as it occurs explicitly, 
ignoring the fact that r depends on x, and by djdx a partial 
differentiation with respect to x taking into account the dependence 
of r on X. 

We shall suppose that d and h and their first partial derivatives 

t Rendiconti del Lincei, 26 (1917), 286-9. 
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are continuous in the volume D bounded by S, and, in the first instance, 
that Pq lies outside /S. Then, by Green’s transformation, we have 


/// "r i 1 1/ ¥“!'*' 


since = (a;— Xo)2+(y-^„)2+(z— Zo)2, 

The X component of the vector G defined by (1.24) is given by 


4«y;, = JI {-v[?,]+,xfc]}^ 


J//'- 


'dxdydz d C C C 

' +¥.JJJ 


. dxdydz 


m 'dq„ dqy dxdydz 
_dy dzj r 




dxdydz 




dxdydz 


after using the complex form of Maxwell’s e,qiiations. Hence 


47riG ~ curLF — F, 


where 


= 1 / 1 ' 


, dxdydz 


(1.31) 


(1.32) 


By a similar argument 


4^Q ^ JJ {A[gJ+M[?„]+v[yj^ 


[divq] 


dxdydz 
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The second term on the right-hand side vanishes by Maxwell’s 
equations, and so 4.0 = div„F. (1.33) 


The equations (1.31), (1.32), (1.33) also hold when is a point of 
Z>, but the proof just given has to be modified since \jr becomes 
infinite at JJ,. In this case we apply Green’s transformation to the 
volume U bounded externally by S and internally by a small sphere 
2 with centre and radius e. We then obtain 


477i(C+G') - curloF'~-F', 

c 


where 


47r(G+a') diVoF', 

4,iC- = jj 


[nxq]_, 

T 


-J/ 


, AS 
[q.n]—, 


[q] 


dxdydz 


The equations (1.31), (1.32), (1.33) follow at once, since 
C' 0, ir -V 0, F' F 

as € 0. 

We now consider the vector 

0 = — gradQQ — -G-f-i curlpC 


(1.34) 


which appears in the Larmor-Tedone formula. By (1.31) and (1.33) 
we have 

4.0 = — gradodiVoF + -curloF-(-.iF+curloCurloF— -curljF 


== (v*— i- ?!-W 

\ c^dty ’ 


since 

Hence 


curl curl X ™ graddiv X— V^X. 


<>= JJJw 


dxdydz 

r 


(1.35) 
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But by the well-known properties of retarded potentials the expres- 
sion on the right-hand side of (1.36) has the value 0 

zero according as Pq is or is not a point of D, This completes the 
analytical j)roof of the Larmor-Tedone formula. 

§1.4. A modification of the Larmor-Tedone formulae 

The case when all the singularities of the electromagnetic field lie 
inside the closed surface S may be discussed by applying to the 
result of § 1.2 an argument which will by now be familiar (cf. Chapter 
I, § 6.1). The result is as follows: 

Let the complex electromagnetic vector q and its first partial derivatives 
be continuous on and outside the closed surface S, and let rq -> 0 a 5 
r 00 . Let n be the unit vector, normal to S, drawn inwards. Then if 

Q = ~grad^fi~-(!j-fi curio C, 
c 

where i-n£l = JJ [q.n]^, 4iTiC = 

s 

the value of Q is — q(^o>^o»^o>0 according as lies 

outside or inside S, 



§2. The failure of the Kirchhoff and the Larmor-Tedone 
formulae 

§2.1. The connexion between the Larmor-Tedone formulae and the 
Kirchhoff formula 

The Larmor-Tedone formula for Qj. is 

+4/1 JJ 

>sr .s 

For definiteness we consider here the case when q is regular inside 
and on S, so that (A,/x,v) are the direction cosines of the inward 
normal. We shall write this formula as 


where /j, /g, /g denote the jiarts of involving A, p, or v 

respectively. Then 




DXq 


JJ« 

s 


^-L_L 



s 





( 2 . 11 ) 



no 


HUYGENS’ PRINCIPLE FOR 


[chap. Ill 


Let djdx denote partial differentiation with respect to x ignoring 
the fact that r depends on x, and let djdx denote partial differentiation 
when we take account of the dependence of r on x. Evidently 


A 

dx 


fc] = 



and so 



1 dr 
cr dx 




dx 



( 2 . 12 ) 


Also 

^ r 1 

1 r . n 1 dr r ^ 

and so 

d [a 

1 dr ^ ^ ^ d l\\ 

dx^lr 



From (2.12) and (2.13) we have 



( 2 . 13 ) 


(2.14) 


Similar formulae hold for differentiation with respect to the other 
variables. 

Applying (2.13) to the first term in and the analogues of (2.14) 
to the second and third terms, we obtain 




[%] 

W 1 

. . 

dz^r W 


+: 


But since == divq — and divq 0, 

8y dz dx 

this reduces to 

(2.15) 

Applying the same argument to 


_0 

dx. 


C C , ^d8 d C C , d C C , .dS 
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But since 
this becomes 


1 3a: dy j 




dS. (2.16) 


111 the same way, we may prove that 

T _ r ^ 1 


r f v(— 

J J ^ . r [ 




If we add equations (2.16), (2.16), (2.17) and denote, as usual, 
differentiation along n by djdriy we find that 

s 

The second term on the right-liand side vanishes by Green’s theorem, 
and so we have reduced the Larmor-Tedone formula for Q^. to 

But this is precisely Kirchhoff ’s integral for the wave-function q^] 
it represents ^TTqJ(XQ,yQ,ZQ,t) or zero according as P hes inside or 
outside S. 

Vectorially, we may write this transformation of the Larmor- 
Tedone formulae in the form 

= (-HS] -c7 "+ 

s 

+ J J (n X \7) X ? dS. 
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The first term is Kirchhoff ’s integral representing the vector 4^q or 
zero, whilst the second term vanishes since S is n closed surface. 
This result is due to E. T. Whittaker. 

§ 2.2. The inadequacy of the Kirchhoff and Larmor-Tedone formulae 
for dealing with diffraction problems 

We have just seen that the formulae of Kirchhoff and of Larmor 
and Tedone are equivalent, in that both express the components of 
the electric and magnetic vectors as integrals over a closed surface; 
the two formulae differ by a quantity which vanishes in virtue of 
Green’s theorem. In the aj>plication to diffraction problems, integra- 
tion is not over a closed surface but over a cap bridging the gap in 
the diffracting body; an application of Stokes’s theorem shows that 
in this case the two formulae give different results. The question 
arises, ‘Which is the correct formula to apply to diffraction prob- 
lems?’ It turns out tliat neither is suitable; for the vectors d and h 
given by either formula when S is not a closed surface do not satisfy 
Maxwell’s equations. 

We demonstrate this by considering the special case when S 
is the half-plane x = 0, y ^ 0, and the given field is that of a 
plane wave, viz. 

dj. = 0, = 0, d. ^ 

A^r=0, A. 0. 

This is the problem of the diffraction of plane monochromatic light 
incident normally on the ‘black’ half-plane x — 0, y < 0. 

In the first instance we use Helmholtz’s formula, the form of 
Kirchhoff ’s integral suitable for dealing with monochromatic pheno- 
mena. The resulting field turns out to be 

= 0, d^= 0, d. ^ U, 

K = K == K 

where U = 

in the notation of Chapter II, equation (3.17). Now U involves t only 
in the presence of the time factor hence, by Maxwell’s equa- 
tions, U should be of the form where ^ is a constant. The 

term is of this form, but is not. Hence the field given by 
Kirchhoff ’s integral does not satisfy Maxwell’s equations. 
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In applying the Larmor-Tedone formulae, we observe that, by 
§ 2.1, the field obtained is of the form 

// 

S 

etc., where the superscript K indicates the field given by Kirchhofl’s 
formula. Using Stokes’s theorem, we liave 

r 

etc., where F is the rim of S, i.e. the axis of z described in the positive 
sense. In our case the formulae reduce to 






.-_i_ f r^i 
' 47r J [ r . 


D = D*' H ~ 
a - jDf , /if. 

The field given by Kirchhoff’s formula we have already found; it 
follows that the Larmor-Tedone formulae give 


where v 



= 0, 

II 

0 

Xt = u*-\-u^. 


V, 

//^ = 

H. = 0, 


1 

f ^ 

CO 

1 ^ikp rosh t -i kci 


Itt 

J r 477 

J 


by the transformation 2 ^j+psiiihr of Chapter II, § 3.1. Using 
the result given by Watson in his Bessel Functions, p. 180 ( 10), we get 

4 


This field can be obtained by adding to the field of geometrical 
optics, namely d = k«*, h = ju* 

which satisfies Maxwell’s equations, the field specified by 
d = ku^, h = iv-\-}u^. 

To prove that the latter field does not satisfy Maxwell’s equations, 
it suffices to consider the two equations 

dz dx c dt ^ dx dy c dt ’ 
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Since and v involve t only in the factor these equations 

reduce to 


du^ •/ 7 ? 

- = iku^, 
dx dx dy 




but these imply that dv/dy — 0, which is certainly not the case. 

Thus neither the Kirchhoff formula nor the Larmor-Tedone 
formulae, applied to the comjionents of an electromagnetic field with 
an unclosed surface of integration, give solutions of Maxwell’s 
equations. 


§ 3. Kottler’s formulation of Huygens’ principle for electro- 
magnetic waves 

It is not surprising that Kirchhoff ’s formula, applied to each of the 
components of d and h, should not give vectors satisfying Maxwell’s 
equations when S is an unclosed surface. For the secondary waves 
emitted by each element of S are not electromagnetic waves; an 
electromagnetic wave is obtained only by the superposition of the 
secondary waves due to all the elements of a closed surface. 

The reason why the Larmor-Tedone formulae also fail is not quite 
so obvious; the example of § 2.2 shows that the secondary waves 
emitted by each surface element are not electromagnetic waves, 
though one would expect them to be in virtue of the physical 
argument by which the formulae were derived. 

We recall that, according to Larmor, each element of S carries 

(i) an electric current-sheet of density 

I = i-n X h, 

4:7T 

(ii) an electric surface charge of density 


(iii) a magnetic current-sheet of density 

K= -Anxd, 

477 


(iv) a magnetic surface charge of density 


T = 



n, 


all in electromagnetic units. 
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Kottlerf has pointed out that, if we wish each element of a surface 
S to emit electromagnetic waves, we must consider, not only the electric 
and magnetic charges and currents on dS, but also, possibly, electric and 
magnetic line charges on the boundary of dS. The eflfects of these line 
charges cancel out when S is a closed surface; but if S is not closed, 
additional terms due to the line charges on the rim of S are added 
to the Larmor-Tedone formula. 

Let us consider then an unclosed surface S with rim F. Unless 
the vector I at every point of F is tangential to F, each element ds 
of F must carry a charge which varies with the time in such a way 
as to provide the charge carried away by the component of the 
current normal to ds. 


At any point O of F, n, the unit vector normal to S, and ds, the 
elementary tangent vector to F, are at right angles, (choose axes so 
that ds is along 0.r., n along Oz. The current-density at 0 has 
components 




1 j 




the former is tangential to T at 0, the latter normal to F. The charge 
carried away from dx in time t is then 

t 

—dx h^dt. 


Hence if the line-density of electric charge on dx is S in electro- 
magnetic units, we have 


I 

^ dx = ~^dx j hj. df, 

t ' 

1 r 

or, vectorially, i] — ds- hdt, 

477 J 


But as 0 was any point of F, this formula gives the line-density of 
electric cliarge at every i)oint of F. Similarly we can show that the 
line-density of magnetic charge is 0 in electromagnetic units, where 

t 

Q ds ~ — ds- d dt. 

477 J 


t Annalcn dcr Phijs. 71 (11)23), 457-508. This paper contains an interest in^;? critical 
discussion ol* the various attempts to provide a rigorous electromagnetic theory of 
diffraction by a black screen of arbitrary form. 
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Naturally these formulae hold no matter how small the surface 8 
may be. In particular, if wo wish each element of a surface to emit 
electromagnetic waves satisfying Maxwell’s equations, we must sup- 
pose it to carry, not only electric and magnetic current-sheets and 
surface charges, but also line charges of electricity and magnetism 
on its boundary. The effects of these line charges cancel when we 
consider the resultant of the electromagnetic waves due to all the 
elements of a closed surface, but they must be retained if we wish 
the secondary sources to behave like real sources of electromagnetic 
disturbances. Moreover, in the theory of diffraction by a black screen, 
we are concerned with integration over an unclosed surface, and the 
line charges on the boundaries of the surface elements are then of 
importance, since they give rise to a non- vanishing integral along the 
rim of the surface. 

The introduction of line charges on the boundary of each surface 
element necessitates a modification of the formulae of § 1.1. The 
field due to the electric current-sheet and the electric charges is now 

= — grad^O— - A, = curl^A, 

c 


where O and A are the scalar and vector potentials 

t—rlc 


<1. = 1 J| [d.nl!^_£ JJ h.ds* (3.01) 

r 


Similarly, the field due to the magnetic current-sheet and magnetic 
charges is 

d,= 


1 


where 


-curloB, hj = — grado'F— - B, 

c 


T = + i /j'd.ds* (3.03) 

s r 



§3] ELECTROMAGNETIC WAVES 117 

The total field is then given by 

d = —graded) — -A— curio B, (3.05) 

c 

h = — gradoT* — -B+curlflA. (3.06) 

c 

We shall refer to the equations (3.01)-(3.06) as Kotiler's formulae. 

We can most simply describe the field specified by Kottler’s 
formulae as being the field obtained by adding to that of Larmor 
and Tedone the terms due to the effect of the line charges on P, viz. 

t—ric 

d = £^grado J J h*ds^, 

1 ' 

I - rlc 

h — — ^-grado J J d*ds — . 

r 

These extra terms vanish when S is a closed surface. 

Example. Provo that the complc'x foi’m of Kottlt'r’s formula is 

q ~ — gradoQ— ic + icurloG, 
c 

whero 

IttO = JJ[q-nl~+«c I" J q-dsy, 

s r 

•tviC -- JJ [nxq]^'. 

s 

§ 4. The diffraction of electromagnetic waves by a black screen 
§4.1. Kottler’s definition of a black screen 

We have seen that the problem of the diffraction by a black screen 
of waves, characterized by a single scalar potential, can be regarded 
either as a boundary value problem or as a saltus problem. Kirch- 
hoff’s formula solves the saltus problem accurately in the case of 
a thin screen, but gives only a first approximation to the solution of 
the boundary-value problem. 

The same difficulty arises in the case of electromagnetic waves, 
and in a somewhat more acute form, since it is very difficult to 
formulate the electromagnetic properties of a black screen. One way 
out of the difficulty is to regard the screen as perfectly reflecting, 
and then neglect the effect of the reflected wave (cf. Chapter IV, 
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§3.1). Alternatively, we can define ^blackness’ in terms of the 
discontinuity across the screen of the vectors d and h. This idea 
has been developed by Kottler (loc. oit.). 

Let us suppose that we have a point-source of light at the point L, 
and that this light is diffracted by a thin black screen in the form of 
a surface S with rim F. We suppose that one face, aS^ say, of the 
screen is illuminated, and that the other, a 9_, is dark. We denote by 
do and ho the electric and magnetic forces when the screen is absent. 
We have to find solutions d and h of Maxwell’s equations, with the 
following properties: 

(i) The vectors d— do and h— ho are regular at L. 

(ii) On [d]^-=do, [hj';^-.ho. 

(iii) At infinity, d and h vanish to an a])propriate order. 

We regard and aV_ as constituting a single degenerate closed 
surface S\ and apply the complex form of Kottler ’s modification of 
the Larmor-Tedone formulae to the volume bounded externally by 
a sphere S of large radius R and internally by S' and a small sphere 
i^o with centre L. Remembering that each surface element carries 
line charges on its boundary, we obtain the result 


where 


q(P) — —gradn — ~ C + curie, 
c 

47rQ = JJ Lq-n]^V J J [q-nj^^ + 


(4.11) 


t - rjr 


+ [J [qo-nl^Vt'« J J qo-dSy, 

r 

= JJ [n><q]^^+ J J [n>cq|^-|- JJ [nxq„]-p 


where n on S means n^. (Hee Fig. 16.) Moreover, we can omit the 
terms depending on S, since they tend to zero as R -> oo, by hypo- 
thesis (iii). An application of the Larmor-Tedone formula to the 
closed surface gives 

q„(P)= -grad O' curl G', (4.12) 

c 


47r£2' = J J [q„.nj 4niC' -• J J [n x qo] 


where 
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q(P) qo(P)— grad(Q— O') — i(C— <!]')4-icurI(G— C'). (4.13) 

c 



Now, by liypothesi.s (i), q—Qo is regular at L. Hence if we make 
the radius of tend to zero, we obtain from (4.13) the equation 


where 



-Co+icurlCo, 

c 

(4.14) 

47rt2o ~ J J 

r r 1 ^^*’ , • 
[qo*n] — 

' T 

1’ 

(4.15) 

s 




r r .ds 

1 [*ixqo]— ; 


(4.16) 




the normal n is drawn from the dark to the illuminated side of S 
and the direction of description of V is related to n by the right-hand 
screw law. Equations (4.14), (4.15), (4.16) solve the problem of 
diffraction by a ‘black’ screen. 

It is often more convenient to use integrals, not over the 
diffraction screen, but over surfaces bridging the gaps in the screen. 
We discuss here the simplest case, when F is a simple closed curve. 
We construct a surface >S\ with F as rim, such that S+S^ is a closed 
surface which encloses L but not P. The direction of the normal 
vector n on is chosen so that n and the direction of description of 
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r are connected by the right-hand screw law; this means that n is 
drawn from the ‘illuminated’ to the ‘dark’ side of >S^^. (See Fig. 17.) 



If we apply the Larnior-Tedone formula to the closed surface 
S-{-Siy bearing in mind the difference in the directions of n on S and 


Si, we obtain 

qo(^) 

= gradfii-f icurlCi, 

c 


where . 

47rQj = I 



47rfCi “ J 

J [nxqo]^- JJ [nxqo]^. 


Substituting for qo(i^) in (1.14:), we obtain 


<i{P) 

= — gradi^ — -C-h^curlC, 
c 

(4.17) 

where ^ 

[ [qo-n]^ + ^c J J qo-ds^, 

V 


47rQ= 1 

Si 

(4.18) 

47riC = 

X 

0 

(4.19) 

which provide the required modification of formula (4.14). 


§ 4.2. The diffraction of plane waves by a black half-plane 



As an examplef we consider the diffraction of plane waves of light 
by a black half-plane. We sup^iose that the half-plane is a; ~ 0, ^ C 0, 
and that that the light is incident normally^ on the screen. 


t For an account of the application of the theory to the diffraction of electro- 
magnetic waves from a Hertzian oscillator, we refer the reader to Kottler’s paper. 

X There is no loss of generality in making the assumption, since, on the present 
theory, diffraction by a black screen is an edge effect. 
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If the incident light is plane-polarized perpendicular to the edge 
of the screen, the electric and magnetic forcics are given by 

(Iq ~ ~ jg-tfcu+fO 

before diffraction — more precisely by the real parts of these ex- 
pressions. Alter diffraction, the field is specified by 

d — dj^-f-d2, h = h]^-j-h2, 

where and are the vectors given by the Larmor-Tedone formulae 
applied to the surface — 0, y ;> 0, and dg and hg are the terms 
due to Kottler’s line charges on the rim F of the screen. 

We have already proved (§ 2.2) that 

where e — 0 or 1 according as P is or is not in the geometrical 
shadow, 

— L r 

477 J cosh T+ cos ^ 


and V = 

4 

We write these expressions in the form 
di -- d^-f-ki^^"^,. hj = 

where d* and h* form tlie field according to the laws of geo- 
metrical optics. 

The terms due to the line charges on F turn out to be 

d2 = 0, 

h2 —"^grad — J//i^^(A:/o)e“*^'^'(icos^-|-j sin^), 


or ho ~ ?e(i cos (/> -f-j sin <f>) 

say. The total effect after diffraction can therefore be written in the 

d = d*+d« h = h*+h» 
where d^ = ka^, 

= i(iJ+w?cos^)-f ^^sin<^). 

To obtain the effect when the incident light is plane -polarized 
parallel to the edge of the screen, so that 


_jg-iAr(d;-fcO^ 


ke -lArU+cO 
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we recall that Maxwell’s equations are invariant under the trans- 
formation d h, h —d. The total effect is in this case therefore 

d = d*+d^, h = h*+h" 

where d^ = — i(«;-f-?/;cos^)— ~ \iu^. 

So far we have been working with complex wave -functions. To 
deal with a real jiroblem, the final step is to take real parts of tliese 
complex functions: but the change from light polarized perpendicular 
to the edge of the screen to light polarized parallel to the edge of the 
screen is still obtained by the transformation d h, h -> ~d. When 
Ave are using real wave-functions the intensity of light is measured 
by It follows that when xdane-polarized light is diffracted 

by a black half-plane, there is no difference in intensity of the 
diffracted fight corresponding to a difference in the planes of polariza- 
tion. In this respect diffraction by a black half-plane differs from 
diffraction by a perfectly refiecting screen, a point to which we return 
in Chapter IV. 


§ 4 . 3 . The behaviour of the diffracted wave at large distances 
To discuss the behaviour of the diffracted fight Avhen p, the distance 
from the edge of the screen, is large compared with the wave-length 
^Trjk, we use the asymptotic formulae 




2yJ{27Tkp) 


tan 




gt'C 

2yJ(27Tkp) 




^ 

2yJ(27Tkp) ' 


where ^ = k{p—ct)-\~l7T, 

The first of these formulae v^as proved in Chapter II, § 3.2, the other 
two are immediate consequences of the asymptotic expansions of the 
Bessel functions. 

It follows that, when the incident fight is plane-polarized perpen- 
dicular to the edge of the screen. 




-k 


tan IJ) .y 
2^(2Trkp)^ ’ 


h« 


(— i sin<^+j cos<^) 


2^J(2TTkp) 
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Thus, at large distances, the electric force in the diffracted light is 
parallel to the edge of the screen, and the magnetic force is in the 
transverse direction perpendicular to the edge of the screen. 

The corresponding result when the incident light is plane -polarized 
parallel to the edge of the screen is obtained by the transformation 
d “> h, h -> — d. 

In both cases the edge of the screen appears luminous, and the 
intensity of the diffracted light is proportional to 


1 

Snkp 


tan2 1^, 


which agrees with the result obtained on the scalar theory in 
Chapter IT. 

The Poynting flux-of-cnergy vector in the diffracted light is 


— d^^ X ~ (i cos sin (f)) ~ 

47t Sirkp 

Hence, at large distances, the flux of energy is directed radially out- 
wards from the edge of the screen. In particular, at points of the 
illuminated face of the screen, the flux of energy in the diffracted 
light is along the s(;reen. Kottler’s definition of ‘blackness’ is thus 
satisfactory in so far as no energy is reflected from the screen. 



IV 

SOMMERFELD’S THEORY OF DIFFRACTION 

§ 1. Sommerfeld’s many-valued wave -functions 
§1.1. Introduction 

In his work on diffraction by a black screen Kirchhoff was attempt- 
ing to solve a boundary- value problem, but he met the difficulty that 
it is impossible to give a correct physical statement of the conditions 
of the problem. He therefore made what seemed to be a reasonable 
assumption, namely that the field on the illuminated part of the 
screen was the same as in tlie incident light and that the field on the 
dark part was null. Unfortunately it turned out that his solution 
did not satisfy his boundary values, which were actually incom- 
patible. The results of his theory agreed quite well with experiment, 
so that most authors regarded Kirchhoff ’s solution as a fairly 
accurate first approximation. 

When we deal with a perfectly reflecting screen, a quite different 
state of affairs arises, since the boundary conditions at the surface 
of the screen are well known from the electromagnetic theory of light. 
At the surface of a perfect reflector the electric force is normal to 
the surface and the magnetic force tangential: and it is quite easy 
to solve such problems as that of the reflection of light from a perfectly 
reflecting plane mirror by constructing solutions of Maxwell’s equa- 
tions which satisfy these boundary conditions. 

Theoretically it should be possible to solve diffraction problems 
relating to perfectly reflecting screens in the same way, but the 
analytical difficulties are considerable. It is only in the case of two- 
dimensional problems that much progress has been made. The 
fundamental idea, due to Sommerfeld, of using many -valued solu- 
tions of Maxwell’s equations has been applied successfully to the 
problems of diffraction by a perfectly reflecting half-plane or wedge. 
In the present chapter we give an account of Sommerfeld’s method, 
and work out as an example the problem of the diffraction of plane 
monochromatic light by a half-plane. 

§ 1.2. The case of an electrostatic field 

As a preliminary to a discussion of the diffraction of electro- 
magnetic waves by a perfectly reflecting screen, the consideration 
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of a similar problem in electrostatics is likely to prove instructive, 
since Maxwell’s equations include the equations of an electrostatic 
field as the particular case when h is zero and d is independent of 
the time t. 

The potential F in a two-dimensional electrostatic problem is a 
solution of Laplace’s equation 


dW dW 

dx^ dy'^ 


( 1 . 21 ) 


which is constant on the surface of a conductor. Certain problems 
in electrostatics can be solved by the method of images. Let us 



consider, for example, the two-dimensional problem in which the 
fieldf is produced by a charge E at where > 0 in the 

])resence of an earthed conductor along the axis of x. This field is 
null in the half-j)lane y <0: but in the region y > 0 the field is the 
same as if the conductor were removed and ap additional charge —E 
placed at the point (a^i, — yj), wdiich is the image of (x^,yi) with 
res])ect to the conductor. A similar method will give the field due 
to a charge placed between two eiirthcd conductors whose equations 
are 0 ~ 0 and 0 = tt/ji in polar coordinates, n being a positive integer. 

’I'his inethod of images fails completely when we try to find the 
two-dimensional potential due to a charge in the presence of an 
earthed condjictor along the radius 0 — 0 or, more generally, that 
due to a charge placed in the angle between tw^o earthed eonductors 
whose equations are 0 = 0 and 0 = Ht7r//i, where w(> 1) and n are 

■j* Striclly wo should dos(*ribo fliis field as being duo to a uniform line-charge 
parallel to an earthed plane conductor. 
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integers without common factor. For the image process introduces 
additional charges into the part of the plane considered. The 
diagrams illustrate the two cases of a charge placed in angles tt/S 
and 27 t/ 3 respectively. 

A similar difficulty arises when we try to discuss the reflection of 
light by a perfectly reflecting half*])lane or wedge: the image process 
introduces additional sources of light into the part of space outside 
the wedge. 

In electrostatics the difficulty may be overcome by means of the 
theory of conformal representation. We recall that, if V is a solution 
of Laplace’s equation (1.21), there exists a conjugate function U such 

t'*^*^* U+iV--^f{z), (1.22) 

where /(; 3 ) is an analytic function of the complex variable z ~ 

If the complex variables z and Z ~ X-{-iY are connected by a relation 

z = <f>{Z), ( 1 - 23 ) 


where <^(Z) is an analytic function of Z, areas in the plane are 
mapped conformally f on the corresponding areas in the Z plane. 
If we apply this transformation to the complex potential U-\~iV we 


obtain 


u+iv^MiZ)}, 


so that V is also an electrostatic potential in the Z plane. Equipoten- 
tial curves are transformed into equipotentials and, in particular, 
earthed conductors into earthed conductors. If F is the potential 
due to a charge E in the presence of a vsystem of earthed conductors 
in the z plane, it can be shown that the transformation (1.23) turns 
V into the potential due to a charge E', not always the same as E, 
at the corresponding point of the Z plane in the presence of the 
transformed system of earthed conductors. J 

As an api)lication of these ideas, we find the potential due to a 
charge E at the point in the plane of the complex variable 

z ~ x-\-iy = re^*, 


there being an earthed conductor along the positive half of the real 
axis. The variable 6, which is usually called the argument of z, is 
undetermined to an additive multiple of 27r; we fix it by requiring 


t A mapping is ‘conformal’ if it leaves unaltored tho angle hotwoen any two 
intersecting curves. 

J See Joans, 2Vte Matheynatical Theory of KlectricUy and Magnetism (Cambridge, 
1915 ), 264 - 83 . 
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that d shall lie between 0 and 27r. The two faces of the conductor 
can then be described as lying along the radii 0 = 0 and d = 27 t, 
Now apply the transformation 

Z (1.24) 

where has its positive value. The restriction on 0 implies that the 
argument of Z lies between 0 and tt. Hence this transformation 
provides a (1, 1) (conformal majjping of the whole z plane, apart from 
the radius 0 = 0 (or 2tt), on the part of the Z plane in which T > 0, 
and the radii 0=0 and 0 = 27r are turned resx^ectively into the 
positive and negative jjarts of the real axis. We have then to find 
the jiotential due to a charge E' at a point Z^ in the region F > 0, 
there being an earthed conductor along the line F = 0. 

Jf the conductor in the Z jjlane were absent the j^^otential would 
be given by ^ ^2iE' \og{Z-Z^y, 

but wlien the conductor is iiresent we have 

F - 0 (F < 0), 

Ij^iV = _ 2iE' log f ~ f 1 (F > 0), 

the extra term being due to the charge —E' at the image point Z^. 

To obtain the potential due to the charge E at Zi in the presence 
of an earthed conductor along 0 = 0, wc ai^xfiy the transformation 
(1 .24). This givesf 

U+iV = -2iE\og (1.25) 

where the arguments (or jihases) of the complex numbers and 
^Zi lie between 0 and tt, and denotes the comjilex number con- 
jugate to ylz^. 

Now consider the function U-{-iV defined in (1.25) as a function of 
the comjfiex variable z, \'z^ and V;:^ being kept fixed. If we relax the 
restriction that the argument of s shall lie between 0 and 2tt, \'z has 
two values whose arguments differ by tt, and so the function £F 
is a two-valued function of z. 

The analytical difficulties of handling a two-valued function may 
be overcome by introducing the idea of a Riemann surface. Take 
two su})erimx)osed coj^ies of the z plane, and cut each of them along 

t That E' = K follows from a eoiiskleratioii of the behaviour of U-\-i V near the 
point z^. 
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the positive part of the real axis. Then bind the lower edge of the 
cut in the upper sheet to the upper edge of the cut in the lower sheet, 
and finally bind together the remaining edges. In this way we have 
constructed what is called a two-sheeted Rieinann surface. If we 
start at any point on the upper sheet and describe a circle about the 
origin, we reach the congruent point in the lower sheet; a second 
circuit about the origin brings us back to our starting-point. 

We distinguish between the two sheets of the Rieniann surface by 
writing 2 ; and requiring that 6 shall lie between 0 and 27 t in 

the upper sheet, between 27r and Itt in the lower sheet. With this 
convention we define ^/z as being equal to the function 

f(r,0) = 

where has its positive value. Then evidently 

f{r,e+27T) = 

so that what we previously regarded as the two values of ^Iz are the 
values taken by /(r, d) at congruent points in the two sheets. The 
function defined in this way is a continuous one-valued function 
of position on the Riemann surface. 

We now return to the function 

V+iV = 

where is the complex number conjugate to With our new 
definition of V 2 ;, U-\-iV is a one-valued function of position on the 
Riemann surface. In our physical jjroblem 0 lies between 0 and 2tt, 
and so z lies in the upper sheet: accordingly we describe the upper 
sheet as physical space, the lower sheet as non-physical sj)ace. 

The function IJ -f-iF has two singular points on the Riemann 
surface, namely z^ and z^* Of these, z^, the point at which the charge 
E is placed, lies in physical space, but does not. For if 

^z^ = (0 < a < tt), 

we have Vzg == 

and so z^ = ^ 

Hence is the point of the lower sheet congruent to 


in the upper sheet. 
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To sum up : in order to obtain by the method of images the potential 
duo to a charge E in the presence of an earthed conductor along the 
line ^ — 0 we have to consider a two -sheeted Riemann surface of 
which only the upper sheet corresponds to xihysical space. The image 
charge —E lies in the lower sheet. 

In a similar manner we can ajiply the method of images to find the 
potential due to a charge placed in the angle between two earthed 
conductors along the lines 0 — 0 and 0 = mTTjn, where m{> 1) and 
n are integers without common factor. Since this angle can be turned 
into an angle Trjn by the transformation 

we shall have to use an m-sheeted Riemann surface on which the 
function is one- valued. 

§ 1.3. The introduction of a ‘Riemann space’ into diffraction problems 

From the analogy with two-dimensional problems in electrostatics 
Sommerfeld saw that, in order to treat a diffraction problem as a 
boundary-value problem, it was necessary to consider many-valued 
solutions of Maxwell’s equations which are single-valued functions 
of X3osition in a ‘ Riemaim space’, this space being an imaginary space 
which bears the same relation to ordinary three-dimensional space 
as a Riemann surface bears to a plane. Such ‘ Riemann spaces ’ must 
not bo confused with the Riemannian s^iaces considered in modern 
differential geometry. 

For suj3X30se we have a screen with rim F. We construct a ‘ Riemann 
space’ out of a numl)er of suxierimposed ordinary three-dimensional 
spaces. The sheets of this ‘ Riemann space ' are joined together along 
the boundary of the geometrical shadow, sa that we pass from one 
sheet to another by crossing the shadow. The rim F x^lays the same 
part as the ]Doint O did in our simpler Riemann surface in § 1.2. 
Physical space near the screen does not all belong to the same sheet 
of the 'Riemann space’, since we can get from one side of the screen 
to the other only by crossing the boundary of the shadow. It follows 
that light incident on the screen leaves physical space and enters 
non-physical space. On the other hand, the reflected light originates 
in non-physical space and can be observed only when it crosses the 
screen into physical siiace. 

On account of the analytical difficulties nearly all the problems 
which have been solved in this way are two-dimensional problems. 

4617 q 
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§ 1.4. Two-dimensional solutions of Maxwell’s equations 

An electromagnetic field is said to be two-dimensional if the effect 
is the same in all planes perpendicular to a given straight line. If we 
take this line to be the axis of 2 , the electric and magnetic vectors 
d and h are independent of z and so Maxwell’s equations reduce to 


-4 

c 
1 


c/,, — 


dxf 
dK 


dx ’ 


-A “ 

dx 




?4 1 ^ — 0 
dx'^ By ' ' ’ 


c 


dd. 


-h 

C dx’ 


-L 

C " 


dd 


M. 

dx 


ddj . 

' Sy’ 


= 0. 

dx By 


(1.41) 


These equations fall into two groups. In the first grou]), namely 


dx 






1 ^ dhy Bliy, 

c " dx By ^ 


0 = 

dx 


dh^ 


only hy, hy, and d, occur; whereas in the second group, 


BJk 

dx 


1 

c 


d 


1 /’ 


(1.42) 




dh._ 
dy~ 

c 


(1.43) 


0 


only u^yy 


d x ' By ^ 

I 

By 


ddy ddy 
dx ' 


dyy and \ occur. 

If hyy hyy uud sutlsfy (1.42), it is evident that the vectors 
d = d^k, h==hyi+hy] 
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satisfy the equations (1.41). Hence the group of equations (1.42) 
specify an electromagnetic disturbance polarized perpendicular to 
Oz. Similarly, if dy, and A. satisfy (1.43), the vectors 


d == h = ^^k, 

sj)ecify an electromagnetic disturbance polarized parallel to Oz. The 
general .solution of (1.41) is obtained by adding solutions of (1.42) 
and (1.43). 

By eliminating h^. and from ( 1 .42) we find that the third equation 
becomes 


and that the fourth ccjuatioii is satisfied identically. Thus satisfies 
the equation of cylindrical waves, and, if we can determine the 
comf)onents hj. and are given by the first two equations of the 
group. Similarly, if wc eliminate and d^, from (1.43), the third 
e(j[uation of the group becomes 


bVi, _ 2 

dif “ ^ dV^ ' 

and the fourth equation is satisfied identi(*nlly. Thus Ik also satisfies 
the equation of cylindrical waves, and a knowledge of determines 
dj. and d^,. 

The problem of finding a two-dimensional electromagnetic field 
reduces in this way to that of solving the equation of cylindrical 
waves, hut the boundary conditions which have to be satisfied at 
the surface of a ])erfect reflector depend on whether the field is 
polarized perpendicular to or parallel to the axis of 2 . 

At the surface of a ])erfectly reflecting 
screen the electric force is in a direction 
normal to the surface, the magnetic force 
in a tangential direction. In a two-dimen- 
sional problem a perfectly reflecting screen 
is a cylinder with generators parallel to 
Oz, It follows that, if the field is polarized 
perpendicular to Oz so that 

dj, ^ d^ ]k 0, Fto. 19 

wc have to find a solution r/, of the equation of cylindrical waves which 
vanishes on the surface of the screen. But if the field is polarized 
parallel to Oz, so that = A,, = 0, the boundary condition is 

dj, cos ifj-{-dy^iin/j = 0 
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in the notation of the figure. By equations (1.43) this reduces to 


. , dh. , dK 
sin^-^— cos^— 
^ dx dy 


that is, 


where djdn denotes differentiation along the normal to the screen. 
Thus our problem is to find a solution of the equation of cylindrical 
waves whose normal derivative vanishes on the screen. 

It will be observed that, when the field is polarized parallel to Oz, 
satisfies the same conditions as the velocity potential of sound 
waves in air in the presence of a rigid boundary . 


§ 1 . 5 . Sommerfeld*s many-valued solutions! of the wave-equation 
We now introduce the jnany-valued wave -functions which Soin- 
mcrfeld used to solve two-dimensional problems of diffraction by a 
wedge or half-plane. If we take the straight edge of the wedge or 
half-plane as axis of z, we have to solve 

^4.^' ^ ^ 

dy^ dt^ 

under the appropriate boundary conditions. Actually it is more 
convenient to use polar coordinates {p,<f>)y in terms of wdiich the 
equation becomes ^ ( 151 ) 

As the problem is two-dimensional we need to consider, not a 
Riemann space’, but only the Riemann surface constructed in the 
following way. We take p superimposed planes each slit along the 
radius ^ tt, where is any given angle between ± 77 ; we then 

bind the edge ^ 0 of each cut to the edge ^ of 

the cut in the plane below, and finally bind the remaining edges of 
the top and bottom sheets. We thus obtain a Riemann surfiice of 
p sheets with the origin as branch -point; wo can pass from one sheet 
to the next by the cross-bridge at — <^'+^ 7 . The diagram shows 
the connexion of the sheets in the case p ^ 4, as seen looking along 
the cross-bridges towards the origin. 


t The relevant papers are Math, Ann, 45 (1894), 203; 47 (1896), 317. Quit, Nach, 
(1894), 338; (1895), 267. Proc. London Math. Hoc. (1), 28 (1897), .395. Zeitschrijt 
f. Math. u. Phys. 46 (1901), 11. See also Caralaw, Proc. London Math. Hoc. (1), 30 
(1809), 121 ; Lamb, ibid. (2), 8 (1910), 422. 
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If we start from any point on this Riemann surface and make p 
complete circuits about the origin, we return to our starting-point. 
In doing so, the polar angle ^ is increased by 2^7r. We specify the 
particular sheet in which the point (p,^) lies by saying that it lies 
in the nth sheet if ^ lies between l)7r. The coordinates 

and specify the same point. 



Fig. 20 


Our first task is to find solutions of the wave-equation which are 
one-valued on this Riemann surface.f We start by considering the 
simplest case p — 1, so that the wave-functions in question are 
periodic in of period 277. The simplest such solution is that which 
represents * plane waves’, viz. 

giA'p cos(^— 

where a is a constant, real or complex. A more general solution is 
obtained by multiplying by a function of a and then integrating with 
respect to a. In particular, 

pikct p pioc 

U ^ - pikpcos{<j>-(x) 

277 J pioc_pi^' 


integrated round any closed contour in the^plane of the complex 
variable a satisfies the wave-equation. If we write a — this 


becomes 


u — 


pikct 

277 






e}^ 


M. 


(1.52) 


The integrand in (1.52) has simple poles at the points 


When the contour L is a simple closed contour about the pole 
^ the solution (1.52) reduces to 

pikpCK^^i^ ikcl^ ( 1 . 53 ) 


f The present discussion is based on that of Carslaw, loc. cit. 
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which represents monochromatic plane- waves. Moreover u remains 
equal to ilq when L is continuously deformed, so long as L remains 
closed and does not cross any of the singularities. In particular, we 
may take L to be the contour in Fig. 21. (Tn the figure it is assumed 
that lies between ±77; if it does not, there is, in any case, one 

pole ^+2^^1yhig between and the same conclusion follows.) 



Fig. 21 


On the line AB we have ^ where ^ varies from a to 

and 0 < a < jS < 77. Since 

cos ^ == cos ^ cosh R-\~i sin ^ sinh R, 
we see that o 


- 277(6^^-!) J ^ 

ix 

P 


as iJ 00. Again, on the line CD, we have ^ where ^ varies 

from 77+/ to 77-f-S and 0 < 8 < y < 77. Hence 


^ " r . 

27r(l-e-«) J 

TT+S 

y 

1 J ^-kpSoinO ^0 


gAr/) sin f sinh H 
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as 8 00 . Similarly, the integrals along EF and OH tend to zero 

as these segments recede to infinity, remaining parallel to the real 
axis. 



It follows from this argument that we can take the contour L to 
be that of Fig. 22. But it is evident that the portion is obtained 
from by translation through a distance 27r and, since the integrand 
has period 2^, the integrals along and cancel. Hence 

piket r pi^ 

2tt J 

For brevity, we shall call the two curves and taken together, 
the path A, 

We next consider the case — 2. The expression 

piket ^ pioclZ 

gjA*pcos((^-a) _ ^ n 54\ 

integrated along any fixed complex path, is a solution of the equa- 
tion of wave-motions, provided that w e can perform the necessary 
differentiations under the sign of integration. It follows that, if A 
is the path just defined, the function 

piket r piCl2 

UM) » - J « (1.55) 

A 

satisfies the wave-equation; for it can be expressed as a sum of two 
expressions of the form (1.54). 
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The wave-function (1.55), regarded as a function of does not 
have period 2tt, For 


'M(P,^+277) = 



^ikpcon^ 

6*^ 




dL 


(1.56) 


where A' is obtained from ^ by a translation to the left through 
a distance 277 (Fig. 2.*]): and as the integrand in (1.55), regarded as a 
function of does not have the period 277, u{p,(f>) and 
are in general unequal. By a similar argument, we can show that 
u(p, (f)) is, nevertheless, periodic in but that its period is 477. Thus 
the wave-function defined by (1.55) is not a one-valued function in 
the ordinary plane, but is a one-valued function on the two-sheeted 
Riemann surface. Moreover, it is finite and continuous for all real 
values of p. 



Keeping the point (p, 0) in the first sheet of our Riemann surface 
we have 


u{p,<f>)+u(p,4>+2iT) = 



A + A' 


fiikpco^l^ 


eiCii 


dc. 
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But if the polygonal path L[ is obtained from by a translation 
to the left through a distance 47r, this equation gives 

pikct /• ^itl2 

«(p.^)+^p,^+2Tr)= - e»fcpco8? 1 at, 

A-\- A' i- 1(1 — LI 

since the integrals along Zr^ and — L'y cancel. Now the only singu- 
larities of the integrand are simple y)oles at the points 

C = (^'--<^)-1-4W77. 

Hence if we deform the contour .4 into a simple closed 
curve L surrounding the point ^ we obtain 

u(p,<f>)+u{p,(f>+2Tr) ~ 

by the calculus of residues, and so 

«(p.0)+m(p.^+27t) = Uo(p,(f,). {1-57) 

This relation connects the values of the function u at corresponding 
points of the two sheets of the Rieniann surface. 



Still keeping the point (p,(f>) in the first sheet, the path A' can be 
deformed into the two-branched path B (Fig. 24) without altering the 
value of the integral (1.5G). Now’ 

|gifcpco8^j _ gArpsin^siuhiy 

T 


4617 
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and sin ^ sinh rj is negative on B, Hence 

u(pf <()-{- 27r) 0 

as /D -> 00 . From (1.67) it follows that, as /o -> oo, 

The same type of argument can be applied to the function 

u= -- e^icpo<mC_ e dCy ( 1 - 68 ) 

'inp j -^)Ip 

A 

where p is any positive integer. This function has the following 
properties: 

(i) it is a solution of the equation of wave^motions which is p-valued 
in the ordinary plane\ 

(ii) on the p-sheeted Riemann surface it is one^valuedy finite^ and 
continuous; 

(iii) as p ao, u 0 in all sheets except the first; 

(iv) on the first sheet, u—Uq Q as p -> oo; 

(v) the sum of the p values taken by u at any point in the first sheet 
and the corresponding points in the other sheets is Uq. 

In the first sheet of the Riemann surface the wave*function (1.58) 
behaves like the wave-function u^ of plane waves, at any rate at large 
distances from the origin. 

The corresponding solution for cylindrical waves expanding from 
a line-source at (p',<^') isf 

A 

where denotes the Bessel function of the third kind J^—iY^, 

§ 1 . 6 . A transformation of the formulae for the case p = 2 
If (p, <f>) is a point in the first sheet, we have 

uM+2n) = j <'•*** 

B 

where B is the two-branched path of Fig. 24. This path can be 
deformed into the two lines R1 $ = —it, R1 $ = — Stt without alter- 

f Cf. Carslaw, loc. cit. 147 et seq. Carslaw also discusses the corresponding 
thi^-dimensional problem with a point-source. 
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ing the value of the integral. Hence if we write $ = —Tr+irj, 
5 — —^TT+irj on these lines, equation (1.61) becomes 


00 

£iHcct ( 

uip,<l>+27T) = 7- 

4 tt J ^te- 




This reduces to 


«(p.^+2rr) 


g-ifepcoshr; 


00 

oikci n 


^-ikpcoshrj 


COS — 117 ) 


COS 


'^cosi((^'— ^+^)) 

^ F g-iifcpco6h^<^Osh Jr; cos 

TT J coshr; + cos(<^'— <^) '* 


Hence we have <^+27r) = Xuq(p,<I>), (1-62) 

where 

CO 

Z = 1 cosh{ 6 '-<{>) f 

^ 2vr r/ J cosh7^+cos(^ — ^) ' 

0 

If we now make the substitution r = sinh ^77, the formula for X 
becomes 


X = ^co 8 M-<)>) j e-«MT’ 


+cos’(<^'-^)/ 2 ) 


T^ + COS^K^' — <^)* 


It follows that 


dX __ 2iA; 

TT 


COS J( 0 ' — ^)e-2<^pco8*((^'-i^)/2 g-2i7:pT* 


the differentiation under the sign of integration being justified by the 
uniform convergence of the resulting integral. Using a well-known 
result, we have 

— = — /(A \ cos i(i'— 

Bp V V-^PJ 


ewi /4 Q 

Bp 


'I;/'-"'"’ 
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By (1.71) we find that 

Z{t) = V7re-^M-®2^ (1.72) 

where . 

§ 2. The diffraction of plane-polarized light by a reflecting 
half-plane 

§ 2.1. Sommerfeld^s solution 

Let us suppose in the first instance that the plane of polarization 
is perpendicular to the edge of the screen. If we choose the axes so 
that the screen is the half-plane a: = 0, y < 0, the incident light is 
sijecified by ^ ^ (2.11) 

where the angle of incidence, may be supposed, without loss of 
generality, to lie between zb 2 ^* We have to find a solution of the 
equation of wave-motions which behaves like when p is large 
compared with the wave-length 27r/^’, and which vanishes on the 
screen. 

At first sight it would appear that the approj)riate solution ought 
to be = Uj where 

u = UQ{p,(f>\ (f)')~UQ(p,(f>; TT— (^'). (2.12) 

But this solution vanishes all over the plane a: = 0 and so solves, 
not the problem of diffraction by a half-plane, but the problem of 
refiection by a perfectly conducting plane. If, however, we replace 
Uq by Sommerfeld’s two- valued wave-function 

T 

u{p,<f>; <f,') = i') J dX, (2.13) 

— 00 

where T = cos \{<f)—(j>')^j{2kp), 

and consider the wave-function 

U(p,<f>] 4>') = Stt— ^'), (2.14) 

then d^= U provides the solution of our problem, it being supposed 
that in physical space ^ lies between —In and ^tt. This solution is 
obtained by applying the method of images to Sommerfeld’s two- 
valued wave-function, just as the simpler solution (2.12) was obtained 
from Uq. There is a slight difference, in that the angle of incidence 
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of the reflected wave must be Stt— not tt— since this wave 
originates in non-physical space. 



Ki(i. 25 

Each of the wave-functions < 56 '), Stt— ^') has its own 

associated Riemann surface, the cross-bridges in physical space being 
(f) ™ 7 r+^' and (f) ^ respectively. These cross-bridges divide 
physical space into three sectors /Sg? shown in the figure, 

Of these and S 2 belong to the upper sheet, to the lower sheet of 
the Riemann surface associated with ^l{p, <f>; but belongs to the 

upper sheet, and and 8^ to the lower sheet of the Riemann surface 
associated with u(pj(f); Stt— ^'). In both cases 8^ and 8^ belong to 
different sheets. This means that the incident light, which is specified 
by u(p,(f>] ^'), at any rate at large distances from the screen, remains 
in the upper sheet on crossing the screen and so disappears from 
physical space; on the other hand, the reflected light specified by 
~u(p,<f>; Stt— <^') originates in non-physical space and enters physical 
space by crossing the screen into the sector 8^, There is, however, 
no discontinuity of U{p,<l>; <f>') in physical space except across the 
screen. 

We have, however, somewhat anticipated the fact that U(p,<f>; </►') 
does provide the solution of our diffraction problem. In the first 
place, is a solution of the equation of wave-motions and evidently 
vanishes on the two sides of the screen, viz. the radii <!> = and 

i> = 
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Secondly, when p is large compared with the wave-length 27rjk and 
{py(f>) lies in the sector /Sg* follows from (1.71) and (1.72) that 

U(p,^\ 0 ') = Uo{p,(f>; (l)')+0(kp)~K 

Thus U(p,<f>\ <f)') gives a field wliich is practically unaffected by tlie 
presence of the screen at points in the se(?tor iSg large distances 
from the origin. 

Next, in the geometrical shadow we have 
U{p,<f>^<j>')=^0{kp)-^ 

when p is large compared with the wave-length. Thus the geometrical 
shadow is practically dark. 

Finally, in the sector 

U{p,<f>; <l>') = <f>')—Uo{p,<^; Tr—<f>')+0{Jcp)-i, 

where p is large compared with 27r/k, Hence in the sector we have, 
in effect, the ordinary field of equation ( 2 . 12 ) above due to the 
incident and reflected light. 

The function U{p,(j>\^') does satisfy all the conditions of the 
diffraction problem. A closer examination of the order-terms 0(kp)-^ 
will give the diffracted light. 

If, however, the plane of polarization of the incident light is 
parallel to the edge of the screen, the original field is 

K = i>')y 

where is the angle of incidence, supposed to lie between i in. 
The diffraction problem is solved by taking — K, where V is a 
wave-function which behaves like Uq when p is large comj)ared with 
the wave-length and which satisfies the boundary condition dV Jdcf} — 0 
on the faces <f) = —in and <f> ~ ^n of the screen. An argument 
similar to that used above shows that 

Vipy<f>\ i>')+'f^{pyi>\ Sn—(f>') (2.15) 

satisfies these conditions. When kp is large, we have 

(i) in S 2 , 

V(p,<l>] <f>') = Wo(p,</»; <f)')+0{kp)~^] 

(ii) in /S 3 , 

V(py<l>y ^') = 0(kp)~~^\ 

(iii) in 8^, 
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These equations can be interpreted in a manner similar to that given 
above, and a closer examination of the order terms will give the 
diffraction effect. 


§ 2.2. The diffracted light 

The function T = cos which occurs in the wave- 

function vanishes on the radius vector (f> = <^'+7^. To 

avoid difficulties which this might cause, we consider only the part 
of physical space which lies outside the parabola = 1 ; this 

involves no loss of generality when we consider the effect when kp is 
large, since we are cutting out only a finite part of each radius vector 
except the critical one. In the parts of and S 2 which lie outside this 
parabola 1 

but in the part of outside the parabola 

T<--' 

€ VtT 

Again, since the function T' = cos — 37r)^(2Arp), which 

occurs in the wave-function u{p,j>\ Stt— vanishes on the radius 




direction of 
incident light 


direction of 


Screen 


vector <l> — <f>', we leave out of consideration the part of physical 

space within the parabola =1. In the part of 8^ which lies 

outside this parabola, we have 

m# ^ 1 


4617 


U 
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but in the parts of 8^ and 8^ outside the parabola, 

1 


T < 


E VtT 


We denote by /Sg, 8'^ the parts of 8^, 8^, and 8^ respectively which 
lie outside both parabolas. 

By using the results of §§ 1 .6 and 1 .7 we readily obtain the following 
approximate formulae: 

, p—ikpi-ikct-lTTi 

uM, i') = )_ (in s\ and SJ), 

^—ikp+ikc(—l7Ti 

U{p, <!>■<(,)=- (in ^s), 

fi—ikp+ikct—iTTi 

p—ikp+ikcl—lni 

u{p,<f}; Stt— ^') = ~ . -T/iT Jj\ ~r \ ^2 8^); 

2m\^{<f)+(f>')yJ{27Tkp) ^ 

the error in each case being less than ^Tre®. We now apply these 
approximations to the diffraction problems considered in the pre- 
ceding section. 

When the incident light is polarized perpendicular to the edge of 
the screen, the total resultant field is given by 

4 = df+d?, 

where d* represents the effect according to geometrical optics*]* and 

^-ikpJfikct-ini 

-{sec|(^->^')+cosec^(^+^')} (2.21) 




2y/(27Tkp) 

gives the diffracted light in the regions 8[, S^y 8^ with an error less 
than ^TTE®. To the same order, the magnetic vector in the diffracted 
light is found to be 

^-ikp+ikct-ini 

hs = -(isin^-jcos^)-^^-^— -{8ec|{^-,^')+cosecK^+f)}. 

When the incident light is polarized parallel to the edge of the 
screen the total resultant field is given by 

K = 

where h* represents the effect according to geometrical optics and 

£^~ikp+ikcl~^iTi 

27(-27rM (2.22) 

t By this, we mean that there is the ordinary incident and reflected light in 
only the incident light in and darkness in S^, 
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gives the diffracted light in the regions with an error less 

than To the same order the electric vector in the diffracted 

light is found to be 


— (isin^~j cos<^) 


^—ikp+ikri—lni 

'2yJ(27r/cp) 


{sec cosec 


The diffracted light rays are all perpendicular to the edge of the 
screen since, in both cases, the diffracted light diverges from the edge 
like a non-isotro]:)ic cylindrical wave. Thus, when the eye is focused 
on the edge of the scTcen, the edge should appear luminous although 
it is not a true source of light. This phenomenon can actually be 
observed in the geometrical shadow where it is not masked by the 
much greater intensity of the incident light. In the region /Sg the 
incident and diffracted light interfere and xiroduce the interference 
fringes observed experimentally near the boundary of the geometrical 
shadow. (Cf. Ch. 11, § 'h*!). 

The amplitude of the magnetic vector in the diffracted light is 
given by the aj^proximatc formulae 




1 

2V(27rA-p) 


(sec i^')+cosec^-(<)i+<^')} 


or A.. = ^2 i^lnk-p ) cosec 

according as the incident light is polarized })erpendicular or parallel 
to the edge of the screen. The ratio of these amplitudes is 

4 ^ = — COt(|7T-|-i^')COt(-j774-i^). 

^ L 

This should bo compared with the equation — obtained in 
the theory of the black screen (Chai)ter 111, § 4.2). 

The angle between the diffracted ray at a point F and an incident 
ray ])roduced through F is called the migle of diffraction, the lAsme 
of these rays the of diffraction. In the present case the angle 

of diffraction at (p,^) is 

8 = (f>~-(f)' 7T, 

the plane of diffraction being the plane of incidence. With this notation, 

A == cot(l7r+-Jf)tana,r+J84-Jf ). (‘2-23) 

We notice at once one imxDortant difference between this theory 
and the theory of the black screen; for the diffraction effect at a 

4617 XJ 2 
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point P due to a black half-plane depends only on the angle of 
diffraction and not on the angle of incidence, whereas with a reflect- 
ing screen it depends on both. 

Although the diffracted light is not plane-polarized, it does behave 
like plane-polarized light at large distances, and so we may still call 
the plane which contains the diffracted ray and the magnetic vector 
the plane of polarization. 

The plane of polarization of the diffracted light is parallel to that 
of the incident light in the two cases we have considered so far. We 
now consider the more general case of plane-polarized incident light 
in which the plane of polarization makes an angle a with the plane 
of incidence 2 : ~ 0. The incident field is then 


d = (isin^'sina— j cos^'sina-|-kcosa)e^*^‘/*‘**''^<'^-‘^'^+^‘^'^^, 
h = ( — i sin cos a-f-j cos cos a+k sin a)e^^P cos(t^~^') vikct^ 
The resultant field| turns out to be 

d = d*-fd^ h = h*+h». 


where d*, h* is the field of geometrical optics, and 


= — {AjLisin<^cosa:— 'A^jcos<^cosf>:-t-A;|ksinQ:}e“’^v> 

is the approximate magnetic force in the diffracted liglit. 
The amplitude of is 


] ik('t-ini 

(2.24) 


A = cos^a-f-A^sin^a} 


= A 


^|cOs2(X + 


tan 2 ( J 77 + 18 + 1(f)') . 2 

tan^ilTT+W) ‘ ‘ 


1- 


The intensity of the diffracted light is proportional to and so 
varies with a. Moreover, since 0 < 8 < 1^—^' in the geometrical 
shadow, the intensity is a maximum when <x ™ ^tt, that is, when the 
incident light is polarized parallel to the edge of the screen. This 
effect does not occur with a black screen since A , ~ A in that case. 

The plane of polarization of the diffracted light (2.24) at the point 
{p,<f>,z') is 


— Aia;sin^sina-f-A||ycos^sinQ:-fA^(2;— 2;')cosa' 0. 


This plane makes an angle 

tan“^ 



t This is obtained by supposing that the incident light is the result of superposing 
a field with a — 0 on a field with a = ^tt. 
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with the plane of diffraction z ~ 2', so that the plane of polarization 
is rotated in diffraction by a perfectly reflecting screen. | Since 

-^11 tan^ (| 7 r -l-S+ 0 O 

A tan 1(^77+^') ’ 

the direction of rotation depends on the angle 8. In the geometrical 
shadow A^^jA j > 1 since 0 < 8 < and so the plane of 

polarization is more nearly perpendicular to the plane of diffraction 
in the diffracted than in the incident light. J 

The [)henomcnon of the rotation of the plane of polarization in 
diffraction was first discovered theoretically by Stokes, jj who found 
experimentally a rotation in the sense opposite to that which we have 
just found — a result of considerable interesttf from the point of view 
of the elastic-solid theory of light with which Stokes worked. The 
experiment was repeated by Holtzmann, L. Lorenz, and Quincke 
with very divergent results. All these experiments were conducted 
with a grating, and, as Poincare i)ointed out, the effect looked for 
may be masked by a rotation due to refraction. The experiments 
of Gouy with a very acute-angled reflecting wedge give a rotation 
in the sense which the present theory indicates. The experiments 
of Wien'i'ill with very sharp steel blades require a more rapid increase 
of ^j| '.Ayy^ with 8 than is given by ( 2 . 23 ), which may be due to the 
finite thickness of the blade or to its not being a perfect conductor. 

§ 3. Diffraction of plane -polarized light by a black half-plane 

§3,1. Voigt's theory 

When ])laiie-polarized light specified by 

is incident on a reflecting screen a’ 0 , y < 0 , the effect is specified, 
as we have seen, by 

d = k{tt{p,<f>-, <{>')— u(p,<l>] 37t— 

t Thoro is no rotation witli a black screen, ns is confirmed by experiments of 
(jiouy. Ann. <fc chint. vt de phys. (6), 8 (1886), 14.5. 

i if tbe incident light is natural light with all possible directions of polarization, 
it is easily soon that the diffracted light is partially polarized. 

II In his paper on tlie ‘Dynamical Theory of Diffraction’, TranJt. Camb. Phil. Soc. 
9 (1849), 1, reprinted in Stokes's Papers^ 2 , 243. 

tt See, for example, Whittaker, Hiatory of the Theories of Aether arid Electricityf 
168. 

Theorie math, de la Lumidre, 2 (1892), 213—26. 

1111 Wied. Ann. 28 (1886), 117 . 
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In this equation the first term represents the effect of the incident 
light, the second the effect of the reflected light. Voigtf suggested 
that, if the screen is black, the effect should be obtained by omitting 
the term due to the reflected light, so that 

d (3.11) 

Similarly, when the incident light is polarized parallel to the edge 
of the black screen, . . / //\ 

There is no theoretical reason for making this assumption, since an 
opaque non -reflecting screen must absorb all the incident light, which 
is impossible if it is very thin. It is, however, interesting to work out 
the consequences of the assumxflion, es])ecially as it gives results 
which differ from those obtained by using Kottler’s modified Larmor- 
Tedone formula. 

In the first place, it suffices to consider only the case when the 
incident light is polarized perpcndicndar to the edge of the screen, 
since the results when tlie plane of polarization is parallel to the edge 
can be obtained by the transformation 

d -> h, h -> — d. 


It follows immediately that there will be no rotation of the plane of 
polarization, just as in Kottler’s theory. 

The field specified by (3.11) can be written in the form 

d - d*+d^, h - h*+h^^ 

where (d*,h*) represents the field of geometrical optics and (d^^h"^^) 
the field of the diffracted light. In the f)art of ])hysical space 
(““ Itt < <56 < 377 ) outside the parabola 


we have 




2kpQOH^ i {(f) -—(/>') ~ 1/(7T€“), 
p ~ikp } ikci-\ni 

2 COS l(<f>~<t>')sJ{27Tkpy 

p kct~ {Tri 


~ —(isin^— j cos<i) 

2 cos \ {<f>~(f>')\j{2TTkp) 

with an error less than \7re^. The diffracted light, as before, diverges 
like a non-isotropic cylindrical wave from the edge of the screen; the 
apparent brightness of the edge again agrees with experiment. 

There is, however, an important difference from the theory of the 
reflecting screen in that the diffracted light depends only on the angle 


t Odtt. Nach. (1899), 1. 



§3] 


SOMMERFELD’S THEORY OF DIFFRACTION 


151 


of diffraction 8, so that dilfraction is here purely an edge effect. The 
amplitude of the diffracted light is 

^ 2sin |8^(277i:p)‘ 

This should be compared with the amplitude 

^ 2 tun lS^J(2Trkp) (3.13) 

given by Kottler’s formulae or by Kirchhoff\s scalar theory in its 
exact form. 

Actually, (3.13) agrees with experiment better than (3.12), so that 
Voigt’s black screen is not black enough. It has been suggested*)" 
that matters could be improved by increasing the number of sheets 
p of Sommerfeld’s Riemann surface; for then light would be less 
likely to enter physical space from non-physical space. In particular 
by making -> oo, we obtain V^oigt’s ‘blackest’ screen. The corre- 
s])onding solution is then 


^ikct 

27Ti 




gj'A'prosf 




dC 


(3.14) 


where A is the two-branched path of Fig. ‘23 on 136. 

By u.sing t'auchy’s theorem it can be shown that, in this field, 

d = d*-fd«, 

where d* is the electric force according to the laws of geometrical 
optics and 


d« 


— — 2ke’^'''' J* 


g-iArpeoshr 


(7T“ 


— 02 i ^2 


d"= -k 


(where d = 0—^') is the electric force in the diffracted light. 

The approximate formula 

2^g-iT*p+iA'd-j7ri 

{7r^-d^)^(27rkpr 

valid when p is large compared with the wave-length 27 t/A’, can be 
deduced from (3.15) by the principle of .stationary phase. Hence the 
amplitude of the diffracted light with Voigt’s ‘blackest’ screen is 

^ (3.16) 


A 


when p is large. 

t See, for example, Sonimorfeld, Zeits, J. Math. u. Phys. 46 (1901), 11-97, §5 
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Formulae (3.12), (3.13), and (3.16) all give the same approximate 
amplitude l/{S-y/(27rifc/o)}, when S, the angle of diffraction, is small, but 
do not agree when the angle of diffraction is large. 

§ 3,2. The unsatisfactory state of the theory of the black screen 

Sommerfeld’s theory of diffraction by a perfectly reflecting screen 
can be regarded as a satisfactory one. The boundary conditions 
satisfied by his solution are easy to formulate on theoretical grounds, 
and the structure of his ‘Rieniann space’ is fixed by the conditions 
of the problem. On the other hand, although each of the theories of 
diffraction by a black screen gives results which agree fairly well with 
the available experimental evidence, not one of them rests on a sound 
theoretical basis. 

Voigt’s work suffers from two disadvantages: the electric and 
magnetic vectors satisfy no prescribed boundary conditions and the 
number of sheets in the ‘Riemann space’ is quite arbitrary. Kottler 
tried to overcome these difficulties by regarding the problem as a 
saltus problem instead of a boundary -value problem. After a critical 
study of Kottler ’s work, Ignatowskyf has concluded that there is no 
definition of a black screen applicable in all cases, and he asserts 
that each black-screen diffraction problem must be discussed inde- 
pendently on its own merits. 

The real difficulty lies in the fact that the black screen is an 
idealization which cannot be attained experimentally and which has 
no precise definition in electromagnetic theory. What is needed is 
a rigorous theory which does not assume that the screen is perfectly 
reflecting but takes into account the properties of the material of 
the screen. An empirical way of finding the effect due to a screen 
which is neither perfectly black nor perfectly reflecting is to suppose 
that the reflected wave has its amplitude and phase reduced in some 
definite way, due to the imperfections of the reflecting power of the 
screen. This amounts .to multiplying the term representing the 
reflected wave in Sommerfeld’s solution by a complex constant. 
The effect of this modification of Sommerfeld’s theory has })een 
worked out in some detail by Raman and Krishnan. J 

t Annalen d. Vkys. 77 (1925), 589-64;{. 
t Proc. R.S. (A), 116 (1927), 2.54. 
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